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1 Introduction 



Since the discovery of quantum black hole thermal radiance by Hawking [1 (extremely 
good reviews and surveys are in [2], [3J, [3], [5] and [5]), it became pretty clear that 
something remarkable concerning the interface of gravity, quantum theory and ther- 
modynamics was at work. In the usual picture, a radiating black hole loses energy 
and therefore shrinks, evaporating away to a fate which is still debated. Many new 
ideas came out from the recognition that quantum field theory implied a thermal spec- 
trum: the most impressive probably being 't Hooft's idea of a dimensional reduction 
in quantum gravity and the associated holographic description (indicating a drastic 
reduction in the counting of degrees of freedom in finite regions, which scale with area 
rather than volume) [7J[S]; and the principle of black hole complementarity aimed to 
reconcile the apparent loss of unitarity implied by the Hawking process with the rest 
of physics as seen by external observers [9]. But there were also other, more practi- 
cal, issues regarding these matters, some of which bewildered scientists since the very 
beginning and that have been only partly resolved. A key issue is that the original 
derivation of Hawking radiation applied only to stationary black holes, while the pic- 
ture above uses quasi-stationary arguments. Actually, an evaporating black hole is 
non-stationary. Nevertheless, a surprising aspect of the semi-classical result is that 
the radiation caused by the changing metric of the collapsing star approaches a steady 
outgoing flux at large times, implying a drastic violation of energy conservation. This 
certainly means that one cannot neglect the back-reaction problem which, however, 
has not been solved yet in a satisfactory way. A sample of other key issues have to 
deal with the final state of the evaporation process, the thermal nature and the re- 
lated information loss paradox, the Bekenstein-Hawking entropy and the associated 
micro-states counting, the trans-planckian problem, and so on. It was to address some 
of these questions that some alternative derivations and descriptions of the Hawking's 
emission process began to appear over the years [TTTHT8] . one of these being the so 
called tunnelling method, the leitmotiv of this review. 

Before introducing it in full mathematical detail, we shall spend the next few para- 
graphs to spell out some general features of this method and how it compares with the 
general properties of quantum fields and states in the presence of a black holfl 

The tunnelling method has to do with the particle interpretation of the emission 
process, which would be the most natural way to explain the loss of energy suffered 
by the radiating black hole, but also poses a problem. It is that while the original 
derivation of Hawking radiation used the theory of quantum free fields propagating in 
a fixed gravitational background, for the case at hand a particle interpretation in the 
near horizon region of a black hole space-time was found to be far from unique and 
therefore inherently ambiguous. This fact is deeply rooted into the physical nature 
of the problem, which is that the general covariance of the physical laws allows many 
inequivalent choices of time, as was explained several times (see, e.g. [^HHH]). It was 
partly to overcome such difficulties that DeWitt launched in the Seventies the "stress- 

For those interested, an up to date resource for papers, books and other stuff about black hole 
physics and astrophysics is the "Resource Letter BH-2: Black holes" by Marolf and Gallo |19| . 
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energy-tensor program" aimed to compute the expectation values of the stress tensor 
and other observables to better describe quantum effects in a curved space-time [22] . 
Thus, it was found that it is not in general possible to divide (I) t!y ) re „ into a real particle 
contribution and a vacuum polarisation part in an unambiguous way. However, it will 
always be possible to choose any coordinate system, as long as it is regular across 
the horizon, and use it to define an observer dependent vacuum relative to which a 
particle definition is feasible. As we will see, with the tunnelling method we are only 
concerned with the probability that such an observer dependent notion of particle be 
emitted from the horizon of the black hole. If this probability is a coordinate scalar 
it will not depend strongly on what particle concept one employs. In fact, we can call 
it a click event in a particle detector if we like, without committing ourselves with 
the particle concept; concept which remains of great heuristic value however, and will 
find its realisation in the choice of the particle action used to compute the tunnelling 
probability. Such an observer dependent notion of particles has been advocated also 
by Gibbons and Hawking in their treatment of de Sitter radiation |23) , on the ground 
that an observer independent definition of particles is not relevant to what a given 
observer would measure with a particle detector. 

The second aspect we would like to mention is the fact that, as it will be seen, the 
tunnelling method uses only the classical action of a single massless, spin-less particle 
and therefore appears to be state independent. It does what it does by relating the 
particle emission to an imaginary contribution to the classical action localised at the 
horizon, which only comes from the local geometry. This can be traced back to the 
work of Damour and Ruffini [24] with contributions of Padmanabhan et al. [25lI28j . 
Massar and Parentani 29j and Kraus and Wilczek, Parikh and Wilczek [3UH12] (see 
also Kraus's Ph.D. thesis [53]). To what extent this is so was recently studied in [M] 
where, for a large class of quantum states with Hadamard short distance singularity, it 
was shown that the two-point correlation function, which is directly related to the tun- 
nelling probability, has a universal thermal appearance whose temperature is, needless 
to say, the Hawking temperature. 

Now all these important findings have the drawback that they have been proved, 
and one can easily understand why, only for stationary black holes. In that case the 
tunnelling picture has not much more to say than was already known, at least for the 
simplest solutions, besides helping to understand black hole radiation from a different 
viewpoint. Where exact calculations can be done, it typically gives a less detailed 
picture of the radiation process since it is mainly related to a semi-classical emission 
rate. The cases where the method is more powerful have to do either with more 
intricate stationary solutions or else with truly dynamical black holes. We can list 
here some of its achievements and properties: 

1. The estimation of the leading correction to the semi-classical emission rate aris- 
ing from back-reaction to the background geometry, whose introduction by Kraus 
and Wilczek, Parikh and Wilczek 30-32 motivated the tunnelling method in a 
form that will be discussed soon; its extension to Anti-de Sitter (AdS) [35], de 
Sitter (dS) [5BH4D] and higher dimensional black holes; the existence of correla- 
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tions among successive emission events [4lTl43| : 



2. The original tunnelling method can be generalised to a Hamilton-Jacobi variant, 
originated with the work of Padmanabhan and collaborators [2M28] and system- 
atically applied either to stationary or dynamical black holes (see e.g. [12,44 49 
for a sample of papers) . For dynamical black holes this was particularly impor- 
tant, since even approximate quantum calculations are notoriously hard. More- 
over, in this variant it manifestly preserves general covariance; 

3. Supplemented with a more precise and more general mathematical definition 
of a local horizon, the Hamilton-Jacobi variant can be applied to any sort of 
horizon within this class, and in particular to cosmological and weakly isolated 
horizons |50j (for definitions see below). It can also be applied to past horizons 
and white holes, in which case a clear notion of temperature emerges in complete 
analogy to black holes; 

4. The tunnelling method strengthens the connection of the semi-classical rate with 
the surface gravity of the horizon even for dynamical black holes, provided that 
opportune definitions of such quantities are employed. In this respect its appli- 
cation to extremal black holes should have something to say; 

5. The tunnelling picture can promptly give indications on what occurs in other 
situations. One example is the WKB approach to Unruh's radiation, reviewed 
in [5T]. It can also be applied to see whether naked singularities are going to 
emit radiation, or to study the decay of unstable particles that in the absence of 
gravity would be stable by the action of some conservation law. This is valuable 
given the great efforts that are normally necessary to analyse quantum effects in 
the presence of gravity. 

It remains to explain how the tunnelling picture works and to what extent it gives 
a consistent picture, at least from a logical point of view. This will cover the rest 
of this review whose structure, we anticipate, will essentially follow the material just 
indicated in listed form above. 

In Section 2, the methodological core of the topic is outlined and applied to various 
stationary solutions. In Subsection 12.11 we present the original method introduced 
by Parikh and Wilczek to study the black hole radiance — namely the null geodesic 
method — which is subsequently applied to verify the thermal properties of a static, 
spherically symmetric Schwarzschild black hole in the scalar sector. In Subsection 12.21 
we comment results pertaining to statistical correlations of successive emission events, 
as displayed by the tunnelling method. In Subsection 12.31 and Subsection 12.41 after 
pointing out some flaws of this approach, we extend it to the more recently conceived 
Hamilton-Jacobi tunnelling method: this approach is examined in great detail in order 
to stress those points that — given the vastness and heterogeneity of opinions in the 
literature — might lead to confusion. In Subsection l2.5l we present the extension of the 
Hamilton-Jacobi tunnelling method to the fermionic sector. The description of the 
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tunnelling phenomenon in non-static, but still stationary space-times is undertaken in 
rest of the section: starting from the Kerr and Kerr-Newman solutions in Subsection 
12.61 we then conclude the section showing in Subsection 12.71 more peculiar scenarios 
— such as non axis-symmetric or lower and higher-dimensional solutions. 

In Section [21 we study the mathematical foundation of the two tunnelling methods. 
The Hamilton- Jacobi method, preferable in view of its manifest covariance, is exam- 
ined in Subsection l3.ll It is found to be equivalent to the null geodesic method, at least 
in its generalised form, in Subsection l3.2l The aforementioned variety of opinions also 
comprehends a number of alternative proposals to the two methods, some of which 
are summarised in Subsection 13.31 

In Section^ we approach the dynamical case but restricting ourselves to spherically 
symmetric non vacuum black holes. Though it may appear strange that stationary 
black holes, being a particular case of dynamical ones, are discussed first, there is a 
reason: in the stationary case there is much more control on the mathematics of the 
tunnelling picture. In particular, it is possible to use both the null geodesic method as 
well as its Hamilton- Jacobi variant to deduce the exact semi-classical thermal spectrum 
of the Hawking radiation, a fact which greatly enhances the understanding of the 
tunnelling picture in a way which is not available to dynamical situations. At the same 
time it provides the foundation for the dynamical treatment and inspires extensions 
to more general cases. The Section is further divided into several sub-sections. In 
Subsection 14. 1 1 we describe various definitions that have been proposed over the years 
to characterise a black hole horizon and corresponding definitions of the surface gravity. 
Again few dynamically inequivalent definitions will be reviewed, putting particular 
emphasis on the one which the tunnelling method privileges. In Subsection 14.21 we 
review material pertaining to spherically symmetric space-times and the form that 
surface gravities and horizons take in this highly symmetric case. In Subsection l4.3l we 
finally come to grasp the tunnelling mechanism. After a careful scrutiny of the escaping 
paths we derive the basic results; the computation of the imaginary part of the action, 
its relation with the surface gravity and the ensuing quasi-thermal interpretation. 

In Section [SJ we discuss how the tunnelling picture generalises to cosmological hori- 
zons, naked singularities and particle decays. We shall see that the first two cases 
are governed by a simple pole singularity while decays correspond to branch cut sin- 
gularities of the classical action. In Subsection 15.11 we review few facts regarding 
Friedmann-Robertson-Walker space-times (abbr. FRW) with emphasis on the exis- 
tence of horizons, their structure and surface gravity. In Subsection 15.21 we apply the 
tunnelling method to cosmological horizons and find their emission rate. In Subsec- 
tions [531 and [52] we analyse its application to particle decays and particle creation by 
black hole singularities. The former was extensively studied in de Sitter space-time, 
where exact results can be obtained, but the tunnelling method is not restricted to 
particular space-time choices. In Subsection 15.51 we review the possibility that naked 
singularities themselves radiate away their (suitably defined) mass. Indications of this 
phenomenon are supported by the tunnelling method. 

We end with Section [5] where achievements, prospects and conclusions are sum- 
marised. 
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As a final note, since we found the literature on quantum tunnelling to be quite 
extensiv^], the references will not exhaust all the existing literature. We struggled to 
cover a lot of papers, and especially those that appeared more meaningful to our eyes, 
being aware that this will inevitably come to reflect some personal taste. A sample of 
the vast literature is here collected for the reader's benefit, beyond those already cited. 
For studies on tunnelling from a general point of view see |196) ; tunnelling of fermions 
and bosons from stationary black holes in asymptotically flat, de Sitter or anti-dc 
Sitter spaces and other exotic objects in several coordinate systems were considered 
in [197] ; tunnelling of bosons and fermions in higher and lower dimensions are covered 
in [198] ; for contributions to dynamical black holes see |199j , for cosmology pOO . Our 
sincere apologies are directed to the many unwittingly omitted. 

We shall use the following conventions throughout: the metric signature is (—,+,■■ ■ , +); 
the constants c, G, ks and h set equal to unity unless explicitly stated; we use an in- 
dex notation for tensor components. Greek indices are used to denote space-time 
components of a tensor; mid-Latin indices as i, j are typically used to denote the 0, 1 
component; late-Latin indices as m,n for purely spatial tensors part. If not otherwise 
stated, an overdot will represent a time derivative, and a prime a spatial derivative in 
the (case by case) relevant spatial coordinate. 

2 Quantum Tunnelling: Stationary Black Holes 

According to quantum field theory, the physical vacuum - the state with no real 
particles - is in general a complex entity populated by virtual particles which constantly 
are created, interact among themselves and are then annihilated in this vacuum. In 
the absence of external fields, the vacuum is usually stable in the sense that virtual 
particles are not able to survive long enough to become real. However, it is now well 
proved that external fields are able to convert short-living particles into real ones just 
supplying enough energy to the process. 

Considering particle creation by a static gravitational field endowed with a Killing 
vector £ a , the energy of a particle is w = — p Q £ Q , where p a is the particle's four- 
momentum. For massive (massless) particles, the momentum p a is a future-directed 
time-like (null) vector. Thus, the energy a; of a particle is always positive in the region 
where the Killing vector is also future-directed and time-like. It follows that particle 
creation in a space-time region with £ Q future-directed and time-like is forbidden by 
energy conservation. Such considerations provide a heuristic argument [5 2) to conclude 
that a static gravitational field can create particles only if the space-time contains a 
black hole. In fact, when a virtual particle pair is created just inside the horizon, the 
positive energy virtual particle can tunnel outside — no classical escape route existing 
— where it materialises as a real particle; in alternative, for a pair created just outside 
the horizon, the negative energy virtual particle, forbidden outside the horizon, can 
tunnel inwards |53] . In either cases, the negative energy particle is absorbed by the 
black hole, resulting in a black hole mass decrease, while the positive energy particle 
escapes to infinity, visible to distant observers (Hawking radiation). 

2 An INSPIRE search gave us more than 500 records about tunnelling from black holes, to date. 
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In 1974 Hawking [I] proved that indeed black holes are classically stable objects 
that, due to quantum fluctuations, become unstable with respect to particle emission. 
Nowadays, several methods are known for deriving Hawking radiation and calculating 
its temperature, but none of them has been able to catch the intuitive picture above as 
the tunnelling method proposed a decade ago by Parikh and Wilczek [32] : it were not 
for the inclusion of back reaction, their method could be considered a semi-classical 
version of Hartle-Hawking's path integral derivation of black hole radiance [54] . 



2.1 The null geodesic method 

The so called null geodesic method — as named by Kerner and Mann |46j — derives 
black hole radiance from few reasonable assumptions, namely: energy conservation, 
detailed balance and the use of non-singular coordinates across the horizon. 
Energy conservation requires fixing the total energy of the space-time before and after 
particle emission. Since black hole mass and volume are linked together, a mass 
reduction due to the emission of a particle translates into a size contraction; so one 
might worry how to deal with quantum fluctuations of the metric originating from such 
contraction. However, this is not a problem as far as we consider transitions between 
zero-spin geometries. In this case in fact, no graviton quantisation is involved or, said 
in other words, passing from different spherically symmetric configurations does not 
produce gravitational waves. As a consequence, the only degree of freedom remained 
in the problem is the position of the emitted particle (actually, a thin shell). Thus, 
to keep things as simple as possible, we can restrict to consider uncharged, static, 
spherically symmetric black holes emitting neutral matter, referring to next sections 
treatment of the most general case. 

Because of Birkhoff's theorem [55], we are dealing with Schwarzschild geometry 
which, written in the Schwarzschild frame, reads 

ds 2 = - ( 1 - T -A dt 2 + ( 1 - T -A ~* dr 2 + r 2 {d9 2 + sin 2 0d<f> 2 ) (2.1) 



r j v r 



where the gravitational radius is r g — 2M . Computing the free fall acceleration of a 
body initially at rest in the Schwarzschild frame, we get 

M 



a:=^a~a^= (2.2) 
Hyl - r g /r 

The geometry (|2.1I) becomes singular at the gravitational radius and the acceleration 
(|2.2p becomes infinite in the same limit. Nonetheless, near r g the space-time is regular 
(i.e. neither curvature singularities or geodesic incompleteness): an indication of the 
fact that the Schwarzschild frame ceases to be defined at r — r g (For r < r g systems 
must be non-rigid in the sense that the metric components g a p must be functions of 
time, cf. [51]). 

It is well-known that the gravitational radius r g represents the Killing event horizon 
of Schwarzschild black hole [55]; thus, in order to deal with tunnelling of a particle 
across it, we need to cover the region r ~ r g with regular coordinates. A possible 
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choice, not the only one, is represented by Painleve-Gullstrand coordinates [56]. In 
order to construct this new regular-across-the-horizon frame, we consider the geodesic 
equations of a family of observers, freely falling (dr/dr < 0) in the geometry (|2.ip . 
along a radial trajectory (9, <f> constant): 



dt E I dr 



l-!±\=E 2 (2.3) 



dr 1 - \dr J V r 

r representing observer's proper time and E = E/m is the observer's conserved energy 
per unit mass. The energy parameter is related to observer's physical three-velocity 
v 2 = g m ndx m dx n I '(—g o)dt 2 : in fact, given the Killing vector — generating the 
time-symmetry transformations of (|2.1I) . the energy reads 



r- , /' V 500 ,.11, 

E = -u M C( t) = -mo = -5oo u = = . (2.4) 



We specialise to the case of observers starting at infinity with null initial velocity so 
that E = 1. For this family of observers, there exists a function t p such that 

V^p = —Up . (2.5) 

Taking (|2.3[) into consideration, we get the Painleve time t p in terms of Schwarzschild 
coordinates 

^ " ' + f ^ (»-?)"'-« + W + r. 1" (£=$f) . ,2,6) 
In terms of (t p , r, 0, 4>), the Schwarzschild geometry is written as 

ds 2 = -(l-^dtl + 2 ^dt p dr + dr 2 + r 2 (d9 2 + sin 2 9d4> 2 ). (2.7) 

At fixed time, the spatial geometry described by (|2.7I) is Euclidean; while at any fixed 
radius, the boundary geometry is the same as that in (I2.1j) . Most importantly, the 
metric is no more singular at the gravitational radius r g ; the space-time is stationary, 
but no more static. The t p coordinate — being the time perceived by infalling observers 
who reach the curvature singularity r = in finite time — remains a valid coordinate 
as far as it remains finite. 

The geodesic for a null s-wave is given by 

with upper (lower) sign corresponding to outgoing (ingoing) geodesic. The basic idea 
behind the approach by Parikh and Wilczek |32] is to study the emission of such a null- 
radial thin shell from the black hole through quantum tunnelling. Imposing energy 
conservation means that the total space-time energy is fixed and one allows the black 
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hole mass to fluctuate. As showed in [30], the motion of a shell of energy w emitted 
from a black hole of initial mass M is described by ea. (|2.8p in the space-time geometry 
of eq.(H21with M replaced by M - uj. 

Differently by tunnelling in quantum mechanics, where two configurations (classi- 
cal turning points) separated by a potential barrier are connected by a trajectory in 
imaginary time, here we have the perception that nothing similar happens. In fact, 
the problem with black hole radiation is that, as soon as a particle is displaced outside 
the horizon, it can escape along a classical trajectory. It is worth to mention how the 
crucial ingredient of the story here is played by energy conservation. As the black 
hole of mass M and area WwM 2 emits a particle of energy ui, its mass reduces to 
M — oj and its horizon recedes to 2{M — oj) — r g — 2u>. This amount of contraction, 
evidently, depends upon the mass of the emitted particle, so that quoting Parikh "it 
is the tunnelling particle itself that secretly defines the tunnelling barrier" 53J . 

Before evaluating the black hole emission rate, we can ask to what extent is justified 
the point-particle approximation. Given that a distant observer detects a wave with 
frequency fl b, this has been emitted with frequency fl em oc (1 — r g /r)~ 1 t 2 VL b. In the 
limit of r — > r g , where the emission process occurs, the wavelength vanishes, making 
the point-particle (WKB) approximation fully reliable. 

In the WKB approximation, the tunnelling probability rate T em turns out to be 
equal to 

r em ~exp(-2 Im/) (2.9) 

where / is the action for the tunnelling process evaluated to the appropriate order. 
We shall derive later, in Subsection 13.11 the exact correspondence between T em and 
the exp-function. Eventually, one would expect to be able to read off the black hole 
radiation temperature from comparison of the the probability rate eq. (|2.9p with the 
Boltzmann factor e~"/ T . The imaginary part of the action for a null s-wave outgoing 
positive energy particle which crosses the horizon outwards from r, n to r out can be 
expressed as 

Im7 = Im / drp r = Im / dr / dp' r (2.10) 



Changing the integration from momentum to energy thanks to Hamilton's equation 
r = dH/dp r and noting that H = M — uj' with constant M and < oj' < to, we have 

Im I = Im f ^ dr — = Im f °" dr f . (2 .H) 

Jv in JM r Jr in Jo r 

Using ea. (|2.8[) and switching the order of integration, 

rTout i 

Im I = Im / f-dw') / dr ; . (2.12) 

JO Jr in l _ . / 2(M-u') 



With the understanding that the particle starts at n n = 2M and materialises at r out = 
2(M — uj), ri n > r out , the integral over r is divergent and needs to be appropriately 
regularised. A tunnelling process corresponds to an energy which is not in the spectrum 
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of the Hamiltonian. We need to continue the energy lo to complex values, that is 
lo' — > lo' + ie, in order to ensure that positive energy solutions decay in time: 



Tout 



Im / = Im J (—dco') J dr 
= — Im in 



r- y/2(M - (7) + ie 
I (-duo 1 ) 4(M - J) = +4ttlj (M - ^) . (2.13) 

Of course, Hawking radiation can also be regarded as pair creation outside the horizon, 
with the negative energy particle tunnelling into the black hole. The calculation pro- 
ceeds as above with slight change of signs due to the fact that anti-particles travel back 
in time (cf. |32j for details). That both particle and anti-particle channels contribute 
to black hole emission is something which only affects the amplitude of the process, 
that is something which enters the proportionality factor of eq. (l2.9[) . In conclusion, 
the emission rate obtained is 

r ~ exp(-2Im/) = e - 87rM "(^w). (2.14) 

The expected Boltzmann factor e _ "/ Tff is recovered only to order 0(d): in this or- 
der of approximation, the black hole temperature perfectly coincides with standard 
Hawking's result, Tjj = (87rM) _1 . The 0(lo 2 ) correction arising from the physics of 
energy conservation, makes higher the effective temperature of the hole as it radiates 
in agreement with the well-known negative thermal capacity of Schwarzschild black 
hole. 

Neglecting the lo(2M)~ 1 term in (|2.14[) and invoking detailed balanc^], the Planck 
distribution is recovered. 

A note - It appears that by neglecting the back-reaction no member of a pair of par- 
ticles created just inside the event horizon can escape it, by reason of causality. In the 



^Consider an ensemble of many identical copies of the same quantum system. Let the energy and 
the number of accessible states of the system be fixed. The probability that a system randomly 
chosen out of the ensemble will be in state i is denoted by pi, with J^V pi = 1. The transition 
probability Ai_>j denotes the conditional probability per unit time of a system going from state i 
to state j, i.e. piXi—,jdt is the probability of a system originally being in state i and going to state 
j within the time interval dt |57| . Prom this definition, • piXi^jdt represents the probability 
of transition from state i to any other possible state of the system in time dt; and Y2j Pj^j—tidt 
represents the probability of transition into state i from any other possible state of the system in 
time dt. It follows that, pi = ^2j(pjXj^i — Pi^i^j) is the net rate of change of pi. The detailed 
balance condition requires that, 

at equilibrium, pjXj—^i = piXi^j , Vi,J . (*) 

In the case at hand, the system is made of "black hole + radiation" . The initial state i represents, 
for example, the black hole with TV — 1 particles (e.g. photons). Through the physical mechanism 
described before, the black hole emits one more particle, so that the state j becomes "black hole 
+ N particles" . Neglecting back-reaction, the internal degrees of freedom of the black hole do 
not change and by (*), we have that Xi^j/Xj^i = r em /T a 5 = pjv/pjv-l = e~^l T . This result 
combined with the classical constraint, r„5 ^ T em = |T;(o;)| 2 , reproduces Planck (minus sign) and 
Fermi— Dirac (plus sign) distributions. Here, T[(uj) represents the transmission coefficient of the 
black hole barrier which in general can depend on the energy u and the angular momentum I of 
the particle. 
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terminology to be used in Subsection [4] this correspond to a type-II tunnelling path, 
the particles created just outside forming instead a type-I path. The inclusion of the 
back-reaction is truly fundamental for the existence of both possibilities. Since in the 
dynamical case we will mainly ignore the back-reaction effect, only type-I path will be 
accessible to the tunnelling probability. 



2.2 Entropy and correlations 

The exact expression for the emission rate ea. (|2.14|) can be intriguingly re- written as 
the change in the hole's Bekenstein-Hawking entropy [58] . 

T ~ e -^ M ^Tt) = e AS BH _ (2.15) 

This is an interesting form since it agrees with Einstein's formula for the probability of 
fluctuations about statistical equilibrium 59 . In this case, it refers to the probability 
rate with which a static Schwarzschild black hole gets out of equilibrium by emission 
of a scalar quantum. This formula is derived to order 0(u>) from the fact that [60] 

r em ~e- 2 Im/ =e-^=exp ^ex P -[S(M)-S(M-uj) + 0(u 2 )] , 

(2-16) 

where the standard relation involving the inverse temperature at equilibrium, ft = 
OS/dM, has been used. It is remarkable that Equation (|2.14l) puts forward this result 
to order 0(uj 2 ). Since the exp- function in Equations (I2.15|) and (|2.16[) can be re- written 
as expASsH — exp Sfi na i/ exp Si n u a u we can re-interpret it as a factor counting the 
internal degrees of freedom of the black hole. It is a long debated question how the 
thermal nature of Hawking radiation can be reconciled with unitarity {information loss 
puzzle). From one side, unitarity is a milestone of classical and quantum physics; from 
the other, Hawking radiation is the most reliable result of quantum gravity derived 
with semi-classical techniques. This is not the place to enter the details of such a 
wide subject: it is enough to say that the discovery of non-thermal corrections to 
the black hole radiation spectrum due to the inclusion of back-reaction effects showed 
by Parikh and Wilczek, and the possibility of writing the emission rate in the form 
()2.15[) gave hope for a possible solution of the problem. For example, let us take into 
consideration the emission of a quantum of energy u> = uj\ + ui2 and, separately, the 
emission of a quantum of energy uj\ followed by the emission of a quantum of energy 
ll>2- The function 

x^ = l °z\r(M^\ (2 ' 17) 

can be considered as a measure of the statistical correlation between the two events 
which is expected to be zero whenever the events are independent. This is indeed 
the case for thermal spectrum, but also for the emission rate (12. 14)) it was argued 
in [21] to be the case. Even the inclusion in Equations (|2.15l) . (|2.16l) and (|2.17[) of 
logarithmic corrections to the black hole entropy, as suggested by direct count of black 
hole microstates both in string theory and loop quantum gravity, seemed not to be able 
to show the appearance of correlations in the spectrum of black hole radiation [62] . 
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However, a more careful statistical treatment due to Zhang et al. in two significant 
papers [41,42 showed that, on the contrary, the statistical correlation ea. (|2.17[) is 
actually non zero and equal to 8itujiU!2, a result later confirmed by Israel and Yun |43| . 
To make clear the point, we display the black hole mass into r(w) and, following Israel 
and Yun, we write the probability that a black hole of mass M emits a quantum of 
energy cj as 

r(M, u) = exp[-47rAf 2 + 4tt(M - uj) 2 ] (2.18) 
which is the tunnelling result. Then this yields the non trivial correlation 

i \ 1 r(M,wi + w 2 ) q , oin -. 

x(wi,w 2 ) = log -— — — r =8tto;iW2 (2.19 

T{M,uji)T(M, u> 2 ) 

which is the result of Zhang et at. In general, one may define a measure of the 
correlation between events x and y as 

x(, I y)=log^^ T =log^ (2.20) 
P{x)P{y) P{y) 

where 

P{x) 

is the conditional probability of y given x and P(x, y) is the probability of both x and 
y. Clearly, if we substitute P{y\x) in place of P(y) into x{ x i y)i we § e t zero correlation 
for any events x and y, because in this way one absorbs the correlation itself into the 
test for its existence. Thus, by the same token, if we put T(M — w\, U2) into xi^i, ^2) 
in place of T(M, u)%), as it might be thought in order to take into account the mass 
loss accompanying the first emission, we find no correlation at all, as can be directly 
verified. 

We see that the tunnelling method leads to a very important conclusion, that the 
radiation does indeed contain correlations that possibly can carry off the information 
content of the hole. 



2.3 The Hamilton-Jacobi method 

Despite the merits of the seminal work by Parikh and Wilczek, we cannot omit to point 
out a couple of unpleasant features of their null geodesic method, as the fact that: (i) 
it strongly relies on a very specific choice of (regular-across-the-horizon) coordinates; 
and (ii) it turns upside down the relationship between Hawking radiation and back- 
reaction. As regard the former point, it should be clear how irrelevant is, in the spirit of 
general relativity, the choice of coordinates: being physical observables invariant with 
respect to the group of diffeomorphisms (the hole temperature is such an observable) , 
there is no reason why Painleve-Gullstrand coordinates should be favourable with 
respect to other (equally well-behaved) coordinates. With respect to the latter, we 
notice that, in the null geodesic description, apparently there cannot be Hawking 
radiation without back-reaction: watching carefully, however, it is the discovery of 
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Hawking radiation that justifies back-reaction and makes commendable the treatment 
of Hawking radiation's self-gravity. The so-called Hamilton- J acobi method can cope 
with both issues. The intent is to give a particle description of Hawking radiation, 
ignoring its self-gravitation, under the assumption that the emitted (scalar) particle's 
action does satisfy the Lorentzian Hamilton- Jacobi equation. Later we will show that 
the null geodesic method can do the same job using instead the reduced action. 
As it will become clear later, this method applies to any well-behaved coordinate 
system across the horizon; it generalises beyond the assumption of spherical symmetry; 
it makes possible to include the study of tunnelling by fermionic particles. Some 
sceptics about Hawking radiation contest the fact that as soon as the black hole starts 
radiating, the same assumptions on space-time stationarity drops down, invalidating 
the whole derivation. Given that the departure from perfect stationarity is ridiculously 
small, nonetheless, as we shall show in Subsection |4l the Hamilton-Jacobi method can 
prove black hole evaporation even for slowly varying, time dependent, space-times. 
Finally, we wish to notice that since the methods of tunnelling are intimately related 
to the physics of (some type of) horizons, we may apply them even to space-times with 
multiple horizons. In standard computations of Hawking radiation, this is typically a 
hard task, if not impossible. In fact, it is well known that, as an example, Reissner- 
Nordstrom-de Sitter space does not possess regular Euclidean section for general values 
of mass, electric charge and cosmological constant |63] . This means that for arbitrary 
values of the parameters, it is not possible to compute the Hawking temperature of 
event or cosmological horizons by Euclidean continuation, simply because it could be 
that no Euclidean section of the given Lorentzian space-time exists. 




Figure 1: Carter-Penrose diagram of the eternal Schwarzschild black hole. H de- 
note future and past horizon; I ± future and past null infinity; r — the 
curvature singularity; v,u represent advanced and retarded null coordinates, 
respectively, abc is a null piecewise continuous path from inside to outside 
the hole, with ~a% running backward in time (classically forbidden trajectory). 

In its minimal formulation, the Hamilton-Jacobi method works according to the 
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following strategy [64]: (a) assume that the tunnelling particle's action satisfies the 
relativistic Hamilton-Jacobi equation, 



gi"d l Jd v I + m 2 = 



(2.22) 



where m 2 represents the invariant mass; (b) reconstruct the whole action J, starting 
from the symmetries of the problem and the partial derivatives d^I, by means of 



where the integration is carried along an oriented, null, curve (like for example abc in 
figH]), to which we refer as the tunnelling path, with at least one point on the horizon; 
(c) split the integration along the null path in two pieces, one along the segment 
crossing the horizon (aS in figHJ), the remaining contribution living in the outer domain 
of the space-time (6c in fig[T|); (d) perform a near- horizon approximation in the first 
integral and regularise its divergence according to Feynman's ie-prescription. 
Once the procedure is done, the imaginary part of the classical action, relevant for the 
tunnelling rate (|2.9p . reads 

Im/=™, (2.24) 

K 

where n and w are the surface gravity and particle's energy, respectively. 

As we have already stressed, Hawking radiation is a semi-classical result of quantum 
field theory in curved space-time. In the spirit of tunnelling, we allow particles to travel 
along classically forbidden trajectories by starting just behind the horizon onward to 
infinity. Thus, particles must travel back in time, since the horizon is locally to the 
future of the static external region. For this reason, we must implement coordinates 
well defined through regions / and //, a requirement which automatically discards 
Schwarzschild-like or isotropic-like coordinates (it is easy to show that, in Schwarzschild 
space-time, isotropic coordinates, (t,p,9,(j)), defined through r(p) := (1 + M/2p) 2 p, 
cover regions I and IV of fig[T|). The classical one-particle action becomes complex 
as an indication of the classical impossibility of the motion and gives, in this way, an 
imaginary part to the emission rate of the hole. In this sense, we can say that, of 
the actual path travelled by the tunnelling particle (e.g. abc), only the infinitesimal 
region across the horizon plays a crucial role in the whole calculation. Even if null 
classically forbidden curves like at do not have anything special with respect to other 
forbidden paths, their choice is preferable to computational purposes. Also, it will 
become clearer later the reason why we specifically regularise divergent integrals by 
Feynman's prescription rather than others. In consideration of our choices, the mass 
term in Equation (|2.22l) is irrelevant to the physics of the horizon (Hawking radiation) 
but, as we shall show in Section[Sl under quite general conditions particles masses play 
an important role in relation to bulk particle creation phenomena. 

To show in practice how the Hamilton-Jacobi method works, we are going to derive 
the emission rate (12.91) and prove the identity (I2.24[) in the case of a scalar, uncharged, 
particle tunnelling from a Schwarzschild black hole, in complete analogy with what we 




(2.23) 
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did in the previous Section. This will give us the opportunity to unveil other debated 
points of the method. Later on we shall free the derivation from any dependence by 
special coordinate systems. 



2.3.1 Painleve-Gullstrand coordinates 

In this case the space-time geometry is described by eq. (l2.7[) . Because of eq.i 
expressing the relation between Painleve and Schwarzschild times, the particle energy 
associated to a Killing observer is expressed by u; = —d tp I (according to (|2.4I) . the 
particle energy is uj = —p ■ £ where p p = — d^I is the four-momentum). Since we are 
dealing with a spherically symmetric space-time, we may neglect without fault the 
angular coordinates. We consider a null direction (At p , Ar), for which 

- (l - ^) At 2 + 2^-At p ■ Ar + Ar 2 = . (2.25) 

In a neighbourhood of the event horizon (e.g. al> in figfl}, the relation between At p 
and Ar is uniquely fixed to be At p = —Ar/2 since, in this limit, the alternative 
solution Atp = — (1 — y/rg/r)^ 1 Ar diverges in contrast to the physical meaning of the 

coordinate t p . It simply means that the segment is part of null direction across the 
horizon. The reduced Hamilton- Jacobi equation 



(d tp i) 2 + 2^d tp id r i + (l - r f) {d r if 



(2.26) 



can be written as 



uj 2 -2cjJ^d r I+ (l-^) {d r I) 2 = 0. (2.27) 



The imaginary part of the action is 

Iml = Im I (d r I dr + d tp I dtp) (2.28) 

J a— s-h-+c 

= Im J dr (d r I + (2.29) 

= Im J|C drd r I . (2.30) 

Passing from (|2.28|) to (|2.29[) , we took into consideration that: (i) only the classically 
forbidden part at of the whole null path abc contributes to complexifying the action 
of the tunnelling particle; (ii) as far as we are concerned with the infinitesimal region 
of intersection between the horizon and the path ao, dt p and dr are related to each 
other in the way specified above; (iii) the particle's energy goes as —dt I. Being uj a 
real positive constant, the second term in the integrand of (12 .29[) does not enter the 
imaginary part of the action so it drops down in passing to (|2.30l) . From now on, 
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integration along the classically forbidden path {e.g. ab or, by extension, abc) will be 
denoted simply by^. 

It is remarkable that also the integration over the "time" coordinate enters, in gen- 
eral, the reconstruction process of the particle's action. In this specific case, we see 
that this "temporal" contribution — being manifestly real — does not affect the re- 
sult. The next subsection will show a coordinate system where instead the "temporal" 
contribution vanishes exactly. It should be clear that names such as "temporal" or 
"spatial" contributions arc coordinate dependent terms in no way fundamentally re- 
lated to the physics of the process, which must be covariant. The general covariance 
of the tunnelling method will be continuously emphasised as one of the main themes 
of this review. 

Inserting the non-manifestly real solution of eq. (|2.27l into (I2.30|) , we get 



and regularising the divergent integral according to Feynman's ie-prescription, the 
imaginary part of the action becomes 



Strictly speaking uj should be computed on the horizon, because this is how the method 
of residues works; however this is the same thing as the energy measured at infinity 
since on the outgoing trajectory u) is conserved. That said, after the known identi- 
fication k — 1/AM, the identity (|2.24p is fully recovered. Weather we had used the 
opposite prescription, r — > r + ie, a corresponding change of sign of the imaginary part 
would have resulted. Its meaning will be explained in Section [3] 

2.3.2 Eddington-Finkelstein coordinates 

Let us introduce another reference frame without singularities on r g originally con- 
structed by Eddington (1924) and Finkelstein (1958), [55]. This frame is fixed to 
radially moving photons. Since no observer can move together with photons, this new 
frame is not, strictly speaking, a reference frame. Nevertheless, this system of test 
photons proves to be very convenient, |52j . According to (|2.ip . the equation of motion 
of a radial incoming photon is dr/dt = — (1 — r g /r). From the viewpoint of a distant 
observer, the photon, starting from r\ at time t\, arrives in r (r g < r < ri) at time 




(2.31) 




(2.32) 




(2.33) 



This expression can be opportunely re- written as v = V\, where 



v :— t + r, 



* ! 



r* := r + r g In 



1 



(2.34) 
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r* is the so-called tortoise coordinate with v± a constant characterising the initial data 
of the photon at (ti,ri). Because of the logarithm in (|2.33[) . r* is defined for any 
r > 0. Chosen a set of photons at fixed t, we may label each photon of the set through 
a number v, which will identify uniquely that photon during its whole motion: v rises 
to the role of a new null coordinate, usually called advanced time. After differentiation 
of (|2.34[) and substitution in (I2.ip . the line element takes the so called Eddington- 
Finkclstein form, 

ds 2 = - (l - ^ dv 2 + 2 dv dr + r 2 (d6 2 + sin 2 9 d(/> 2 ) . (2.35) 

Because of (|2.34j) . the Killing vector field £ p = dt transforms to = d v and particle's 
energy as measured by such (natural) observer is simply u) — —d v I. Favoured by the 
spherical symmetry of the problem, we may consider only radial motions. Expanding 
the trajectory along a null direction in a neighbourhood of the horizon {e.g. ~a% in 
fig IB), gives 

= 2 Av ■ Ar . (2.36) 
However, along the null path aS, the v coordinate is constant, so that 

Av| N =0 (2.37) 
is the right solution. The Hamilton- Jacobi equation, 

2d v Id r I + (l-^(d r lj 2 =0 (2.38) 

can be re-written as 

' 1 - ^ ) d r I = 2u . (2.39) 



r 

According to the reconstruction assumption, 



ImJ = Im / (d r I dr + d v I dv) , (2.40) 

J a— >b— >c 

and thanks to ((2371) and (EH, 



1- 



Im/ = 2Imt* 7 ft + / (■•■)• (2-41 



6— ¥C 



Neglecting the real contribution coming form the classically allowed path be and fol- 
lowing the same procedure as in (|2.3ip . we end with the same result as in Equations 
(|2~M1) and (12321 . 

Among the coordinate systems of a certain importance covering the space-time re- 
gion across the gravitational radius there is the so-called Lemaitre frame. In this 
frame the Schwarzschild space-time appears as truly dynamical. We shall see how the 
Hamilton Jacobi method deals with it in Section 14.2.11 

We end this section pointing out the role played by the local observer. As discussed 
in |123j . the particle's energy u> = —dtl given in the two previous examples is obviously 
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the conserved energy as measured by an observer living at infinity. It follows that the 
Hawking temperature of a Schwarzschild black hole, T = l/8irM, is actually the one 
measured by a Killing observer at infinity. The particle's energy as detected by a Killing 
observer at position x Q b is given instead by cu b = u/ y— goo( x ot>)- Comparing the 
tunnelling probability derived from (|2.24p with the thermal distribution as measured 
by the observer /3 & w & = 2ttuj/k, we obtain 

Tob V-3oo(x b) = ^ — constant , (2.42) 
which expresses the expected result given by Tolman years ago [66] . 



2.4 More general static solutions 

Up to now we have only considered the Schwarzschild solution in order to keep the 
discussion as simple as possible and to illustrate the principles involved without un- 
necessary complications. However more general static solutions are of interest for a 
variety of reasons. In order to include a broader class, the metric can be written in a 
diagonal gauge as 

ds 2 = -V(r)dt 2 + —— + CirfidO 2 + sin 2 6 d<j) 2 ) , (2.43) 
W(r) 

which describes what Visser [57] termed dirty black holes. Black hole solutions are 
defined by functions V(r) and W(r) having simple and positive zeroes. This is only 
a necessary condition to have a black hole; we must also require that the domain of 
outer communication be "outside of the black hole" i.e., it should correspond to values 
of the radial coordinate larger than the horizon and extending up to spatial infinity. 
Interesting black holes described by such metrics can be obtained in the Einstein- 
Maxwell-dilaton coupled system. An example illustrating the feature is the following 
two-parameter family of solutions 

ds 2 = - (l - dt 2 + (l - — ) 1 (l - — ) 1 dr 2 + r 2 {d0 2 + sin 2 6 dcf) 2 ) . (2.44) 

The dilaton is exp(2tp) = (1— r_/r) -1 / 2 ; the hole has magnetic charge q m — 3r+r_/16, 
horizon radius r + = 2M (which defines M) and can be extended to a non singular, 
geodesically complete solution with horizons and asymptotically flat infinities. 

As long as V ^ W in cq. (|2.431 it is possible to set C(r) — r by an opportune radius 
redefinition. But for V — W, this is not always the case as can be seen from the 
following example 

ds 2 = - 1 ~™ ,T dt 2 + ^ , dr 2 + VAr 2 (d6 2 + sin 2 0# 2 ) . (2.45) 
1 — 2M/r 



Here C(r) := A 1//2 r with A = 1 + 2Mr 1 sinh 2 7 and 7 a real constant. The metric 
is asymptotically flat with an event horizon at r = 2M and an electric field due to a 
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charge Q = Msinh27/2. This solution can be obtained from dimensional reduction 
of a boosted Schwarzschild solution in Kaluza-Klein theory. 

The metric ea. (|2.43[) admits also a Painleve-Gullstrand representation associated 
with a redefinition of the time coordinate, 



t v = t± I dr< / 1 , . (2.46) 

' V V(r)W(r) y ' 



A simple computation leads to 



ds 2 = -V(r)dt 2 p ±2 J 1 ■ V(r) dr dt p + dr 2 + C(r) 2 (d6 2 + sin 2 9d<p 2 ). (2.47) 

Turning our attention back to ca. (|2.43|) with C{r) = r, by computing the Einstein's 
tensor one sees that the stress-energy tensor must have the following form 

T» = diag (-g(r), P(r), S(r), S(r)). (2.48) 

Thus, in general, the metric ea. (|2.43[) has continuously distributed sources and is not 
asymptotically flat. If the space-time we have in mind is such that 

V(r) = { — Tl) ' ' ' (r " Ts) 



W(r) = fr-O-fr-r,) g ^ g (2 49) 

where n < T2 < • ■ ■ < r a < r q , then 

y(r)-Wir)>0{l r >^ S s = l : (2.50) 

which means that outside the outermost event horizon r = r q , V(r) and W(r) are 
both positive functions. In this case from Einstein's equations it is possible to prove 
that [69] 

V(r) = W{r) g + P = through space-time (2-51) 

which is satisfied only in vacuum space-times with cosmological constant. This condi- 
tion is completely equivalent to the fact that 68 in a static space-time, 

g tt g rr = -l «=> T^n v = Q, (2.52) 

for all radial null vectors . Another important fact regarding such dirty black holes 
is that as a consequence of the dominant energy condition and the Einstein's equations, 
if V{r) has a simple zero at some r = rn then W(r) also has a simple zero at the same 
point (for a proof of this and other properties see [69] ). 

The metric ea. (|2.43|) is of interest also because it admits two inequivalent definitions 
of conserved energy: one is the Killing energy, u = —p^^ = —dtl; the other is 
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the Kodama energy which uses the vector field K — y/W/V dt- This vector field 
comes about because, in spherical symmetry (as spelled out extensively in Section [4]), 
it has the amazing property that both K M as well as J u = —K^T^ U are conserved: 
V^-ftT^ = 0, V M J' 1 =0. As a consequence it is possible to define two different notions 
of temperature, depending on which energy one is using. 

As we shall prove in detail in Section [31 even if the metric eq. (I2.43[) is singular on 
the horizon, the tunnelling method works. Using for example the Hamilton-Jacobi 
version, the massless Hamilton-Jacobi equation for radial motion reads 

-(cV) 2 + W(r){d r I) 2 = 0. (2.53) 



V(r 



The classical action is given by (plus/minus sign corresponding to outgoing/ingoing 
particles, respectively) 

r dr 

I± = -ujt ± / u — , (2.54) 



y/V(r)W(r) 

where w = —dtl represents the Killing energy. Assuming the near horizon expansion 
(as noted above, the occurrence of the same zero in V(r) and W(r) is a theorem, given 
Einstein's equations and some energy condition) 

V(r) = V'(r H )(r - r H ) + . . • , W(r) = W'(r„)(r - r H ) + . • • , (2.55) 

and Feynman's prescription, we readily obtain 

Im J dl+ - Im Jdl- = — (2.56) 

where 



^y/V'{r H )W'(r H ) (2.57) 



is the Killing surface gravity. One may use the Kodama energy, ojk = \/W/V w, which 
on the horizon takes the form ljjc = y/W (rjj)/V' (th)oj via de l'Hopital rule; in this 
case one obtains the so called Hayward's surface gravity 

*H=\w\r H ) (2.58) 

and the corresponding temperature Th = Knjlix. As we shall see in Section 01 this 
is the one that can be generalised naturally to dynamical situations, where there are 
no more Killing vectors in general. Notice however that uj/k — lok/i^h, so that the 
tunnelling probability is invariant. 



2.5 Tunnelling of fermions 

At this point of the discussion, one might ask what particles are to be found in the 
Hawking radiation spectrum of a black hole. Since a black hole has a well defined 
temperature, in principle it should radiate all the standard model particles like a 
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black body at that temperature (ignoring grey body factors) . The emission spectrum, 
therefore, is expected to contain particles of all spins, in particular fermions as well. 
The basic reference here is the seminal paper by Kerner and Mann 70} , to which 
others followed for both stationary and dynamical black holes [7TH75] . The fact that 
masslcss fermions are emitted with the same temperature as massless bosons is not a 
trivial result, given the fact that fermionic and bosonic vacua are in general distinct. 
It is important to mention that these expectations have been recently extended to 
include spin-1 bosons; and that the Hawking temperature does not receive higher 
order corrections in h beyond the semi-classical one [7TJ[771[75] as originally proposed 
in [75] (see also [50H55] for further extensions). 

Now, what do we take as the action of fermionic particles? Undoubtedly the most 
convenient choice would be to take the phase of the spinor wave function which, as is 
well known, satisfies the Hamilton-Jacobi equation as for spin-less particles. However, 
one may also consider the action from another point of view. Virtually, all known 
variational formulations of the motion of spinning particles have an action of the 
form [SUES] 

If = Iq + (spin corrections) 

where Iq is the kinetic term equal to the classical action of scalar particles as considered 
in previous sections, and the correction terms contain the coupling of the spin degrees 
of freedom with the spin connection of the manifold. These can either be derived from 
covariance considerations or more directly by the semi-classical treatment of the Dirac 
equation itself. Some of these terms give additive corrections to the particle four- 
momentum but in no case the remaining terms contain horizon singularities, as they 
are only responsible for spin precession effects. In the following course, therefore, we 
shall ignore them. Next, we neglect any change of angular momentum of the black hole 
due to the spin of the emitted particle. For zero-angular momentum black holes with 
mass much larger than the Planck mass (mpi ~ 10~ 5 gr) this is a good approximation. 
Statistically, as many particles with spin in one direction will be emitted as particles 
with spin in the opposite direction, producing no net change in the black hole angular 
momentum (although second-order statistical fluctuations will be present in general). 

As an example, we consider the tunnelling of fermionic particles across the event 
horizon of a spherically symmetric, uncharged, static black hole described by the metric 
(|2.43p with C(r) = r; and compute the emission rate in a singular Schwarzschild-likc 
coordinate frame. We refer to the original paper (70) for analogous discussions in 
generalised Painleve-Gullstrand and Kruskal-Szekeres coordinates. 

Temporarily re-introducing H, the Dirac equation in curved space-time reads, 



where 



(^D^ + j)v(t,r,9, ( f>)=0, 

£ Q/ 3 = jbcclp] — 



(2.59) 

(2.60) 
(2.61) 
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The 7 M — matrices satisfy the Clifford algebra, 

where I is the (4 x 4)-identity matrix. For this case, we pick the 7 matrices 



(2.62) 



1 
1 



7'' = #Wl .3 



a 3 
a 3 



7 



1 / 



a 1 



f 



r 



r sm ( 



a 2 
a 2 



where cr's are the Pauli matrices satisfying usual relations, 

o l <jj = 15ij +ie ljk vk, i,j,k = 1,2,3 
The matrix for 7 s is instead 



(2.63) 



5 ■ t r <t> 

7 = 7 7 7^ = li 



! W(r) 1 



1 

1 



V(r) r 2 sin ( 

For the spin-up (positive r-direction) Dirac field, we employ the ansatz 



( H(t,r,9,cf>) \ 


Y(t,r,9,4>) 




exp 



(2.64) 



We will only solve the spin-up case explicitly since the spin-down (negative r-direction) 
case is fully analogous. Employing the ansatz (|2.64l) into the Dirac equation (|2.59[) 
results in 



iH 



iY 



d t I t + Yy/W(r)d r It + Hr 



Y 

r 



d <t>h 



: d t I t - H v / W(r)d r I t +Ym = 



-— ( d e If 



d <t>h 







to leading order in h. As we assume that the action takes the form 

I t = -Ljt + T{r) + J{6,(l)) 



(2.65) 
(2.66) 
(2.67) 
(2.68) 

(2.69) 
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these yield the set of equations 
/ i uj H 



- Yy/W(r)T'{r) +mH = (2.70) 



r \ sm# 



H^/W[r)T'(r) ]+Ym = (2.72) 



•— (je + ^-zJ* ) = ( 2 - 73 ) 



Regardless of H and Y, Equations ((2T7T|) and (|2~73"|) result in J e + i(sin6») _1 = 0, 
implying that J(9, <ft) must be a complex function. The same solution for J is obtained 
for both the outgoing and incoming cases. Consequently the contribution from J 
cancels out upon dividing the outgoing probability by the incoming probability as in 
Equation (|3.6p ; and we can ignore J from this point. Equations (|2.70[) and (|2.72[) (for 
m = 0) have two possible solutions: 

H = -iY, T'{r) = T' ou , 



H = iY, T'(r)=Tl n (r) = - 



^V(r)W{r) 



where Font,™ corresponds to outward, inward solutions. The overall tunnelling prob- 
ability is r ~ exp —2 Im (T ou t — Tin), with 



Toutir) = -T in {r) = %dr . (2.74) 

J y / V(r)W(r) 

Let us remind readers that under the dominant energy condition and Einstein equa- 
tions, the functions V, W have the same zeroes. Therefore near r# we have, to first 
order, 

V(r)W(r) = V' {r„)W' {r H ){r - r H f + ■■■ 

and we see that there is a simple pole with a determined coefficient. Using the Feyn- 
man's prescription we obtain 

Im (T^ - T in ) = Im i dr ^ = — (2.75) 
J \/V(r)W(r) k 

where the surface gravity eq. (l2.57[) (see also eq. (|3.36l) is recovered, namely 

K=\^V'{r H )W'{r H ) . (2.76) 

ea. (l2.75l is a special case of a general identity we shall prove in Section [3] Alternatively, 
regularising as fully explained in Equation (|3.34[) we end with the same outgoing action 
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of (|3.35[) and the same surface gravity. Solving equations (12.70[) and (|2 . 72[) for H and 
Y in the case that m ^ leads to the result: 

fH\ 2 _ -iuj+ ^V(r)m 
\YJ ~ iuj + ^/V{r)m 

and approaching the horizon, we see that lim r ^ rjf (H/Y) 2 = — 1. Following a pro- 
cedure similar to what was done above, we obtain the same result for the Hawking 
temperature as in the massless case. The spin-down calculation is very similar to the 
spin-up case discussed above. Other than some changes of signs, the equations are of 
the same form as the spin-up case. For both the massive and massless spin-down cases 
the Hawking temperature k/2tt, with k provided by (|2.76|) . is obtained, implying that 
both spin-up and spin-down particles are emitted at the same rate. This is consis- 
tent with the initial assumption that there are as many spin-up as spin-down fermions 
emitted. 

With this calculation, we basically end the whole story about tunnelling in static, 
spherically symmetric black hole space-times. We have seen how the tunnelling picture 
arose and developed in a variety of different proposals; in Section [3] we shall confront 
with more technical aspects as the equivalence between so called null geodesic and 
Hamilton- Jacobi methods and their mathematical foundations. For now, however, it 
seems us to be more important to move on and see how the tunnelling method works 
in conditions where, for example, time starts to play an active role. 

2.6 Axis-symmetric stationary black holes 

The generalisation to a less symmetric scenario than the spherical one has not only ob- 
vious motivations, being the static spherically symmetric situation not so realistic, but 
also of fundamental relevance. Accretion processes occurring naturally in astronomical 
stages are able to spin up a black hole [571155] as well as to determine the growth in 
time of the black hole itself. In this subsection we are not concerned with the latter 
effect (being substantially a non-stationary process, it will be treated in Section 2]), 
but rather we would like to point out some efforts that have been made to extend the 
tunnelling approach to stationary black hole space-times. 

The most immediate generalisation of the Schwarzschild spherical symmetry is the 
axis-symmetric one, i.e. the Kerr solution for a rotating body with mass M and 
specific angular momentum a = JM . The subsequent inclusion of electric charge 
brings us to the Kerr-Newman solution, which still observes the axial symmetry and 
can be seen as a rotating Reissner-Nordstrdm black hole. We will apply the tunnelling 
method in order to retrieve the emission probability for scalar and fermionic particles 
in these cases, highlighting some features that Kerr and Kerr-Newman solutions show. 

2.6.1 Kerr black hole 

The Boyer-Lindquist form for the metric of a stationary rotating black hole — corre- 
sponding to the choice of a time-independent reference frame which is asymptotically 
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a Lorentz frame at infinity — is given by 

ds 2 = - (1 - 2Mr ) dt 2 _ 4Mras^e dt ^ + S ^2 + £ ^2 + 



sin 2 6># 2 (2.77) 
where we have defined, 

E := r 2 + a 2 cos 2 , A := r 2 + a 2 - 2Mr . (2.78) 



The roots of the equation A(r) = 0, r± = M ± V M 2 — a 2 , identify the event (outer, 
plus sign) horizon and the Cauchy (inner, minus sign) horizon. 

In order to study the tunnelling process of massless particles from the rotating black 
hole, however, we consider a metric where the world lines are those of photons with 
energy uj moving at infinity with constant and whose projection of the angular 
momentum on the rotation axis of the black hole is L z = aw sin 2 0. In these Kerr 
ingoing coordinates |89j . obtained from the transformation 

dv^dt+^—^-dr , d(j) = dj> + ^-dr (2.79) 
the line element (|2.77l) reads 

ds 2 = - (1 - 2ft) dv 2 + 2dv dr + S dff 2 - 4°M™n a 6 d(j) dy 



-2a sin 2 9 d<j>dr + ir2+a2)2 ~/ A sin2 6 sin 2 9 dcf> 2 . (2.80) 

We recall that in a rotating black hole the static limit surface r st , that is the external 
boundary of the ergosphere given by the equation g tt = 0, does not coincide with the 
event horizon. Some remarks regarding the validity of the semi-classical approach have 
been raised |90j because the geometrical optical approximation is reliable in proximity 
of r st , not the event horizon, where the tunnelling is supposed to occur. In some 
works [90191) the problem is circumvented thanks to a co-rotating change of coordinates 
(j) cp — £l H t, where Vl H — a/ (r+ + a 2 ) is the angular velocity of the horizon. In the 
following we will show that this is not necessary and the right result is recovered 
automatically. 

As fully explained in the previous sections, in the semi-classical approximation, 
the relevant quantity that plays a role in the calculation of the transition probability 
through the horizon is the classical action / satisfying the relativistic Hamilton- Jacobi 
equation. Being the metric eq. (|2.80l) independent of v and (f>, a suitable ansatz for the 
form of the action is given by I = —ojv + J <p + J 7 (9,r). We are going to see that 
the imaginary contribution is due to the r-dependent part of the action, which will 
produce a pole in correspondence of the horizon. Expanding the Hamilton-Jacobi 
equation with the ansatz for / above, we obtain 

a 2 {in esc 6 — wsinfl) 2 + 2 [a 2 m — [a 2 + r 2 ) w] T r + 

+ [a 2 + r(-2m + r)] T 2 +F$ = (2.81) 



2G 



where subscripts r,9 are for d r> e- Solving for F r {9, r), we immediately note that the 
function to be integrated is also ^-dependent: 



X(r) ^X(r) 2 -A(r)[a 2 (Mc S c9-u JS m8) 2 +^ e (9,r) 2 } 
^ r) = -A^) ± AW ( } 

where X(r) = a 2 M — (a 2 + r 2 ) u). A way to deal with this dependence is to simply fix a 
constant value 9 — 9q and show that eventually the result has no effective dependence 
on the choice of 9o 05]. Actually there is no need to fix 9 because in order to apply 
the method we consider the regime in which the tunnelling occurs near the horizon, 
A(r+) = 0: in this situation the term that brings the ^-dependence drops out and the 
function depends only on r. The root of this result is the complete separability of the 
Hamilton- Jacobi equation in Kerr space-time. Recasting all the expressions in terms 
of r + and r_, so that X(r) = (r+ +r_)a 2 /2 — cj(r 2 +a 2 ) and A(r) = (r — r + )(r — r_), 
we get 



immr) = - Im I m±Vm -(r-r + )(r-r-)(...) dr 

J (r-r_J (r-r+) 

where (...) contains the whole ^-dependence. Regularising the integral by Feynman's 
prescription, 



ImJ = -2tt 



X(r + ) 

r + — r_ 



{r\+a 2 ) a 2 (r + 



r_) (r+-r_) 



(2.84) 



Reintroducing J = a(r + + r_)/2 and the angular velocity of the horizon tt H , we can 
rearrange the terms and obtain the tunnelling probability rate: 

exp [-2 Im/] = cxp [-/3 Kerr {cu - n H J)] (2.85) 

where the inverse temperature is /3k err = 47r(r 2 _ + a 2 )/(r + — r_ ). In the limit a — > 
we recover the result for the Schwarzschild black hole. 



2.6.2 Kerr-Newman black hole 

The tunnelling of scalar and fermionic, electrically charged, particles in the family 
of Kerr-Newman space-times does not present new features with respect to previous 
discussions, the black hole charge Q and specific angular momentum a parameters 
being treated on the same footing. In the following, we are going to consider first the 
tunnelling process for a scalar particle of charge q and energy oj, then for a fermion. 
In both the examples, we implement singular Boyer-Lindquist coordinates (t,r,9,<f>) 

ds 2 = _[ A-.W» j dt 2_ 2(rW-A) «sin^ ^ ^ + g ^2 + g rf g2 

sin 2 9 dcj) 2 (2.86) 
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where both incoming and outgoing trajectories are expected to contribute to the emis- 
sion rate, according to our master Equation (|3.4p . The symbols in (|2.86[) are only 
slightly changed with respect to above: A := r 2 + a 2 + Q 2 — 2Mr. The static limit 
surface is now located at r st = M + y M 2 — Q 2 ~~ a 2 cos 2 9, while the event horizon, 
still the greater root of A(r) = 0, is r + = M + \J M 2 - Q 2 -a 2 . 

Tunnelling of scalar particles As in the previous case, the action can be separated 
as / = — ujt + J(f> + J- {9, r) and further separation of the last terms occurs near the 
horizon. We explicit the metric functions inside Hamilton-Jacobi equation 

gT (dpi - qA^j (d v I - qA v ) = (2.87) 

where the vector potential one- form is given by A = — Q r Y<~ 1 ^dt — as'm 2 9d(j)j, to 
obtain 

„ _ [(r 2 +a 2 ) 2 dtl+a d^I-qQr] 2 (d^I+a sin 2 9 d t lf 

U ~ 2EA 1 2S sin a 8 + 

+|(Vf + ¥ (2.88) 
Making use of the ansatz for / and after some manipulations, we solve for J>(#,r), 

J \ (a 2 +r 2 )uj- a 2 M - qQrf - A(r)(. . . ) 

d r T in ,out(e, r) = ±^1 1 J (2.89) 

A(r) 

where the plus (minus) sign corresponds to outgoing (incoming) trajectories. Again, 
the (...) under the square root contains all the ^-dependence and near the horizon it 
will be negligible. Equation eq. (|2.89p has to be integrated along the path crossing the 
horizon, so let us write A(r) = (r — r + ){r — r_) in order to show explicitly the pole 
in the integrand and choose the plus sign to select a path that comes out of the black 
hole. The imaginary part due to the pole, rearranging a little bit the terms, is given 
by 

v f^ + a ^f— -^— 2 Ma - (2.90) 



2 \r + - M 

where we easily identify the angular velocity Sl H , the angular momentum parameter 
J and the term arising from the presence of a charge parameter: this last term is 
expressible in terms of the electric potential $ = Qr^-/(r 2 _ + a 2 ) of the black hole. 
In order to retrieve the final expression, we take into account also the contribution 
coming from the ingoing trajectory, which is given by the very same procedure as 
before applied to the function Equation eq. (|2.89[) with the minus sign. The imaginary 
contribution is found to be equal and opposite in sign so the tunnelling probability 
will be given by 

^ = exp [-2 Im {T^ - Fin)] = e~ 4 . (2.91) 

t ab 
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A more rigorous derivation of this formula will be given in the next section. The final 
expression for the tunnelling probability can be written as 



from which the temperature can be read /3g N = Tkn = (?"+ — r_)/47r(rl + a 2 ). It is 
easy to see that two subsequent limits Q — > and a — > lead us first to the Kerr and 
then to the Schwarzschild solutions. 

Tunnelling of fermions As a further extension, it is interesting to consider the tun- 
nelling of gravitinos (spin-3/2 fermions predicted in supergravity theories 92 ) in Kerr- 
Newman black hole space-times [7H[7S]. The result is analogous to the scalar case, 
suggesting again evidence for the universality of black hole radiation, and follows the 
same line of Section 12.51 where now only little more shrewdness in the choice of the 
representation for the Dirac matrices is needed. 

The equation of motion of spin-3/2 particles, namely the Rarita-Schwinger equation, 
can be written as 



where = W^t, is a vector valued spinor of charge q and mass m; matrices 7^ satisfy 
the Clifford algebra; and the covariant derivative is defined as in (|2.60[) and (|2.61[) . 
The Rarita-Schwinger equation presents as a Dirac equation applied to every vector 
index of the field W and a set of constraints preventing ghosts propagation. In the 
present context, the anstatz for the wave function is — (a^, 6 M , c M , d M ) T e 11 ^, where 
T is for transposition. It turns out that the Dirac equation can be solved for the 
action / independently from the components a M , d M ; conversely, 7^Vp = yields a 
set of constraints for the components a^, ...,d M , independently by the action /. Thus 
we may conclude that to calculate the Hawking temperature of gravitinos emitted 
from the Kerr-Newman black hole, we only need to solve for the Dirac equation in 
precise the same way as in Section 12.51 and that, being the action unaffected by the 
constraints, fermions with every spin will be emitted at the same temperature. We 
refer the interested reader to the original paper [74] for the precise form of the Dirac 
matrices suitable to the geometry at hand. 

2.7 Further extensions 

In the literature, the tunnelling method has been extensively applied to generalizations 
of the stationary black holes. Among the plethora of models that have been considered, 
we choose to focus our attention on few interesting cases, comprehensive of non axis- 
symmetric solutions as well as higher and lower dimensional cases. Our main interest 
in considering such extensions is to point out features that would be difficult to get by 
quantum field theory techniques. 




(2.92) 



%-f (Dv + iqA^V^ + j-V^O 
7^ = 
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Rotating and accelerating black holes The metric 76, 93, 94 can be written in the 
form 

dr 2 

ds 2 = -f(r, 9) dt 2 + ——- + E(r, 9) d9 2 + K(r, 9)dcf> 2 - 2H(r, 9) dtd<p (2.93) 
9(r, 9) 

where 

f(r,9) := Q - P /J m2e , g(r,6):=^, B(r,0) := ^ 

K(r 9) ;= ^^WfV sin^]^ H{r 9) :=2asin 3gP(r-+g- Q 

and further, 

ft :=l-arcos9 , Q := [r 2 - 2Mr + (a 2 + e 2 + g 2 )] (1 - a 2 r 2 ) 
£ :=r 2 +a 2 cos 2 6» , P := 1 - 2aM cos 9 + a 2 (a 2 + e 2 + g 2 ) cos 2 9 

This metric is described by the parameters M, a, a, e and g - respectively the mass, 
the angular momentum, the acceleration, the electric and magnetic charges of the 
black hole. The vector potential for these solutions is given by A = —^ 1 [er(dt — 
a sin 2 9d(f) + g cos 9(adt — (r 2 + a 2 )d<p)] . Inspecting £1 = Q = 0, we see that, apart from 
the inner and outer horizons located at r± = M ± (M 2 — a 2 — e 2 — g 2 ) 1 ^ 2 , this black 
hole has two acceleration horizons r a ± = ia" 1 and r a = (acos#) _1 . The calculation 
of the tunnelling probability from the outer horizon gives a temperature 



T = 1 (l- a r + f{l + ar + ^M 2 -a 2 -e 2 -g 2 ) 
+ 2tt 2M(M + y/M 2 - a 2 - e 2 - g 2 ) - e 2 - g 2 

while for the acceleration horizon, given by the same expression with the substitution 
r + — > r a , one gets a vanishing temperature, a confirmation that "the acceleration 
horizon is not a black hole horizon essentially" |76j . 

Kerr-Godel space-time A second example comes from the Kerr solution embedded 
in five-dimensional Godel universe, that is the Kerr-Godel solution [48,95 . Among 
several useful ways to write down the metric for this solution, we can start from the 
lapse-shift form 

ds 2 = -N 2 dt 2 + g(r) ( a - ^jdt) 2 + *1 + \ r 2 (, 2 + ,|) (2.95) 

and the vector potential A — >/3j r 2 q/2 where 

N 2 := ^ , 9 (r) := -jV + \(1 - 8Mf)r 2 + f£ 
v{r) ;= 1 _ iM^efM^Mi + 2_Mf a{r):=jr 2 + m j 

the Cj are the right-invariant one-forms on SU(2) 

<Ti = sin 4> d9 — cos <j> sin 9 dip , ft = cos <j> d6 + sin sin d<0 , ^3 = d0 + cos 
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and the parameters j and I are respectively the Godel parameter (determining the 
rotation of the space-time) and the black hole rotation parameteiQ: when I = the 
metric reduces to the Schwarzschild black hole embedded in Godel universe; when 
m = I = it reduces to the five-dimensional Godel universe; for j = we recover the 
five-dimensional Kerr black hole. Peculiarity of this solution is the existence of closed 
time-like curves inherited from the Godel universe. In fact, in the suitable range of 
values for which g(r) < 0, we have < and becomes time-like: in this case the 
curve with constant t, r and 9, being (f> a periodic coordinate, becomes a closed time-like 
curve. The main global properties of this space-time are the (possible) existence of: (i) 
a closed time-like curve horizon defined by the condition g(r) = (being N 2 divergent 
there, it means that nothing can cross the closed time-like curve horizon starting from 
the region without closed time- like curves); (ii) black hole horizons located at the 
positive roots of V(r) = 0; (hi) the ergosphere given by f(r) = 0, that is r erg — \JlM. 

We consider the parametrisation J = j^/8M with J 2 < 1: in this case the closed 
time-like horizon will be the outermost surface beyond which closed time-like curves 
occur. In this region of parameter space, the event horizon of the Kerr black hole can 
be thought simply as embedded in a Godel space-time 

We want to show the result for the tunnelling in this solution using the null- geodesies 
method, so it is more convenient to cast the metric in a Painleve form [35] : first of all 
we move to co-rotating coordinates, choosing geodesies with zero angular momentum 
(j) = Q t and, for convenience, also d9 = dip — 0; then we pass in Painleve-like form 
with the transformation 



t t p = t-2ldr ^ ( ?- Vir)) (2-96) 
' rV(r) 



The final form of the reduced line element is 



Ag(r) p g(r) 



The radial null geodesic equation to be used in order to calculate the imaginary part 
of the action is given by 



(±1 - y/l-V(r)} (2.98) 



where the plus (minus) corresponds to outgoing (incoming) trajectories. The expres- 
sion develops a divergence at the horizon, leading to a pole in the integral for the 
reconstruction of the action /, with residue 4y^(r) JthV 1 {th) on the horizon. The 
tunnelling probability turns out to be 



r em ocex P u) (2.99) 

\ r H V'(r H ) I 



4 The five-dimensional Kerr solution has usually two rotation parameters h,l2 I96| : in the case 

considered here we have ii = '2 = 
5 It is not clear what the parameter region J 2 > 1 represents physically. This case, when the closed 

time-like horizon is the innermost surface, is analysed in |48| . 
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from which the temperature can be read off 



M 



-KG 



•%(l-8fM-4Ml) 



21' 



Trr%yJ-Aj 2 r% + (1 - &j 2 M)r% + 2M1 2 



(2.100) 



For I = it reduces to the Schwarzschild-Godel case [ST]- The application of the 
method to the closed time-like horizon gives a vanishing temperature: in fact, from 
Equation eq. (l2.98[) we see that r — > oo when g(r) — > 0, so the integral has no poles 
there. This means that there is no tunnelling through the closed time-like curve 
horizon, a very interesting result that would be difficult to obtain by quantum field 
theoretic methods. 



Taub— NUT— AdS black hole As a third example, we recover the temperature asso- 
ciated to the horizon of a Taub-NUT-AdS black hole [46] . whose general form for the 
metric can be written as [581 



ds 2 



-F(r) 



dt + An 2 f 2 [ - ) d(j) 



dr 2 

+ -FT-,+^ + n 2 ){de* + fl{eW) (2.101) 

F{r) 



where 



F{r) := k{r 2 - n 2 )(r 2 + n 2 )" 1 + (r 2 + n 2 ) - 2Mr + l 2 {r 4 + 6n 2 r 2 - 3n 4 ) 



A = — 3 [ 2 is the cosmological constant, n is the so-called NUT charge and the 
parameter k determines the explicit form of the function fk(&), 



fk(0) 




if k = l, 
if k = 0, 
if k = -1 



The Taub-NUT solution represents a generalisation of the Schwarzschild space-time 
(for details see [55j 98 ) and, as in the previous example, also these solutions possess 
closed time-like curves for some ranges of the parameters involved. Here, we shall 
explicitly refer to the case k = — 1 and 4n 2 < [ 2 , that is a hyperbolic subclass that does 
not contain closed time-like curves. In this case, in order to apply the Hamilton- Jacobi 
tunnelling method we expand the metric functions near the horizon, F(r) ~ F'(r + ) (r— 
r_|_). The metric ea. (12.101j) can now be written in the near horizon approximation as 



ds 2 



-F'(r + )(r - r+)dt 2 



dr 2 



F'(r+)(r-r+) 



G{r+, 



(2.102) 



The Hamilton-Jacobi method can be applied straightforwardly, leading to a ^-independent 
expression for the tunnelling probability T em ~ exp (— Attuj / F' (r + )) and hence a tem- 
perature Ttn = -F 1 '( r +)/47r. This result can be compared to the one calculated by 
Wick rotation [98j|99] (to compare the two results, it is sufficient to Wick-rotate the 
NUT charge n 2 = —N 2 ). 
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BTZ black hole As a further application of the tunnelling method in stationary 
black hole space-times, we consider the Bahados-Teitelboim-Zanelli (BTZ) solution 
[100H102] . which represents a (2 + l)-dimensional rotating black hole with negative 
cosmological constant. The line element is written as 



ds 2 = -N 2 dt 2 



dr 2 
N 2 



r 2 (N^dt + d^f 



(2.103) 



where 



N 



I 2 



M ■ 



with the cosmological constant A 



J 2 J 
4r 2 ' 2r 2 ' 
-I -2 . The line element shows two horizons in 



correspondence of the zeros of N , which can be written as 



r% = -Ml 2 
± 2 



1 ± 1/1- 



J 2 
M 2 \ 2 



(2.104) 



)(r — r_)/i r . We maintain the condition Mi > J, which 



so that N 2 = (r 2 - 

ensures non-extremality. The two Killing vectors of the space-time, d t and allow 
the separation of the action as I = —uit + J<\> + T(f). The Hamilton-Jacobi equation 
in this case reduces to 



J 



(w + JN^f 
N 2 



+ N 2 F r (r) 2 = 



from which the function J-'(r) is easily recovered 



F' ± {r) = ± 



^/(ruj + JrN't' - JN)(ruj + JrA* + JN) 
rN 2 



(2.105) 



(2.106) 



The calculation of the imaginary part of the action goes on taking into account that, 
being the metric singular on the horizon r + , we have to consider both outgoing and 
ingoing trajectories. From the considerations already made, we know that it is suffi- 
cient to calculate the contribution coming from the outgoing path and then consider 
—2 Im(J 7 4- — J--) = —4 ImJ r + . Using Equation eq. (|2.104[> to show explicitly the pole 
of the integrand, the result is found to be 



exp [—4 Im {J-+)] ~ exp 



2nl 2 r 4 



2r 2 



(2.107) 



We recognise the angular velocity of the horizon il H = J/ir 2 ^ and infer that the BTZ 
temperature is Tbtz = (r+ — r 2 _) /2ir[ 2 r + . 



Topological black holes In the simplest form these are asymptotically Anti-de Sitter 
metrics of the form ea. (|2.43p with C(r) = r and 



V(r) = W(r) = k- — 
r 



— k 



0,-1 



(2.108) 
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except that the transverse manifold is either a torus (for k — 0) or a compact Riemann 
surface with higher genus (for fc = — 1) |103til06j : the case k = 1 is the Schwarzschild- 
AdS solution with spherical symmetry and spherical horizon. In all cases the horizon 
inherits the topology of the transverse manifold but, independently of this, the horizon 
pole responsible for the imaginary part of the action is given by the outermost zero of 
the function V(r), or equivalently as a root of the cubic algebraic equation rV(r) = 0. 
The derivation of the tunnelling probability given for spherical symmetry does not 
actually depend so much on the symmetry of the transverse manifold, what really 
matters being the warped form ds 2 = dj 2 + r 2 da 2 . In this way one recovers the 
Hawking temperature, Tjj = V'(rji)/iw, from the tunnelling picture for AdS black 
holes of any topology. It is worth to mention at this point that the local character of 
the tunnelling method makes the temperature and the rate independent of the effects 
of "large isometries" , when they are used to construct new solutions by the quotient 
method (identification of points). Thus, for example, the temperature of the RP 3 geon 
is the same as the Schwarzschild black hole. The universality of the Hawking effect is 
all the more evident in the tunnelling picture. 

Extremal black holes It may happen that for certain values of the macroscopic black 
hole parameters the corresponding surface gravity vanishes; such solutions are called 
extremal black holes since, as a rule, these parameters are on the line separating black 
holes from naked singularities. Such is the case of the Reissner-Nordsrtdm solution, 
which describes a black hole if and only if M 2 > Q 2 , or of the Kerr solution for values 
M 2 > J 2 /M 2 , or of the more general Kerr-Newmann family. These are the extremal 
black holes of general relativity without cosmological term A. There are many other 
examples either in the presence of A, or in theories involving coupled gravitational, 
electromagnetic and scalar fields [107,108 , but we will not cover all such solutions since 
for certain values of the dilaton coupling (the term exp(— 2a(p)F 2 in the Lagrangian), 
the thermodynamic description breaks down [10911110] . As a warning, the extremal 
black holes considered here must be distinguished from the extremal solutions which 
are obtained by means of limits taken on the parameters of the metric (such as the 
limit Q — ¥ M in the Reissner-Nordstrom solution) further accompanied by coordinate 
changes becoming singular in the limit (see, e.g. |lllj and references therein). Such is 
the case of the charged Nariai solution [112H114] describing de Sitter black holes with 
the maximum mass for a given charge, which are finite temperature objects. 

As a rule the horizon of extremal black holes is at infinite spatial distance, cor- 
responding to a double pole in the metric component g rr , but there are exceptions. 
One example is the family of extremal black hole solutions in dilaton gravity with 
electromagnetic coupling exp(— 2atp)F 2 , a > 1, which reads [108] 

ds 2 = _ (l _ + (l _ dr 2 + r 2 (i _ ^ 3£ d tf . 

This metric has zero surface gravity, an event horizon at finite spatial distance, but 
also a curvature singularity at r = To. 

The tunnelling method has been originally applied to extremal black holes in [4l)ll9"l"] , 
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where several problems were pointed out. For example, consider an extremal Reissner- 
Nordstrom solution 



which can be either magnetically or electrically charged, with charge Q = ±M in 
both cases. Kerner and Mann [46] showed that there is a breakdown of the WKB 
method in the null geodesic approximation, invalidating the method. The fact is that 
the quantity r out , so crucial for the null geodesic method, actually does not exist since 
upon emitting a neutral particle the extremal hole becomes a naked singularity. They 
also showed that both methods yield a divergent real part to the action, which is 
suggestive of a full suppression of particle emission since then there is no trace of an 
imaginary part. Further, taking into account self-gravitating effects, they were able to 
show that the temperature approaches zero the closer the original black hole (that is 
prior to the emission of a particle) is to extremality. Angheben et al. [SH] found the 
same result using the proper spatial distance to treat covariantly the horizon pole. On 
the other hand, using the radial coordinate one can easily find a prescription giving a 
finite result. From the metric one is led to study the integral 



Note the double pole, still present when using regular advanced time coordinates. 
Making use of the formula 



one gets a non vanishing imaginary part. Using instead the principal value, the integral 
will be real. How should one interpreted these results? It is known that an extremal 
Reissner-Nordstrom black hole will not radiate away its mass by emission of neutral 
particle^], as can be seen either by taking the extremal limit of the emission rate 
of non extremal black holes, or by computing the spectrum starting directly with the 
extremal solution. Therefore the temperature vanishes and one has to use the principal 
value prescription in order to be consistent with the (free) quantum field theory result, 
or alternatively a regularisation of the double pole giving a null result. The formal 
reason is easy to understand: for a simple pole the horizon is reached logarithmically 
fast as log(r — Th), so to pass the horizon inward necessarily generates an imaginary 
part ±?7r, while for a double pole the behaviour of the integral is (r — m)" 1 which can 
assume negative real values. The physical reason is that revealed by Kerner & Mann: 
there is no horizon's interior in a naked singularity. Of course one has to use the same 
prescription even for the emission of charged particles, so these too will not be emitted. 



6 It has been argued in 1311 that the emission vanishes even for charged particles if self-gravitation is 
taken into account. 




(2.109) 
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It is of interest that the hole still emits charged particle in the super-radiance channel, 
with total luminosity [11511116] 



1 



L = --f- V(2Z + 1) / Ti(uj)ujduj (2.110) 

271 1>Q Jm 

where q and m are the particle's charge and mass, respectively; = Q/M is the 
electrostatic potential at the horizon; and F;(w) the absorbitivity of the hole. Note 
that (i) is the sign of the hole charge since M — \Q\ in the extremal case; and (ii) 
the range of integration is within the super-radiance region, u> — q&+ < 0. Therefore 
the hole will preferentially emit charges of the same sign of the hole charge, thereby 
driving it away from extremality. We can interpret this result as cosmic censorship at 
work. 

Under further conditions toM 1(1+1), mM ^ 1 (in Planck's units) the absorbitiv- 
ity can be computed in the WKB approximation, which gives 



r,H = -exp^ — j, qQ>0. (2.111) 

Similar results and identical conclusions can be obtained for extremal rotating black 
holes. Here the sum in eg. (|2. 1101 is extended to cover the azimuthal quantum number 
m and the upper limit of integration is to be replaced with mil H , so that the hole will 
not emit radiation in the s-wave channel. 

We conclude that the tunnelling methods, at least in the present form, cannot de- 
scribe the super-radiance regime of black hole evaporation but are otherwise consistent 
with their stability as predicted on other grounds. Perhaps this means that with super- 
radiance we are dealing with a bulk phenomenon that should be treated by analogy 
with the Schwinger's pair creation effect by an electric field. 

Higher order corrections As our last point we report on a case which has been put 
forward in [TSUHO] and further extended in [5TW83"] . regarding corrections of higher 
order in h to the Hawking temperature. The argument starts from the semi-classical 
expansion 

I = I + J]ft n /„ (2.112) 

and the assumption that the field obeying the Klein-Gordon equation in the black hole 
background is as usual tp = exp (il/h). It follows from this that all /„ obey the same 
partial differential equation, giving rationale to the further assumption that /„ oc Iq 
for every n. Therefore, writing /„ = c n Io, the action takes the form 



l h n h (2.113) 



where Iq, to fix ideas, obeys the Hamilton- Jacobi ea. (|2.53l and {c„}„ remain unpre- 
dictable quantities. It is clear that the tunnelling probabilty V = exp (—2 Im7) will 
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predict a modified temperature 



T=! r — r (2.114) 



where 



T o = ^VV'(r H )W'(r H ) = — (2.115) 

is the uncorrected black hole temperature associated to the Killing surface gravity. 
Obviously this is not changed by simply passing to Kodama-Hayward surface gravity 
and energy. The procedure looked very suspicious to several authors [7Tl[77lfTT7) . 
who did not find any higher order correction neither from scalar emission nor for 
bosons or fermions. The contradiction was settled by Yale [75], who noted that the 
result eq. (|2.114[) assumes an incorrect definition of energy, namely loq = —dtlo, while 
Hamilton- Jacobi theory would require lu = —dfl instead. Then, from ea. (|2.113[) 



ui 



Cnh n wo (2.116) 



so that 



and all supposed corrections disappear. Of course, one should also allow complex 
solutions of the array of semi-classical equations, so even the assumption /„ oc Iq for 
all n is really unjustified. 



3 Analytic Continuation Arguments 

In the previous sections, we have described the Hamilton-Jacobi strategy in a list of 
four steps, from (a) to (d). It is clear, however, that, at least at first sight, not all of 
them stay on equal footing: besides some irremissible (e.g. postulation of Hamilton- 
Jacobi equation) or very natural requirements (e.g. trajectory splitting as abc into 
~ah + be), we find other less tolerable points. Why, in fact, should it be that particles 
traveling along classically forbidden trajectories from inside the black hole to outside 
must follow null paths? And which fundamental principle suggests us to regularise 
divergent integrals according to one prescription rather than others? By the end of 
the day, we shall show that — in contrast to the common sense even the Hamilton- 
Jacobi equation is an accessory requirement. A certain experience in the field tells 
us how these points can result in some sense cryptic to the same experts. With the 
purpose of clarifying some of the points mentioned above, we are going to outline the 
foundations of the Hamilton-Jacobi method in order to point out what is fundamental 
and what is only an additional assumption; and comparing, if possible, alternative 
solutions attempted in literature with present proposals. 
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3.1 Foundation of Hamilton-Jacobi method 



Let us consider the motion of a scalar particle from region II to region / in the eternal 
version of Schwarzschild black hole (at in figH]). This motion is classically forbidden 
since the particle should travel back in time to follow it. Notice that, in general, nothing 
we can say about the causal nature of the forbidden path. However, if the coordinates 
of the starting point a are displaced to complex values, then such an allowed path exists. 
Thus, the reader must think to fig|2]as a section of the complexified Schwarzschild plane 




Figure 2: A section of the complexified U, V plane of the eternal Schwarzschild black 
hole. 9, (f> coordinates are suppressed, abc is a path from mside to outside 
the hole, with a% a classically forbidden trajectory and a'b representing its 
reflected trajectory. 

with coordinates (9, cj>) constant and suppressed. By analytically continuing the point 
a in the complexified Schwarzschild space, the amplitude to propagate to c from a real 
point a in region II can be related to the amplitude to propagate to c from a reflected 
point a' in region III [54] . This latter process is just the time-reversed of absorption 
of a particle by the black hole. In this way, the emission probability for a black hole 
is related to the probability for it to absorb. All that we need in order to prove the 
case, is to take as a fundamental assumption the analyticity of the particle action in 
the complex (U, V) plane as a function of the space-time coordinates of point a |118j . 

In region 77 (III), the Kruskal-Szekeres coordinates U, V are both positive (nega- 
tive) definite. So, let us rotate U and V from positive to negative values by posing 

U = e M U, V = e~* 9 V, i?g[0,tt], (3.1) 

a choice which clearly preserves the invariant product UV = (1 — r/r g )e r ' r « . The 
integral of dl over the tunnelling path a — > b — > c will now be replaced by the integral 
over the path (say 7) consisting of the semi-circle a — > a with i9 ranging from to tt 
(over which \U\, \V\ are constants) plus the integral over the path a — > b — > c which 
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crosses the past horizon (over which is constant). In formulas 
dl = I (dfjl dU + dyl dV) 

(d v IdU + d v IdV)-i [ dd{V d v - U <%)/ , (3.2) 
Jo 

where, this time, the upward arrow denotes integration along the classically allowed 
trajectory a'bc. Since in region I, U = —e~ KU and V = e +KV , where u,v are retarded 
and advanced time coordinates and k is the horizon surface gravity, k = \/AM, we 
have that d t = n{Vdy ~ Udjj)- Hence we obtain 

dl= fdl -^9,1 = ^1+^, (3.3) 

with uj = —p ■ £ = —dtl the conserved Killing energy of the emitted particle. Taking 
into consideration only the imaginary part of the action, which is the relevant part for 
tunnelling purposes, we get 

-2 Im (J^dl- j»dl) • (3-4) 

Since the geometry is static, it is left invariant by time inversion T: 

Im Pdl = f Im f dl = Im / dl (3.5) 

J J Jb^a 

where b — > a is the time inverse of the path a' — > b (figH]). With a justified abuse of 
terminology, £ will denote also the integration along the time-reversed path. From a 
physical point of view, this is the path followed by an incoming particle absorbed by 
the black hole, so that exponentiating (|3.4I) . we can write 

r em = T abs e " , (3.6) 

in agreement with the result of Hartle and Hawking |54) derived by path-integral 
methods. 

The rotation of the (U, V) coordinates in the complex plane, Equation (13.11) . has been 
chosen in the form given because it corresponds to a Wick-like rotation of Schwarzschild 
time in the lower half complex plane, namely t t — i'd/n, which is the analyticity 
region of positive energy solutions of field equations. In fact: 



v v 



id 

-2k 1 1 

K 



(3.7) 



One might legitimately ask for the meaning of an anti-Wick rotation, t — > t + id/ n. In 
this case, one rotates U (V) clockwise (counter-clockwise) with consequent change of 
signs in Equation (|3.1[) . In turn, Equation (13.41) changes in favour of 

- 21m (ridl- fdl\=^ . (3.8) 
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The plus sign at the right hand side of this equation, tells us that the object we are 
talking about is not a black hole, but rather its time reverse (white hole), for which 
absorption is exponentially suppressed with respect to emission, T a i, s = T em e~~^ . 



The Reissner— Nordstrom solution It is easy to generalise the present discussion to 
include charged, spherically symmetric black holes. In this case, Einstein-Maxwell 
equations can be solved in terms of the Reissner-Nordstrom metric, 

ds 2 = -(\-— + %\dt 2 + -, — r + r 2 (d9 2 + sin 2 9 dtf) (3.9) 

V TV 2 ) U_w. + q\ v Y) 



and the Maxwell one-form, 



A=^-dt (3.10) 



where M and Q represent the black hole mass and electric charge, respectively. With- 
out entering the details of the solution, for which we refer the reader to standard texts 
[e.g. [55,119 ), it is enough to say that Reissner-Nordstrom black hole possesses an 
outer (inner) event horizon r + (r_ ) in correspondence of which we can still define a 
tortoise coordinate r*; null retarded and advanced coordinates u, v; Kruskal like coor- 
dinates U, V; and an outer (inner) surface gravity k + (k_ ). This time, the action to be 
integrated on the tunnelling path is dl = dlo + qA, where 1$ is the free action and q is 
the electric charge of the tunnelling particle. The form A is ill-defined at the horizon, 
for this reason one usually makes a gauge transformation to a form A A = A + dA which 
is regular there. In our case the analytic continuation takes A away from the horizon 
so this will actually be unnecessary. Using as above complex (U, V) coordinates, we 
easily obtain 

Im j£ dl = Imj* dl + ^-(u - g$) (3.11) 

where cj = —dtlo is the mechanical energy and $ = Q/r + . The quantity u— g$ is gauge 
invariant and conserved along the outgoing path from the past horizon to infinity. This 
is because for any Killing vector field £ such that the Lie derivative (£g A)^ = ^"d v A IM + 
d^ v A v — 0, the action I is invariant, to first order in a, under transformation x — > 
x + a^. The Nother charge associated to this symmetry transformation is represented 
by —^(Pn — qAfj,), given that is the (mechanical) momentum of the particle. All 
these considerations are valid, in particular, for the special choice £ = £m = dt- Finally 



2 Im 



(dl +qA)- f (dlo + qA) 



2-7T 

= —(w -«<!>), (3.12) 



which reproduces the standard result 



(w — q<E>) 

- em i abs ^ 



(3.13) 



Note that this derivation applies equally well to particles with non vanishing angular 
momentum, since adding a term d^I d<j> to the differential dl would not change the 
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imaginary part. Hence equations such as ea. (|3.6l) and ea. (l3.13|) are actually true for 
any angular momentum, though of course the implicit grey body factors do depend on 
it. 

The Kerr solution The observation that J d^I d<f> is real will not be true for ro- 
tating black holes. We can extend the calculation to the Kerr solution by noticing 
that throughout the complex manifold the azimuthal angle must also be rotated in 
order to keep the metric regular. More precisely, one has to shift <p — > <p — iVt H d/K, 
where £l H indicates the angular velocity of the horizon. Then, adding the term d^Idcf) 
to the differential dl would produce an imaginary term after a 7r-rotation, equal to 
—iirSl, H d^I I k. The outgoing trajectory from the past horizon is a classical solution so 
d^I = J, the conserved angular momentum. We thus obtain the result 

Im dl = Im f dl+ ^(w - fl H J) (3.14) 

or 

r em = r abse -^- n,,J > . (3.15) 

Of course in quantum theory the angular momentum is quantised. As is well known 
the emission and absorption probabilities for particles with energy u) and angular 
momentum j are related to the Bogoliubov /3-coefficients, whose computation is a 
classical problem involving only the relevant field equations. Unitarity in the space 
of classical solutions relates them to the transmission coefficient T u i m through the 
potential barrier around the horizon 

Tabs ± Tern = |^uiim| (3.16) 

where the plus (minus) sign is for fermions (bosons). Together with ea. (|3.15|) it gives 
the Fermi-Dirac (Bose-Einstein) spectrum. 

Equation (|3.6[) . or its cousins for the charged and Kerr black holes, turns out to 
be perfectly consistent with equations (|2.13[) and (|2.32l) previously obtained, but it 
has been derived under the unique assumption of the analyticity of the particle action 
in the complexified Schwarzschild/Reissner-Nordstrom/Kerr geometry, and nothing 
more. The simplicity of this result seems to us particularly important and probably 
deserves further discussion. In first place, it can be considered as the analogue of 
the well known "path integral derivation of black hole radiance" due to Hartle and 
Hawking [55], derived uniquely from the classical action I. In second place, it helps 
us to understand better how the Hamilton-Jacobi method works. At this regard 
we immediately see that reference to the Hamilton-Jacobi equation and insistence 
on considering only null trajectories (on which performing near-horizon expansions) 
are additional ingredients of the method which have nothing to do with the physical 
phenomenon of tunnelling. Even if they are not absolutely necessary from a conceptual 
point of view, their role turns out to be fundamental for computational purposes. As 
already said, most of the times the analytic continuation does not even exist, since 
only for special manifolds there is a complexified extension; while considering arbitrary 
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forbidden trajectories as ab in figH] would rule out any exact computation. 
As mentioned above, the regularisation scheme typically implemented in the tunnelling 
method is the Feynman's prescription, e.g. (|2.13D and (|2.32[) . At the beginning, we 
introduced it with no justification, either physical or mathematical, but now it is clear 
that Feynman's prescription 



where Q is any interval such that <E Q, has to be the preferred since it is the only one 
consistent with analytic continuation results, when available. From now on, Feynman's 
prescription will be adopted systematically and its justification given for granted even 
in those cases where there is no analytic proof of its validity. 

As a last remark, in the previous sections we repeatedly insisted on the importance of 
using well-behaved coordinates across the horizon. However, even this restriction seems 
to be less compelling in the light of analytic continuation. Equation (|3.4p , which is the 
best achievement of this section, expresses the remarkable fact that it is the difference 
of Im ^ dl and Im J* dl which provides the correct ratio between the emission and 
absorption amplitudes. As far as we are concerned with regular coordinates, only 
the integration over abc trajectories carries a non- vanishing imaginary part, being the 
absorption amplitude completely real. Insisting on singular gauges as Schwarzschild 
coordinates, however, nothing prevents the entering trajectory from picking up an 
imaginary contribution. It is exemplary, in this sense, the (Sg for) Schwarzschild 
gauge where both the escaping 



pick up an imaginary contribution providing in this way the correct result. This 
statement, if true for the Schwarzschild gauge, cannot be taken as a general rule 
valid for any singular gauge, as the reader may easily convince himself by taking into 
consideration isotropic coordinates. 

3.2 Equivalence of null geodesic and Hamilton-Jacobi methods 

At this point of the discussion, it seems natural to ask what relation exists between 
the aforementioned null geodesic and Hamilton-Jacobi methods. In literature, (see 
for instance [46]) we find many specific examples suggesting the essential equivalence 
between the two methods, but nothing resembling a mathematical proof. To this aim, 
notice that in the null geodesic approach one starts with 





(3.17) 



and the classically allowed ingoing trajectory 




(3.18) 




(3.19) 
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which looks non covariant. However we can use the full Liouville one-form w = p^dx* 1 
and write the more general expression, 



Im 



(3.20) 



which, without the " Im" , is nothing but the reduced action. We can perform the 
analytic continuation of this integral just as we did for the complete action, first by 
writing w = pijdU + pydV, then rotating ([/, V) from zero to 7r and finally integrating 
along the rotated curve. In this way, the imaginary part will be m{Upy — Vpu): but 
this is —iituj/k, where lo — ~p t is the Killing energy as measured at infinity. In all we 
get 

■KLO 



Im 



Im r vj H 



(3.21) 



As another example we consider the line element of the Kerr-Newman solution, 



ds 2 - 

ab KN — 



G 2 a 2 sin 2 i 



(r — r_)(r + r + ) 
(r 2 + a 2 )(r 2 + + a 2 ) 

G(r-r-) 
2k 2 _S 

-(2„2 ;„2 , 



s 2 , 



(r 2 +a 2 ) 2 (r 2 + +a 2 ) 2 
G 2 a 2 sin 2 9 (r + r + ) 2 



(U 2 dV 2 + V 2 dU 2 ) + 
dUdV + 



4k 2 + E 
Ga sin 2 9 
" K+ E(r 2 +a 2 ) 



{r\ + a 2 



(UdV-VdU) 2 + 



[£+(r - r_) + (r 2 + a 2 )(r + r+)] (UdV-VdU)dfa 



-Y,d9 2 



[(r 2 +a 2 ) 2 -a 2 sin 2 9 A] sin 2 



(3.22) 



where U, V are Kruskal-like null coordinates, defined as, 

\U\ = exp«(r* — t) , |V| = exp«(r* + 1) , 



(3.23) 



rfr* = (r 2 + a 2 )dr/A and the surface gravity k can refer to n + or k_ according to the 
space-time region of interest. The co-rotating angular coordinate 4>+ is defined as 



(3.24) 



with £l H the horizon angular velocity and G :— (r—r + )/UV. For all other symbols, the 
reader is referred to Section [2.61 where a complete treatment of axis-symmetric space- 
times has been given. This metric is particularly important since in the appropriate 
limits it comprises the maximal analytic extension of Kerr-Newman, Kerr, Reissner- 
Nordstrom and Schwarzschild black holes [52"] . 

The analytic continuation of the coordinates in the Euclidean section requires 



t -tt — i §/k 



and 



4>^ cf)-in H tf/K 



■d E [0,tt] 



(3.25) 
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where from now on, we consider k = k+. We have 137 = pijdU + pydV + p^d4> and 
therefore 

lmj^w= Im / to - ^—P<P + k(U Pu - Vp v ) ■ (3.26) 

But Upjj — Vpv = —K pt = k cj is the energy and p^, = J the conserved angular 
momentum. We end with our main formula, Equation (|3. 41) . 

2 Im Q|£ w - ^ w \ =— (uj-n H J) (3.27) 

which proves the equivalence of a generalised null geodesic method with Hamilton- 
Jacobi ansatz for all stationary, axis-symmetric, charged black holes. In retrospect, 
this could not be otherwise because the Hawking effect is an energy conserving process, 
so that the reduced action is all one needs in a stationary geometry. 



3.3 Alternative proposals 

The tunnelling method has been explored by many authors and counts a considerable 
literature which is all but monolithic in presenting the subject. The large number 
of viewpoints is probably the best indication of the non-trivial nature of the topic. 
Conversely, it has been also the origin for a certain confusion |120j which made difficult 
to export its ideas outside the community of experts in the field. 

In this sense, it is paradigmatic to consider the so called "factor-two problem" |121) , 
actually a fake problem consisting in the fact that a naive application of the tunnelling 
method may end with a black hole temperature which is twice the expected result. 
Roughly speaking, the troubles typically arise when one tries to adopt singular coor- 
dinates (like Schwarzschild) and uncritically assumes that the absorption amplitude is 
real. Then, as it is clear from Equation (|3.17p . one ends with a temperature equal to 
l/47rM. As we have seen in Equation (|3.18[) , this is simply due to the fact that we are 
missing an equal and opposite contribution coming from the absorbed particle which, 
entered the master Equation p.4[) , provides indeed the correct temperature. 
Mitra, for example, circumvents this bug in the following way |122j . The Hamilton- 
Jacobi equation is a differential equation for the action of both the ingoing and 
outgoing particle. Because of the staticity and spherical symmetry, we can solve 
the action in terms of the standard ansatz for the separation of variables to get, 
Iout,in = —u>t ± iirujrg + . . . It is worth to stress that in general, from reconstruc- 
tion of the action can arise some constant of integration C. On the other hand, it was 
argued, the classical theory of black holes tells us that an incoming particle is absorbed 
with probability equal to one. In order to match these apparently contradictory facts, 
we use the freedom of choosing the integration constant C to impose the classical con- 
straint on the absorption probability, that is C — irujr g , so that Im I out = 2irr g uj and 
T em (x exp(— 8ttMuj), as expected. 

However, if we try to repeat the same argument for the Schwarzschild geometry in 
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isotropic coordinates, 



+ p 2 (d0 2 + sin 2 0d(j> 2 )} , (3.28) 



where this time the horizon is located at the minimal two-sphere po — M/2, the in-out 
particle action derived from the Hamilton-Jacobi equation reads 

I in = (-cot + C) - iAttMuj 

I out = (-at + C)+ iAnMuj , (3.29) 

and the classical constraint lmli n — fixes ImC = AttMuj implying a Hawking tem- 
perature of l/167rM which, this time, is a half of the correct, well-known, result. 

A relevant achievement in the discussion of the "factor-two problem" has been pur- 
sued by Angheben et al. [91] and Stotyn et al. (1 23] . who showed that correctly using 
the theory of distributions in curved static space-time it is possible to recover the cor- 
rect result for the Hawking temperature even working in singular gauges. In fact, the 
non locally integrable function 1/r does not lead to a covariant distribution l/(r±?e); 
where instead, Feynman's regularisation of an invariantly defined distance makes un- 
necessary to work with the complex action of ingoing particles. The natural choice 
in substitution of 1/r is represented by the proper spatial distance, as defined by the 
spatial metric, da 2 = g mn dx m dx n (m,n = 1,2,3), a quantity that is invariant under 
re-definitions of the spatial coordinate system and, more generally, with respect to the 
subgroup of the gauge group consisting of time re-calibrations and spatial diffeomor- 
phisms, defined by 

t' = t'(t,x l ) 
x' j (x l ) . 

We consider the previously discussed spherically symmetric metric in diagonal gauge, 

dv 2 

ds 2 = -V(r)dt 2 + — — — - + r 2 (d6 2 + sin 2 9d<j> 2 ) , (3.31) 
W[r) 

where we suppose V(r) and W(r) to have both a simple pole in correspondence of the 
horizon th, 

V(r) =V'(r)(r -r H ) + ... , W(r) = W'(r)(r - r H ) + • ■ • ■ (3.32) 
The proper distance reads, 

dr 2 



/./ _ ii . ... , (3.30) 



y/r - r H , (3.33) 



so the classical action of the outgoing massless particle is 

lout = -urt +^Ldr ; = = -ut + 2 to — . (3.34) 

J y/V(r)W{r) y/V'(r H )W'(r H ) J <* 
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The integral is still divergent as soon as a — > 0, namely the horizon is reached, but 
now the prescription corresponding to the Feynman's propagator 1/a — > l/(a — ie) 
selects the correct imaginary part for the classical action, 

lnio i 

I out = — ; + (real contribution) . (3.35) 

VV'(r H )W'(r H ) 

From comparison of (|3.35[) with (|2.24[) . we soon read off the horizon surface gravity k 
for the dirty black hole (|3.3ip . 



y/V'(r H )W'(r H ) 
k = , (3.36) 

which is amazingly in agreement with the result obtained through the conical singu- 



larity method [1241 . From a qualitative point of view, this coincidence can be ascribed 
to the fact that both the conical singularity and the proper distance methods refer 
the energy uj to the same Killing observer. However, we notice that, without any fur- 
ther assumption, a dirty black hole space-time is not asymptotically flat, so it is not 
completely clear how to normalise in general the Killing vector £^ . 

Can we say that the factor- two problem has been solved? We have seen that, at least 
confining to the evaluation of the Hawking temperature of a Schwarzschild black hole 
in Schwarzschild coordinates, there are two satisfactory ways to circumvent the factor- 
two problem arisen because of the use of singular charts: the former, essentially based 
on the classical requirement that the absorption probability be unity; the latter, related 
to an appropriate way of dealing with distributions in curved space-times. However, 
the combination of these techniques as made in [123] has the deleterious drawback of 
making the temperature again a half of the correct result. According to these authors, 
the point is that to apply the Feynman prescription, the proper distance should be 
extended to negative values, which seems geometrically impossible. The authors then 
note that to recover the standard result Tjj = l/8irM, one ought to use an alternative 
contour of integration when dealing with proper distance formulation, J da / 'a, namely 
the quarter-contour of figEl (left) rather than the standard semi-circular contour of 
figH (right). 



Figure 3: Standard semi-circular contour of integration in the complex plane of J da /a 
and the quarter-contour proposed in [123 . 
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Another point of view is the one advocated in (1 25] where the authors are among 
the first to point out the possible role of a "temporal contribution" to solve the factor- 
two problem, at least in the context of singular gauges. In fact, temporal and spatial 
contributions are on equal footing in the Hamilton-Jacobi method as formulated in 
Section 12.21 since it is covariant. For example, the master equation (|3.4p does not 
distinguish at all between "time" or "space" coordinates, something which makes per- 
fectly sense according to the principles of general relativity. On the other hand, the 
null geodesic method as originally formulated does not involve a time integration, but 
appears non covariant and also non canonically invariant. Along the line of thought 
initially formulated by [1261 , Singleton el al. |125j make use of closed integration paths, 
necessary to ensure invariance of the emission rate under canonical transformations. 
To all the extents, the integration over the closed circuit § p r dr corresponds to tak- 
ing into consideration the imaginary contribution related to the action of both the 
incoming and outgoing particle, exactly as in p,18p : <f — The null geodesic 
method as formulated in Section [3~!?l evidently, is covariant as well as invariant under 
canonical transformations in phase space. Hence, even acknowledging a certain loss 
of importance of canonical transformations in the passage from classical to quantum 
theories, this is certainly an interesting point. On the other hand, since not all the 
canonical transformations correspond to unitary transformations in the quantum the- 
ory, perhaps a more urgent issue is not the canonical invariance of the decay rate, but 
the covariance of the result - as demanded by general relativity - and the unitarity 
of the process - as demanded by quantum mechanics - something that nobody knows 
how to solve. 

As stressed above, the Hamilton-Jacobi method was born also, surely not only, by 
the desire of disentangling Hawking radiation from back-reaction, in the sense that it is 
the existence of Hawking radiation that founds the problem of back-reaction, not vice 
versa. Given that the Hamilton-Jacobi method can prove black hole radiance without 
taking into account the back-reaction, it is natural to ask how the back-reaction can be 
incorporated in this method. As shown by Medved and Vagenas (12711128] . a tunnelled 
particle of instantaneous energy ui effectively "sees" a space-time metric of the form 
((3.31[) in which the horizon radius rjj = ru{M — uS), M being the original Misner- 
Sharp black hole mass. Because of quantum uncertainty principle, it seems too crude 
the approximation of a discontinuous jump from M to M — uj for the black hole mass. 
Rather, the particle action corrected by energy conservation, turns out to be given by 

m \\j . (3.37) 
Jo V V h( M ~ u'W H (M - uj') 

where the subscript H is a reminder that we are evaluating this quantity at the black 
hole horizon. The quantum corrected tunnelling probability deduced from (j3.37|) turns 
out to be that expected from statistical mechanics and in perfect agreement with the 
null geodesic method of Section [2. 1[ namely T em ~ exp ASbh- A somewhat different 
treatment of this problem with equivalent results, within the Hamilton-Jacobi method, 
can also be found in [129] . 

In this section we have taken into consideration alternative proposals to the null 
geodesic and Hamilton-Jacobi methods outlined in Section 12.11 and 12.21 Without 



I' out = I dJl n 



M-s-M-t, 
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claim for completeness, we have tried to point out their salient features, referring to 
original papers for details. We cannot deny that all of them contain small pieces of 
the complete picture, but at the same time the persistent lack of manifest general 
covariance is to be considered a defect, demanding a suitable reformulation. Exactly 
this was our desire in writing this section. 

4 Quantum Tunnelling: Dynamical Black Holes 

In the previous sections we saw that further insight into the Hawking's effect can be 
obtained by the use of the tunnelling methods. We showed that in some circumstances 
it is even possible to recover exact results by analytically continuing the integral of 
the classical action from a tunnelling path to a classical path, and how this possibility 
provided a foundation to the tunnelling method. But of course everything is strictly 
true only for stationary black holes. 

In this section we shall extend this work to what we shall aptly name spherically 
symmetric dynamical black holes, epitomised by the Vaidya solution |130] and its gen- 
eralisations by Bardeen |131j and York (132) . As a matter of fact, it has long been felt 
that the usual semi-classical treatment of stationary black holes had to be extended 
to cover at least slowly changing black holes. By this expression, we mean black holes 
that can be still described in terms of few multipole moments such as mass, angu- 
lar momentum and the charges associated to local gauge symmetries, except that the 
parameters and the causal structure are allowed to change with time. Although a tech- 
nical definition of a "slowly varying black hole" can be given in some cases, an example 
being the Booth-Fairhurst slowly evolving horizon |133j . in general this depends on 
the actual physical processes involved. For example, in the case of Hawking's evapora- 
tion, conditions for slowness in the presence of a near-horizon viscous fluid have been 
given by Brevik |134) in an interesting attempt to generalise 't Hooft's model of the 
self-screening Hawking atmosphere (quantum corrections to this model can be found 
in [135] ). In general, it is understood that the black hole temperature has to be much 
smaller than the Planck mass, while in order to study the effects of the expansion, the 
Hubble rate H" 1 should be comparable with the black hole emission/absorption rates. 

Now, even for the case of a slow evolution, it was pointed out by Fredenhagen and 
Haag long ago [136] that by letting the mass of the black hole to change with time, 
the radiation would originate from the surface of the black hole at all times after its 
formation. This poses the question: what and where is the surface of a dynamical 
black hole? The issue baffled scientists since the beginning of black hole physics and 
produced several reactions during the Nineties, which eventually culminated with the 
notion of outer trapping horizons by Hayward 137H139] and the isolated and dynamical 
horizons of Ashtekar and co-workers [140H143] (an extensive review is in |144j ). 
One is concerned to study at first the dynamical version of the event horizon and to 
provide a mathematical definition which is able to capture a useful local notion of it, 
encompassing the time-lasting textbook definition given in Hawking-EUis renowned 
book [55] . Armed with a precise notion of horizon, we shall proceed to study the 
instabilities occurring near the horizon of the changing black hole. This question looks 
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non trivial since a changing horizon is typically embedded in a dynamical space-time 
and it is not even expected to be a null hyper-surface, though it is still one of infinite 
red-shift. Thus we shall start by reviewing what has been done on this question for 
spherically symmetric metrics. We can anticipate that the arrival point will be the 
very important local notion of future trapping horizons and their associated surface 
gravity, as defined by Hayward. 

4.1 Horizons and surface gravity 

Not long after the classical definition of the event (EH) and the apparent horizons 
(AH) (boundaries of trapped 3-dimensional space-like regions within partial Cauchy 
surfaces), several quasi- local notions of dynamical horizons were proposed in the lit- 
erature (another nice review is in |145j ). perhaps starting with the notion of perfect 
horizon due to Hajicek [146 . But this only applied to equilibrium black holes while the 
apparent horizon, being tied to a partial Cauchy surface, only represents a localisation 
in time. Moreover, it has proven not possible to formulate thermodynamic laws for 
AH similar to those holding good for the event horizons. 

The first successful attempt to go beyond the limitations imposed either by the 
instantaneous character of apparent horizons or by the global, teleological nature of 
event horizons is due to Hayward. His concept of a future outer trapping horizon (to 
be abbreviated as FOTH) then evolved either into less constrained definitions, like the 
Ashtckar-Krishnan dynamical horizons (DH), or more specialised ones, like the Booth- 
Fairhurst slowly evolving FOTH. Although the horizon as defined by Hayward will be 
central in this review, for sake of completeness we shall give an updated list of locally 
or quasi-locally defined horizons which appeared over the years, each playing some role 
in the problem of understanding dynamical black hole (more precise definitions will 
be given soon). At least four types of horizon have been defined over the years: 

1. non expanding and perfect horizons (Hajicek 146 ); 

2. trapping horizons (Hayward [13711139] ): 

3. dynamical horizons (Ashtekar and Krishnan 142,143 ); 

4. isolated and weakly isolated horizons (Ashtekar et al. [140, 141 ); 

5. slowly evolving horizons (Booth and Fairhurst [133] 1 . 

Most of these newly defined horizons have very desirable properties: they do not 
require a space-like hypersurface, no notion of interior and exterior and no conditions 
referring to infinity, like asymptotic flatness for example (all are non local conditions). 
Moreover, they are not endowed with teleological features (they do not anticipate the 
future, so to speak) and, given a solution of Einstein equations, one can find whether 
they exist or do not exist by purely local computations. Finally, unlike EH they are 
related to regions endowed by strong gravitational fields and are typically absent in 
weak field regions. 
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We recall that the expansion 9 of a bundle of null rays is the rate of change of area 
transverse to the bundle 



where v is a parameter along the rays. All quasi-local horizons rely on the local concept 
of trapped or marginally trapped surface: this is a space-like closed two-manifold S 
such that 9+9- > 0; if £± are the future-directed null normals to S, normalised to 
9nu^+i- = £+■£-= —1, then 9 + , 6L are the respective expansions or optical scalars, 
that is the expansions of the two bundles of null rays orthogonal to S. It is further 
assumed that £+ is associated to an outgoing null geodesic beam, so that in a region 
of not too strong gravity 9+ > and the beam is expanding, as for example within 
the exterior of a black hole at a safe distance. It follows then that £- is associated to 
an ingoing null geodesic beam, with < and the beam contracting along the way. 
It is always possible to choose double null coordinates such that 



where r is the areal radius, defined so that a metric sphere has area A = Airr 2 . To 
cover black holes rather than white holes it is further assumed that both expansions 
are negative (or non positive) on a trapped (marginally trapped) surface. If 9 + 9- < 
the surface is untrapped and marginal if 9+9- — 0. In spherical symmetry with radial 
coordinate r this means the co-vector dr is temporal, spatial or null respectively. A 
further subdivision may be made: a trapped surface is future if 9± < and past if 9± > 
0. A marginal S with 9 + = is future if 9- < 0, past if 9- > 0, bifurcating if 9- = 0, 
outer if d-9+ < 0, inner if d-9 + > and degenerate if cL0+ = [T5^IT3WTT7l[ng] . 
Equivalently, the vector field g r ^dfj, is future causal (it means non-space-like) or past 
causal, respectively. For a better understanding of these terms we may perhaps note 
that 



So, for example, on a future trapped surface with 9± < the radial increasing vector 
is past directed, confirming the intuition that one cannot escape the trapped surface 
without moving in the past. 

So much far for black holes, using the mental picture that we, the observers, are 
outside it at some radius larger than the radius of the trapped region. For white holes 
one encounters also marginally trapped surfaces with 9- — and 9 + > 0. In cosmology 
one may encounter again both possibilities, marginally trapped surfaces with either 
9- — or 9 + = 0, but we defer their description to Section [5] 

The abstract definitions will become much more simple and intuitive when referred 
to spherically symmetric space-times in Section 14.21 where only few of them will be 
actually used. Here we take the opportunity for few more definitions. The optical 
scalars can also be defined as follows: the induced metric on each S is 




6± = - d±r 



(4.1) 



r 




(4.2) 
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for in fact q^i^ = 0. Let q^ = g^ + t+£ v _ + not the inverse of q^. Then 

q^ v is the projection tensor to T»(S), the tangent space to S. Associated to the 
null vector fields £± are the projected tensor fields (0±) lll/ = (ffqfflal±p and their 
decomposition into symmetric, anti-symmetric and trace part. They are tensors on S 
because (0±) tlI y£± — (0±) liu ^ = 0. The twists (anti-symmetric parts) vanish since the 
geodesic beams are normal to S (we have to assume that S is part of a continuous 
family of surfaces St ) . The expansions are then given by the traces 

0+ = q^n\, 6- = q^Vr . (4.3) 

Finally, the shear is the trace-free symmetric part (note that g^q^v = q^q^v = 2) 

off = ±I1V - \q^± ■ (4.4) 

As before, we indicate the Lie-derivative along a vector field X by the calligraphic 
symbol Cx and C± = £e ± . Thus, for example, C-9^ u = l a Ld fi ilv -\-d il ,l a iQ av -\-d v l t.9 ll0l , 
while on scalars it acts as an ordinary partial derivative. Let us describe the listed 
horizons in turn, adding comments where appropriate. A black triangle down T will 
close the definitions. 

Perfect and non-expanding horizons — A perfect horizon is a smooth three-dimensional 
null sub-manifold H of space-time with null normal such that its expansion 0£ = 
on H and which intersect space-like hypersurfaces in compact sets. T 

If in the last clause H is topologically K x S 2 and moreover the stress tensor 
is such that —Tjfl v is future causal for any future directed null normal V*, then H is 
called a non-expanding horizon .T 

All stationary horizons are perfect, but the converse is not true. These horizons are 
not very relevant in the present context since they apply only to equilibrium black 
holes. 

Future outer trapping horizons — A future outer trapping horizon (FOTH) is a 
smooth three-dimensional sub-manifold H of space-time which is foliated by closed 
space-like two-manifolds S t , t € R, with future-directed null normals £± such that: 
(i) the expansion 9 + of the null normal £+ vanishes; (ii) the expansion 9- of is 
negative; (iii) £_9+ < 0. A future inner trapping horizon (FITH) is as above except 
that on it £ + 9_ > 0. ▼ 

This proved to be the most important definition since forms the basis for almost all 
other definitions. Condition (i) requires strong fields since certainly 9+ > in weak 
fields. The condition (ii) is related to the idea that H is of the future type (e. g. a 
black hole rather than a white hole); (iii) says that H is of the outer typqll, since 
a motion of S t along I- makes it trapped. It also distinguishes black hole horizons 
from cosmological ones, which for an expanding universe are FITH. Note that no 
reference is made to any space-like hypersurface, nor to infinity. Trapping horizons 
are locally defined and have physical properties such as mass, angular momentum 
and surface gravity, satisfying conservation laws 142, 149 . They are a geometrically 
natural generalisation of Killing horizons, which are stationary trapping horizons. A 

7 For example the Cauchy horizon in the Reissner-Nordstrom solution is of inner type. 
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non-stationary trapping horizon is not null, but still has infinite red-shift. Unlike event 
and apparent horizon they do not require asymptotic flatness. 

One can always find a scalar field C on H so that the vector fields 

V = i\ - Ct_ and W = l\ + Gt , (4.5) 

are respectively tangent and normal to the trapping horizon. Note that V • V = 
—N-N = 2C. Hayward 137,139 showed that if the null energy condition (abbr. NEC) 
holds then C > on a FOTH. Thus, the horizon must be either space-like or null, and 
it is null if and only if the shear aj^y and T^ v ^i u + both vanish across H. Intuitively, 
H is space-like in the dynamical regime where gravitational radiation and matter are 
pouring into it and is null when it reaches equilibrium. Conversely, on a FITH is 
either null or time-like. 

It is worth mentioning that the second law of the mechanics of trapping horizons 
follows quite easily from this apparatus. Taking the Lie-derivative of y/q, which is the 
area density corresponding to the metric q^ on the cross-sections of H, we get 

C v ^q=-C^e^ . (4.6) 

By definition is negative on H and we have just seen that, barring violations 
of the null energy condition, C also is non-negative. Since V is future directed we 
obtain the local form of the second law: If the null energy condition holds, then the 
area element ^Jq on a FOTH is non- decreasing along future directions. 
Integrating over St the same law applies to the total area of the trapped sections. As 
long as the null energy condition is maintained it will be non-decreasing, reaching a 
constant value if and only if the horizon becomes a null hypersurface. 

The main difference between an apparent horizon as defined in Hawking-Ellis and 
the trapping horizon of Hayward is that the AH represents the instantaneous surface 
of a black hole, i.e. it needs a (partial) Cauchy surface E and it is very sensitive to the 
choice of E. To compute the AH one needs only its metric and the second fundamental 
form of E, namely the initial data for Einstein's equations. Hayward's horizon instead 
is a null hypersurface H which is insensitive to a choice of E and does not refer to 
spatial infinity. The trapping horizon is a foliation of H. 

The causal character as well as the area law required the validity of the NEC; an 
evaporating black hole violates NEC, therefore the area law will also be violated and 
the horizon will be time-like. How it could be that tunnelling along a classically 
forbidden path is still possible will be seen soon and represents the real possibility of 
the radiation process even in the temporary absence of a global event horizon. 

Next come Ashtekhar and co-workers: they observe that key results, such as the 
area increase, do not depend on the sign of C-9 + . Hence the following weaker notion 
was introduced: 

Dynamical horizons — A smooth three-dimensional, space-like sub-manifold H of 
space-time is a dynamical horizon (DH) if it can be foliated by closed space-like two- 
manifolds St, with future-directed null normals £± such that: (i) on each leaf the 
expansion 9 + of one null normal i + vanishes; (ii) the expansion of the other null 
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normal is negative. A submanifold H such that 9+ = and 6- < is also called a 
marginally trapped tube (abbr. MTT) T 

Like FOTHs, a DH is a space-time notion defined quasi-locally, it is not related 
to a space-like hypersurface, it does not refer to infinity, it is not teleological. A 
space-like FOTH is a DH on which C-6+ < 0; a DH, which is also a FOTH, will 
be called a space-like future outer horizon (SFOTH). The precise properties of such 
horizons are fully discussed in [142, 143 . Suitable analogues of the laws of black hole 
mechanics hold for both FOTHs and DHs. We only note that DHs cannot describe 
equilibrium black holes since they are space-like by definition. Likewise, DHs cannot 
describe evaporating black holes since for them the trapping horizon is time-like, but 
are better suited to describe how a black hole grows in general relativity. In general, if 
a MTT is space-like it is a DH. If a MTT is time-like it is a time-like membrane. It is of 
interest that in the Oppenheimer-Volkoff dust collapse the unique MTT on which each 
marginally trapped surface (MTS) is spherical is time-like [150 ! This shows that even 
in spherical symmetry the existence of DH with spherical sections is far from obvious. 
However, for perfect fluid collapse and spherical scalar field collapse the MTTs arc 
space- like. 

If X, Y are tangent to a non-expanding horizon we can decompose the covariant 
derivative 

V X Y = D X Y + N(X, Y)l + + L{X, Y)t- 

wherc Dx is the projection of the vector XI xY onto the spheres St in H. If X is 
tangent to the spheres then Dx is the covariant derivative of the induced metric q^u, 
and if X is tangent to H one may regard the operator V x = Dx + N(X, ■)£+, acting 
on vector fields, as a connection on H . If this connection is "time independent" then 
the geometry of H is time independent too and we have Ashtekar et al. notion of a 
horizon in isolation. 

Isolated horizons — A non-expanding horizon with null normal such that [Ci, Vjt] = 
along H.J 

These horizons were intended to model black holes that are themselves in equi- 
librium but possibly in a dynamical space-time. For a detailed description of their 
mathematical properties we refer the readers to Ashtekar-Krishnan's review [144] . 

Slowly evolving horizons — No matter or radiation can cross an isolated horizon, so 
the meaning of the first law for them cannot be treated in full generality. Booth and 
Fairhurst [133] established this law for slowly evolving FOTHs as defined by Hayward 
by introducing dynamical notions of surface gravity and angular momentum. For this 
purpose, the concept of a slowly evolving FOTH was then defined. 

All the horizons just introduced have their own dynamics governed by Einstein 
equations. There are for them existence and uniqueness theorems |151) . formulation 
of the first and the second laws |137U144lll52) and even a "membrane paradigm" 
analogy. In particular, they carry a momentum density which obey a Navier-Stokes- 
like equation generalising the classical Damour's equations of EHs, except that the 
bulk viscosity (foth = 1/16-7T > [15311154] . The newly introduced horizons are also 
mostly space-like or null, therefore the role they may play in the problem of black hole 
quantum evaporation is unclear. In this connection the following notion can be useful. 
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Time-like dynamical horizon — A smooth three-dimensional, time-like submanifold 
H of space-time is a time-like dynamical horizon (TDH) if it can be foliated by closed 
space- like two-manifolds St, with future-directed null normals t± such that: (i) on 
each leaf the expansion # + of one null normal 1+ vanishes, (ii) the expansion 9- of the 
other null normal is strictly negative. T 

It must be recalled that a FOTH for which the NEC does not hold can be time-like 
as well. What about non spherical MTS? The Vaidya metric does not admit other, non 
spherical horizons which also asymptote to the non expanding one. This is general, 
so a DH is unique, that is a space-like 3-manifold cannot be foliated by two distinct 
families of MTS each endowed with a structure of a DH. This should also hold for a 
TH but we do not know whether there is a proof. According to Ashtekhar et al. a 
better control of uniqueness is the most important open issue of the classical theory 
of dynamical horizons, be trapped or in isolation. 

Surface gravity The surface gravity associated to an event horizon is a well known 
concept in black hole physics whose importance can be hardly overestimated. Surpris- 
ingly, a number of inequivalent definitions beyond the standard one appeared over the 
last 15 years or so in the field, with various underlying motivations. We have collected 
the following (we rely on the review of Nielsen and Yoon [155] ) : 

1. the standard Killing surface gravity (Bardeen et al. [156] . textbooks); 

2. a first definition given by Hayward in [137 ; 

3. the effective surface gravity appearing in Ashtekar-Krishnan (144) : 

4. the Fodor et al. definition for dynamical spherically symmetric space-times |157) : 

5. the Visser [TJ] and Nielsen- Visser surface gravity; 

6. still one more definition by Hayward 139 , using Kodama theory of spherically 
symmetric space-times [158] . 

In addition there are some more technical definitions due to Mukohyama and Hayward 
|159j and to Booth and Fairhurst, the latter related to their notion of evolving horizons 
133 . Except for the last item, which is what the tunnelling method leads to, the 
remaining definitions will be less relevant in this review, so we reserve to them only 
some brief considerations. 

The Killing surface gravity is related to the fact that the integral curves of a Killing 
vector are not affinely parametrized geodesies on the Killing horizon H, where the 
norm £ 2 = 0. Hence 

defines the Killing surface gravity k on H, where = means evaluation on the horizon. 
The Killing field is supposed to be normalised at infinity by £ 2 = — 1. The definition 
can be extended to EHs that are not Killing horizons, by replacing £ with the null 
generator of the horizon. However there is no preferred normalisation in this case, and 
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this is one reason of the debating question regarding the value of the surface gravity 
in dynamical situations. 

Hayward's first definition was motivated by the desire to get a proof of the first law 
for THs. It is defined for a future trapping horizon without appeal to un-affinity of 
null geodesies, as 

^y/-n»Vp0 e . (4.7) 

This quantity is independent on the parametrisation of the integral curves of the vector 
field since the evaluation is on a marginal outer surface where n-£ = — 1 and 9f = 0. 

Given a weakly isolated horizon H, Ashtekar and Krishnan showed that for any 
vector field t a along H with respect to which energy fluxes across H are defined, there 
is an area balance law that takes the form 

5E* = t^-^,5A s + work terms 
with an effective surface gravity given by 

1 dr 
K ~ 2RdR ■ 

R is the areal radius of the marginally trapped surfaces, i.e. As = AirR 2 , the function 
r is related to a choice of a lapse function and finally E t is the energy associated with 
the evolution vector field t a . For a spherically symmetric DH a natural choice would 
be r = R so R = 1/2R, just the result for a Schwarzschild black hole. To illustrate the 
naturalness of this definition, consider a slowly changing spherically symmetric black 
hole with mass M(v), where v is a time coordinate. Defining the horizon radius at 
each time by R = 2M(v) and As — AnR 2 , we can differentiate M 

2 2R 8tt 

so as to obtain 8M = k5As/8it. One recognises the usual area law in differential 
form with surface gravity R = 1/2R = 1/4M. Consider, however, the more general 
possibility where the horizon is the solution of the implicit equation R = 2M(v, R), as 
it happens for example in the Bardeen-Vaidya metric. The same computation leads 
to M = A s (1 - 2M')/(16tt j R), leading to 

K ~ 7T7 0- ~ 2M') (4.8) 
AM 

a prime denoting the radial derivative. The surface gravity here deduced does not con- 
form to Ashtekar et al. definitions, suggesting that its value depends on the definition 
of the black hole mass one is adopting. 

The definition of Fodor et al. looks like the Killing form of the surface gravity in 
that kI v = i^V^i", where now is an outgoing null vector orthogonal to a trapped 
or marginally trapped surface. This is because, as a rule, such null vectors are not 
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afBnely parametrized, although they can always be parametrized so that k — 0. Fodor 
et al. choose to fix the parametrization so that 

k = -u^fV^ 

with n M affinely parametrized and normalised to n- 1 = —1 at space-like infinity, where 
t a is the asymptotic Killing field. Note that this definition is non local but looks like 
a natural generalisation of the Killing surface gravity. 

As far as we know, the Visser and Visser-Nielsen surface gravity is only defined 
in dynamical Painleve-Gullstrand coordinates, so there is no guarantee that it is a 
geometrical invariant. Indeed, we shall see that for this class of metrics it is different 
from the invariant surface gravity computed within the general Hayward formalism. 

The last item is a local geometrical definition of the surface gravity for the trapping 
horizon of a spherically symmetric black hole [139] . so we postpone a full discussion 
to Section 14.2.41 Basically, one introduces local null coordinates in a tubular 
neighbourhood of a FOTH and for any sphere of radius r one defines the quantity 

« = \ (g + -d-e + + l - g +-e + e^j . (4.9) 

Evaluated on a trapping horizon, 9 + = 0, it will be positive precisely when the horizon 
is of outer type (cL# + < 0, recall that with our conventions g + ~ < 0) and k = if 
degenerate. The definition may look somewhat artificial, but in fact it can be put in 
a form that strongly resembles the Killing surface gravity of stationary black holes. 
To see this, we anticipate a result of the next section according to which, following 
Kodama [158 , any spherically symmetric metric admits a unique (up to normalisation) 
vector field such that V(K^ G^ u ) — 0, where G^ v is the Einstein tensor; for 
instance, using the double-null form, one finds 

K = -g+-(d + rd--d-rd+). (4.10) 

The defining property of K shows that it represents a natural generalisation of the time 
translation Killing field of a static black hole. Now consider the expression K^V^K^y. 
it is not hard to see that on H it is proportional to K v . So one defines the dy- 
namical surface gravity as K^V^K^ = —nK u . For a Killing vector field V U K^ is 
anti-symmetric so the definition reduces to the usual one. 

4.2 Spherically symmetric fields 

In this section we discuss the general time-dependent spherically symmetric metric 
that will be the arena of our dynamical tunnelling computations. We have a twofold 
intent in doing this. If Hawking radiation proceeds by emission of discrete energy 
quanta, a continuous description of the changing metric would only be possible when 
their number is so large so as to simulate a continuous streaming of energy. In such 
a case, a dynamical excreting black hole may be modelled by a continuously diffcrcn- 
tiable solution of Einstein's equations for most of its history. Given that, it would be 
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quite ironic if there were not evidence of a continuously operating quantum emission 
process. In other words, if the tunnelling method only worked for the event horizons of 
stationary black holes then there would be a problem with the use of such dynamical 
metrics as models of black hole evaporation. 

The second aspect has to do with the region where the radiation originates. The 
tunnelling calculation suggests the the outgoing radiation is emitted from the trapping 
horizon, not the global event horizon. And it further suggests that the semi-classical 
probability is related to the horizon surface gravity as defined in eq. (|4.9l 139 '.. Thus 
we have first of all to learn how to compute these things in a spherically symmetric, 
time dependent metric. 



4.2.1 Forms of the metric 

The paradigm of a spherically symmetric metric is of course the Schwarzschild vac- 
uum solution of Equation (|2.1[) . To introduce dynamics one could think to make the 
replacement M — > M(t); however the resulting metric has a curvature singularity at 
r = 2M(t); for instance, the scalar curvature is 

_ 2r(r - 2M)M + 4M 2 
~ (r - 2Mf 

so only for special values of the mass function is 1Z — 0; other invariants will diverge 
though. A more general "Schwarzschild gauge" can be written down 

ds 2 = (i - ™M) dt 2 + (i - 2 -MM) 1 dr * + Ade * + sin2 m 



r 

(4.11) 

which was studied in depth by Visser in his classic "Dirty Black Holes" [B7j . We shall 
not make use of his metric in the following. 

It was discovered by Vaidya [160] that by taking the Schwarzschild metric in ad- 
vanced Eddington-Finkelstein coordinates (v, r, 9, <j>) (or retarded, (it, r, 8, <f>) for other 
purposes) 

ds 2 = - (l - dv 2 + 2dvdr + r 2 {d9 2 + sin 2 0d(j>) (4.12) 

and by making the replacement M — ¥ M(v), one obtains a regular solution of Einstein's 
equations except at the origin r = 0, the shining star solution of Vaidya. The metric is 
still not the most general one allowed by the symmetry: in fact a spherically symmetric 
metric can depend at most on two arbitrary functions of the coordinates (v,r). One 
form is due to Bardeen |131j 



ds z = -e zq> I 1 - — dv 1 + 2e*dvdr + r z (d6 z + sin 2 6d(j)) (4.13) 

where $ and M are functions of v, r. Another important form is obtained by passing 
from Schwarzschild time to the proper time of a radially infalling observer 

t v =t + 2V2Mr + AM log ( - V^ W\ 

& \^+V2Mj 
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The metric is the Painleve-Gullstrand metric we encountered in Section [2] 



ds 2 = - (l - dt 2 p + dt p dr + dr 2 + r 2 [d9 2 + sin 2 Od(f>) ■ (4.15) 

Taking now M = M(t p ) we get a dynamical black hole long studied by Lindesay and 
Brown [161H163] and others |164| . There is a more general form whereby M can also 
depend on r and a second function $(i,r) is introduced, say 



dg 2 = _ e 2$( t » ^ _ d p + 2e *(L^^M didr + dr 2 + r 2 dn 2 5 (4J6) 

where dfl 2 := d6 2 +sin 2 9d(f>. This metric is locally diffeomorphic to the metric ea. (|4.1ip 
under a change of time t — > t(t, r), with the mass treated as a scalar field, i.e. M(t, r) = 
M(t(i,r),r). It was thoroughly studied by Visser and Nielsen in In this rather 
impressive work the dynamics and decay of evolving horizons were investigated mainly 
using a dynamical version of Painleve-Gullstrand coordinates. The metric ea. (l4.16p is 
also locally diffeomorphic to the metric ea. (|4.13p . under a change v — > t(v,r), where 
now M (t, r) = M (v(t, r), r). 

All these spherically symmetric metrics (though not only these) are special cases of 
metrics that can locally be expressed in the warped form 

ds 2 = y, ul.r'd.r' + R 2 (x l )dQ 2 , i,j € {0, 1} , (4.17) 

where the two-dimensional metric 

dy 2 =y lJ (x i )dx i dx j (4.18) 

is referred to as the normal metric (in the space normal to the sphere of symmetry) , 
x l are associated coordinates and R(x l ) is the areal radius, considered as a scalar field 
in the normal two-dimensional space. The two-dimensional metric y-y has only one 
degree of freedom and the function R{x l ) is the second one; however the two functions 
on which the metric depends can be moved both to dy 2 using the coordinates freedom. 
That is we can rename x 1 = r and put R(x l ) = r. The coordinate r is known as the 
areal radius, and is a geometrical invariant being the quotient of the area of space-like 
spheres by 47r. This will be referred to as the r-gauge. Sometimes the normal metric 
can be further written, somewhat redundantly, in the expanded form 

dy 2 = -E(r,t)dt 2 +2F(r : t)dtdr + G(r : t)dr 2 : y = -(EG + F 2 ) (4.19) 

especially if one wants to discuss certain limits, but we will always try to get metrics 
regular across the trapping horizons. If we wish, we can also locally write the line 
element in a double-null form which is conformally flat 

ds 2 = -2e f dx + dx~ + r 2 (x + ,x~)dn 2 (4.20) 
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with / = f(x + ,x ). For instance, starting with ea. (14. 131 we can introduce null coor- 
dinates via 

e f dx - =Ul- 2m ( r > V ) ) e 2* dv _ e * dr dx + = dv 



where / is an integrating factor. The normal section of ea. (|4.13|) then takes the form 
cq. (l4.20[) . One may use one spatial and one temporal direction on putting dx + dx~ — 
dt 2 — dr 2 , but there is no unique choice of such directions. The remaining coordinate 
freedom consists of conformal diffcomorphisms (in the language of two-dimensional 
metrics) 

x ± ^x ± (x ± ). 

The double-null form will be also referred to as the conformal gauge. 

Another form we would like to describe is the metric in the so called synchronous 
gauge. Let us consider the Schwarzschild space-time in coordinates (T, r, 0, </>) such 
that the line element can be expressed as 

dr 2 

ds 2 = -dT 2 + + (r g B) 2 dfl 2 , (4.21) 

B 



where r g — 2M is the usual gravitational radius, and 



B(T,r) 



2r g 



(4.22) 



We shall refer to these coordinates as the Lemaitre-Rylov gauge. This is indeed an in- 
teresting (time-dependent) gauge since, contrary for example to isotropic coordinates, 
(T, r) extend beyond the gravitational radius, r < r g . Now considering B(T, r) as an 
arbitrary function and replacing r g B with a general function R(r,T) we obtain the 
metric in synchronous gauge 

dr 2 

ds 2 = -dT 2 + — . + R 2 {T, r)dil 2 (4.23) 
B(T, r) 

in which the metric is diagonal but the areal radius is a function of r and T. 

The last form we would like to mention is the isotropic gauge in which the spatial 
part of the metric is conformally flat 

ds 2 = -A(t, p)dt 2 + B(p, t){dp 2 + p 2 dn 2 ) . (4.24) 

An important example of this latter form is the McVittie solution describing in author's 
mind a point mass in a FRW flat cosmology |165| . It reads as in eq. (|4. 241 with 
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When M = it reduces to a spatially flat FRW metric with scale factor a(t); when 
a(t) = 1 it reduces to the Schwarzschild metric with mass M in isotropic coordi- 
nates. This solution had a strong impact on the general problem of matching the 
Schwarzschild solution with cosmology, a problem faced also by Einstein and Dirac. 
Besides McVittie, it has been extensively studied by Nolan in a series of papers [166 . 
To put the metric in the general r-gauge form, the coordinate transformation 

transforms the metric in the so called Nolan gauge, in which it reads 

ds 2 =-(A s - H 2 {t)r 2 ) dt 2 - drdt+^-+ rW (4.27) 



2yfT s A s 

where H(t) = a(t)/a(t) is the Hubble parameter and, for example, in the charged 
four-dimensional case, A s = 1 — 2AI/r + e 2 /r 2 . Unlike the Schwarzschild case, A s = 0, 
or r = 2M in the neutral case or p = M/2a(t), is a curvature singularity rather than 
a global event horizon. In fact, it represents a Big Bang singularity very similar to 
r = in FRW models. When H = one recovers the Schwarzschild solution. For 
constant H, it reduces to the Schwarzschild-de Sitter solution in Painleve coordinates. 
As we shall see, the McVittie solution possesses in general black hole and cosmological 
trapping horizons, and the space-time is dynamical. 

Other coordinate systems which can be used to include dynamics in the Schwarzschild 
metric are discussed in [167] . We may note that ea. (|4.23l contains as a special case 
the FRW family of metrics. We stress that all forms discussed above are completely 
equivalent ways to describe the space-time structure of a spherically symmetric field, 
the use of any particular form being dictated only by computational convenience. We 
also note that in the warped form eq. (|4.17p the sphere of symmetry can be replaced by 
anyone among the two-dimensional manifolds with constant curvature which appear 
in the static topological black hole solutions embedded in Anti-de Sitter space. Only, 
the rotational isometry group must be replaced by the appropriate isometry group of 
the surface, which in the static case are torii or compact Riemann surfaces with higher 
genus. 

4.2.2 Trapping horizon 

To compute the trapping horizons we shall start from the Bardeen-Vaidya form eq. (|4.13l) . 
We can take the two null, future directed congruences normal to spheres of constant 
radius as 



£_ = [0,-e - *0,0] (4.28) 



where as always £+ ■ I- — — 1 . A simple computation gives the optical scalars 

U l-™)e*, 0_ = -^. (4.29) 



r \ r 
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Trapped or marginally trapped spheres have 6* + 6L > 0, therefore there are no such 
surfaces in the region r > 2M(r, v). We also see that < along the surface defined 
by 6+ — 0, which is therefore a trapping horizon of the future type. Its defining 
equation is r = 2M(r,v), which dchncs a line r = r(v) in normal space and therefore 
a hyper-surface with topology R 1 x S 2 in space-time. 

Consider now where the symbol cL denotes the directional derivative along 

l_. One obtains, using M' = d r M, 

r z \ r 

On the trapping horizon d-8 + = —(1 — 2M')/AM 2 will be negative if and only if 
M' < 1/2. As anticipated in Section [4.21 this is the condition which ensures the 
positivity of the surface gravity. Thus in this case r = 2M is a FOTH. The areal radius 
of the horizon will be denoted by th from now on. The signs of 9± are geometrical 
invariants, but their actual values are not because the null directions are defined up 
to an overall scale. An invariant combination is \ — 2g ± 8 + 8- or, using ea. (l4.17[) . 

X = l ij diR djR = gT d^R d v R . (4.30) 

X = is then the condition for a trapping horizon in Hayward's sense. For example in 
any coordinate system where one coordinate is the areal radius R = r, the condition 
is g rr = 0. The trapping horizon of McVittie solution can also be computed using 
ea. (|4.30p : using Nolan form it is a solution of the equation \J~A S = Hth , which in 
turn implies l — 2M/rn = H{t) 2 r 2 H . As M is a constant this a cubic algebraic equation 
with a priori more than one real root and in any case at most two positive roots. The 
situation is similar to de Sitter space, except that the horizon radius is here a function 
of time, rn(t). For positive H both horizons are of the future type, but as a rule one 
is outer (meaning that there is a black hole) while the larger root corresponds to an 
inner horizon, hence to a FITH, if and only if 



M „ H 

H 



H<r H -->0. (4.31) 



For M = 0, i.e. for homogeneous cosmology, only the inner cosmological horizon 
survives. This is the case, for example, in de Sitter space-time. We will make use of 
these results in Section [5j 



4.2.3 Misner-Sharp-Kodama energy 

One special feature of spherically symmetric space-times is the absence of gravitational 
radiation. This feature makes it possible the existence of a special, privileged notion 
of energy, the Misner-Sharp mass, which for spheres with areal radius r is the same 
as the Hawking mass j!68j (a general reference for energy in GR is the review [169 ). 
The energy may be defined by [148] (we recall that r is the areal radius) 

3 

E=~- rg+-d+rd-r = ^~^J (4-32) 
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and interpreted as the energy inside a sphere of radius r; so, by definition, a metric 
sphere is trapped if and only if E > r/2, marginal if and only if E = r/2 and untrapped 
if and only if E < r/2. Note that ea. (|4.32j) is a special case of the Hawking mass 

£(s) = ^)( 1 _^^- +M _ <fS ) (4 . 33) 

and is a geometrical invariant. Two very important properties of E were proved by 
Hayward 148J: in an asymptotically flat space-time, E coincides with the Bondi-Sachs 
scalar energy at null infinity, and with the Arnowitt-Deser-Misner (ADM) mass at 
spatial infinity. Using ea. (|4.13p the energy takes the implicit form 

g^d u rd v r = g rr = 1 - — (4.34) 
r 

which shows the relation with the Schwarzschild mass. Then E = M and the FOTH 
is at r — 2E. Many important properties of E can be displayed using also the Einstein 
gravitational field equations. We report two of them which seem remarkable, leaving 
the interested reader to the literature. Using double-null coordinates the variation of 
E as determined by Einstein's equations is 

d±E = 2Tre- f r 3 (T + -6±-T±±9 T ) (4.35) 

where we remember that e - -' = — . These field equations can also be written in 
the Bardeen-Vaidya form, using eq. (|4.131 where it is seen that E = M and 

d v M = 47rr 2 T; , d r M = -Aitr 2 T v v . (4.36) 

Thus in vacuo E is a constant. It can then be shown that the solution is locally isomet- 
ric to a Schwarzschild solution with energy E. This is an improvement of Birkhoff 's 
theorem. 

The second result is essentially the area law: if the NEC holds on a FOTH then 
E = rn/2 is non- decreasing along the horizon. 

The question arises naturally whether E is the charge associated to a conserved 
current. It was discovered by Kodama [158] that in spherical symmetry there is a 
vector field K^ 1 such that V M (G ^ V K V ) = 0; by Einstein equations it follows also the 
conservation equation 

V M {T^K V ) = (4.37) 
and a corresponding charge. If we define the two-dimensional Levi-Civita skew tensor 

_ a 13 

where n, are tangent vectors to constant radius spheres, then the Kodama vector 
may be defined by 

K" = eTdyT ■ (4.38) 

It can easily be seen that 

V^K^ = 0. (4.39) 

Let us give few examples: 



G2 



(a 
(b 
(c 
(d 

(e; 
(f 
( 

(h 



for static, non-dirty, black holes it is the Killing field; 



for dirty black holes, en. (|P3)l or (HOTI) . K = y/W/V d t : 
in Bardeen-Vaidya gauge K = e _ *9„; 

for the metrics ea. (|4.15D and ca. (|4.16[) . K = d tp or K = e~*<9j , respectively; 
for the metric eq.gJl]) X = (EF + G 2 )~ 1 l 2 d t ; 
in conformal gauge K = —g^ (d+rd- — 9_r9 + ); 
in synchronous gauge K = >/~B(R'd T — Rd r ); 
for FRW flat cosmology K = dt — Hrd r . 



In every case it follows that 

, 2E 

K 2 = 1 

r 

so that K is time-like, space-like or null if and only if r > 2E, r < 2E or r = 2E, 
respectively. Let us also define the current 

f = -T^Ky . (4.40) 

We have just seen that V M j M = 0, W^K^ = 0, and therefore there exist two conserved 
charges 

Qj = - f J^n^V (4.41) 



Q K = - K^dPV (4.42) 
Js 

where S is a space-like three-dimensional surface with fixed boundary at some constant 
r and future pointing time-like normal n^. Here "conserved" means independent on 
the choice of £. The charges as defined will be positive in regions where j 1 *, are 
both time-like. Using ea. (|4.35l) one can easily see that 

j = ^(8+Ed. - d.Ed+) . (4.43) 

To compute the charges we pass to a synchronous gauge by choosing coordinates (t, £) 
adapted to S, with d T normal and 8q tangent to E. In these coordinates we can always 
write the metric in the form 

ds 2 = -dr 2 + e x d( 2 + r 2 (r, Qdtt 2 . (4.44) 

From ea. (|4.38p and ea. (|4.40l) we obtain 

K" = e- x ' 2 (/ ) -r,0 i 0), ^=^^(£',-£,0,0) (4.45) 
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where f = dr/dr, r' = dr/d(. Following Hay ward, we shall say that the point r = 
is a regular centre if it is a boundary point of the normal space and E/r — >■ as the 
centre is approached. Otherwise it is a central singularity. We now assume that £ 
extends from a regular centre to some r > 0: from eq. (|4.45p and the metric we have 
the normal n M = Sq and 

p-A/2 

Therefore integrating over £ with the invariant measure dV — e x / 2 r 2 sin dd^dOdcj), we 
finally obtain 

A<rr 

QK = yr 3 , Qi=E. (4.46) 

We conclude that E is indeed the charge associated to a conserved current. The 
charge Qj is the definition of energy of Kodama. For static solutions the Kodama 
vector coincides with the Killing field, the generator of the time translation symmetry. 



4.2.4 The surface gravity 

The properties of the Kodama vector field discussed so far prompt for a natural defi- 
nition of the surface gravity of a trapping horizon. We have seen that K becomes null 
precisely on a trapping horizon and space-like within. The time-like integral curves of 
K are in general not contained within the horizon. Consider the quantity K^V ^K^y. 
it can be seen that on a trapping horizon it is proportional to K v 

K»V {V K A = -kK v . (4.47) 

The function k is, by definition, the horizon surface gravity of Hayward. For static 
black holes K is the Killing field so V \iK v is anti-symmetric and the definition reduces 
to the usual one. A formula to compute k efficiently was found by Hayward [139 . 
Working in double-null coordinates we have K + = —d+r, K- = cLr (see point (f) of 
the examples list) therefore 

d+K_ = -d-K+ = d-d+r = U y r (4.48) 

where O y = 2g^ V+V_ is the two-dimensional Klein-Gordon operator acting on 
scalars. Computing the left hand side of ea. (|4.47l) then gives the wanted formula 

K = ~n y r = \ r (g+-d-6 + + lg+-9 + 6^j , (4.49) 

where the last form is obtained by using 0± — 2r~ 1 d±r. That is, k is the "Box" of r; 
when evaluated on the trapping horizon, where 9+ = 0, it is the surface gravity and 
is positive if and only if d-9 + < 0, that is if the horizon is of outer type. Its invariant 
character is manifest. Let us give few examples with comparison to other definitions. 
For the Bardeen-Vaidya metric one obtains 

«=^(1-2M') (4.50) 
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which is also the Visser dynamical surface gravity as defined in [12] in a Painleve- 
Gullstrand frame. The first Hayward's definition (see the list in Section |4~TT) . invoked 
without appeal to un-affinity, would give instead 

1 



k = y/l - 2M' 

AM 

which is not even correct for the static Reissner-Nordstrom solution^. The definition 
of Fodor et al. gives 

k= — (l-2M') + $. 
For the dynamic Painleve-Gullstrand metric ea. (|4.15p with M = M(t p ,r), we obtain 

K=^ M (l-2M' + 2M) (4.51) 

while the Visser and Visser-Nielsen surface gravity for this kind of metrics would 
give the same formula without the time derivative term. Incidentally, this shows that 
the dynamical Bardeen-Vaidya metric is not diffeomorphic to a dynamical Painleve- 
Gullstrand metric ea. (l4.15p with the same mass function, as it would happen in the 
stationary case, and therefore represent physically different gravitational fields. For the 
dynamic Painleve-Gullstrand metric eq. ([4.16[) . with mass function m(t,r), we obtain 
instead 

k= — (1 - 2m' + 2me - *) (4.52) 
4m 

again different from the Visser-Nielsen surface gravity for the same mass function, 
which has no time derivative terms. One may also mention the "effective surface 
gravity", re e // = l/2rn discussed, for example, in |144j and |170j . 



4.3 Tunnelling from trapping horizons 

We now come to review what the tunnelling method has to say about dynamical, 
spherically symmetric black holes. We shall start by identifying the dynamical version 
of the tunnelling path which was displayed and discussed is Section [2J see figH] This 
will be accomplished by using a specific and convenient form of the metric, which we 
start to review a little more than already done. In absence of analytical techniques 
we shall make explicit use of the Hamilton-Jacobi equation, the null geodesic method 
being quite inconvenient to treat truly dynamical metrics (but see Clifton |164j V The 
covariance of the method will be stressed throughout. 

4.3.1 Metrics to be used 

We shall use for the time being the Bardeen-Vaidya (BV) metric eq. (|4.13|) . which we 
recall here 

dg 2 = _ e 2Hr,v) ( ! _ 2M(r,v) \ ^ + 2e Mr,v) dvdr + r 2 dQ 2 (4 53) 



For which M' = q 2 /2(2M 2 — q 2 + 2M^/M 2 - q 2 ). 
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The sphere of symmetry will not play any role here, though. To illustrate the covari- 
ance of the results we shall occasionally make use of the metric in Painleve-Gullstrand 
(PG) form, either eq. (|4.15p or eq. (|4.16|) . Few things about metric (I4.53[) will help with 
the understanding. The field equations read 

™ f =4»re*T- rl (4.54) 

and the stress tensor can be written as 

M M' 
T »» = V/.«V V « - V^rV^ . (4.55) 

If M only depends on v, it describes a null fluid obeying the dominant energy condition 
for M > 0. For the excreting black hole M < so the null energy condition will also be 
violated. We already know that r — 2M(r, v) is a trapping horizon, that is a FOTH, 
if and only if 2M' < 1, which we shall assume from now on. Putting r — 2M into the 
metric gives 

a 2$ ( 4e~*M \ 2 

so we conclude that the horizon is certainly time-like if M < 0; we know from general 
results that it will be space-like or null if M > 0, even if it is not evident from this 
expression. The Hayward surface gravity is 

«=^(1-2M') (4.56) 

and is positive under our assumptions. As we extensively explained, the Misner-Sharp 
mass, or energy for short, is the value E(v) taken by M(v, r) on the trapping horizon, 
i.e. 

E(v)=M{v,r H {v)) = T -^-. (4.57) 

Using ea. (l4.54p one can show that an observer at rest at r 3> rjj sees a quasi-static 
geometry with a luminosity L = —dE/dv. We shall not assume, initially, that L has 
the Hawking form L — HJ\fE~ 2 , with AT a constant proportional to the number of 
massless species radiated from the black hole. 

For sake of completeness, we note that a cosmological constant can be introduced 
via 

ds 2 = _ & mr,v) (l _ 2M M _ dv 2 + 2e *(r,v) dvdr + r 2 dn 2 ^ 

We shall admit that A > 0; if M is constant and $ = one recognises the Schwarzschild- 
de Sitter metric in advanced coordinates. Let us call C the expression in round brackets 
above; the null normals and their expansions are 



U = 



,-e*C,0,0 



9+ = ie*C (4.59) 
r 



I- = [0,-e-*0,0], 0_ = -^-^. (4.60) 
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Hence the horizons are located in correspondence of the roots of the equation C = 0; 
the tractable case is M = M(v). Then this becomes a cubic equation which for 
< 9AM 2 < 1, as is well known, admits precisely two real positive roots r c , r^, with 
r c > r& by definition. We see that + vanishes at both roots and 6L < 0, therefore 
the horizon spheres r — r&, r — r c , are marginally trapped surfaces of the future type 
which foliate a black hole and a cosmological trapping horizon, respectively. One can 
easily show that 

d-8+ = --d r C. 
r 

Computing the radial derivatives at both horizons we see that this is negative at rt 
and positive at r c . The cosmological horizon is therefore an example of a trapping 
horizon of inner type, the black hole horizon at r = Tf, remaining of the outer type. 

4.3.2 Rays tracing 

The most important features of a dynamical black hole of the kind discussed here are: 
(i) the existence of the irremovable space-like singularity at the origin r = of the 
coordinate system; (ii) the possible existence of a global event horizon (H in figHJ); 
and (iii) the time-like future trapping horizon (TH in figHJ). 




Figure 4: Evaporating black hole: the trapping horizon is at points where the radial 
outgoing rays have zero speed (dotted line). The event horizon % is repre- 
sented by a dashed line. Red lines represent curves of constant r. 
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We have seen that the BV form and the PG form are locally diffeomorphic to each 
other with the mass function transforming as a scalar field, so the causal structure 
of both solutions must be the same. A visual picture of the trapping horizon for an 
evaporating black hole is displayed in figH] The horizontal line represents the space- 
like singularity r — 0. Consider the radial outgoing null rays: in BV form they obey 
the differential equation 

• dr 1 a / 2M\ 

r= — = -e* 1 4.61 

dv 2 V r J y ' 

while ingoing rays simply are v = vq, for some constant vq. For comparison, in the 
PG form eq. (|4.15p the same equation takes the form 



r = ±1-A/ e* (4.62) 




where m(t, r) — M(v, r) and the plus (minus) sign referring to outgoing (ingoing) rays. 
Returning to the BV form, the event horizon is an outgoing null surface, 



If 2M\ 



r EH = 7j 1 e* < (4.63) 

so that teh < th = 2M . The acceleration close to the trapping horizon is 

r = -e*7TTT + ttt(1 - 2M'V ~ L + —(1 - 2M')f > . (4.64) 

2M 4M V ; r H 4M V 1 y ' 

We see that in both cases at the TH, which we recall is the line r = 2M(v, r) for BV or 
r = 2M(t, r) for PG, photons are only momentarily at rest, subsequently escaping on a 
dynamical time scale where k is the surface gravity ea. (l4.56|) . Therefore outgoing 
photons (massless particles) have r < before reaching the trapping horizon and r 
must decrease on going to the horizon. This can be seen from the lines of constant r in 
figdl Implicit in the Figure is also the idea that the Misner-Sharp mass has a non zero 
limit as v — > — oo. More complicated dynamics are possible in general depending on 
the mass function, including the possibility of horizon formation at finite times (some 
of these are discussed in Kodama's original paper |158j V 

The fact that f = at the trapping horizon and only there is one of the most 
important facts at the root of the quantum tunnelling phenomenon. Now we would like 
to consider the case where M(r, v) goes to a constant at very large negative advanced 
time and the trapping horizon extends to past infinity, since this is the region when 
the black hole does not radiate yet. In this case all radial null geodesies emerging from 
the trapping horizon seem to do so at v — —oo, because the speed vanishes at r — 2M 
(or 2M). For particles with non zero angular momentum the trapping horizon is not 
a surface of momentarily zero speed. For instance, in the equatorial plane the velocity 
is given by 

i / ■> i r \ 7- 

(4.65) 
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where J is the conserved angular momentum per unit mass. The essential point is that 
r < at r — 2M, so that for a rotating particle it will even be harder to escape quantum 
mechanically than to a non rotating one. We will see in the next section that the 
presence of angular momentum will not affect the main conclusion that the radiation 
originates close to the trapping horizon, not the global event horizon. Actually, the 
global event horizon cannot even be computed without knowledge of the dynamics near 
the evaporation end point. A practical definition was given by Bardeen [131] : one looks 
for photons that can only reach large r in a time comparable to the evaporation time 
~ E/L ^> E. Thus, the event horizon is located by the unaccelerated photons with 
f = 0. From eq. ([4.6ip and eq. (]4.64l) it follows then that 

f EH ~ -2L/r HK , r EH ~ 2E ( 1 - — ] ~ 2E(1 - 4L) (4.66) 

should be a solution to eq. ([4.631 to first order in the black hole luminosity L (which is 
very small for large black holes) and ignoring contributions due to $ and M 1 which 
are known to be negligible to first order in L [132] . This is the result of Bardeen- York. 

4.3.3 Tunnelling paths 

We propose now to compute the probability amplitude for a quantum particle to 
cross a trapping horizon. We do this by making use of the Hamilton-Jacobi method. 
Recall that, in the spherically symmetric dynamical case, it is possible to introduce 
the Kodama vector field K such that (K a G a p)'^ = 0, actually its defining property. 
Given the metric (|4.17[) . the Kodama vector components are 

K\x) = -JL=e ij djR , K 8 = = K* . (4.67) 

where e JJ is the numerical skew tensor with e 01 = 1. The Kodama vector gives 
a preferred flow of time and in this sense it generalises the flow of time given by 
the Killing vector in the static case (see [171] for a detailed study of the space-time 
foliation it determines). The conserved charges associated to K, eq. (]4.41l and eg. ([4.421) . 
are respectively the volume and the Misner-Sharp mass of space-time. We may also 
use the Kodama flow to define the invariant energy associated with a particle by means 
of the scalar quantity on the normal space 

w = -K i d i I, (4.68) 

where I denotes the classical action of the massless particle, which we assume to satisfy 
the reduced Hamilton-Jacobi equation 

Y"r),li),l 0. (4.69) 

Thus, for example, ui = —e~^d v I will be the particle energy on an extremal in a BV 
gauge; we may note that this gauge is only fixed up to diffeomorphisms v — > v = v (v), 
under which the field $(u,r) transforms as a conformal field, $(u,r) — > $ (v , r) = 
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<3?(i>, r)+log \dv /dv\. Therefore uj will not depend on the choice of advanced time which 
respect to the BV form of the metric. Similarly, in Painleve-Gullstrand gauge the in- 
variant energy is u = e^^d^I, and is gauge invariant under time re-parametrisation 
£->■? = t(t). 

To illustrate the method we shall work in s-wave and omit the angular dependence 
of /. We stress the importance to have at disposal an invariant definition of energy, 
eg. (14.681 certainly satisfies this requirement if the action is a scalar. In the following 
our aim will be to show that there is a precise invariant prescription to deal with the 
imaginary part of the action, in case there is one, which is valid for all solutions in all 
coordinate systems which are regular across the horizon. The task is to compute 

Imj^dl (4.70) 

along a tunnelling path, say 7. A priori there are two types of such tunnelling paths, 
which we shall call type-I and type-II, because a pair may form on either sides of the 
trapping horizon. A type-I path corresponds to a pair forming outside and therefore 
will be a backward null ray coming out from the future singularity at r — 0, crossing 
the trapping horizon and subsequently escaping to infinity, as illustrated in figlH The 
figure is the dynamical analogue of the corresponding figure in Schwarzschild space- 
time, except that it is smooth at r = 0, namely at the trapping horizon. A type-II 
path corresponds to a pair forming in the trapped region, so we may take a backward 
null ray crossing the event horizon from the singularity till some interior point, which 
subsequently escapes to infinity along a null ray crossing the trapping horizon at some 
other event. In some cases such a path (namely an outgoing light ray emerging from 
the trapping horizon) seems to emerge from the trapping horizon at v = —00 in ad- 
vanced time, therefore in this case, type-II paths are those with exactly one point, the 
starting point, lying on the horizon. Type-II paths are absent in static geometries, 
since there is no region in between the two horizons. The segment of type-II paths 
crossing the horizon outward seems nevertheless an allowed classical path so one may 
wonder whether they can contribute to an emission amplitude: does the pole disap- 
pear? We will see indeed that on such paths the radial momentum vanishes, basically 
for the reason that the speed is zero on the trapping horizon — and only there — so 
in fact they do not contribute semi-classically. 

We shall illustrate these results in some of the gauges we discussed above for spher- 
ically symmetric space-times, having two main purposes: to illustrate the intrinsic 
covariance of the result and to display how contributions to the tunnelling amplitude 
may or may not come from all terms in the differential form dl — dildx 1 . We also 
have the opportunity to test the formalism and to see how different can be the prosaic 
description of the same physical effect from the perspective of other coordinate sys- 
tems, a fact well known in special and general relativity. 
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The BV gauge 



Here the metric is 

ds 2 = -e^Cdv 2 + 2e*dvdr + r 2 dfl 2 (4.71) 

where M = (1 — C)r/2 is the Misner-Sharp mass discussed in Section T4.2.3I The 
Kodama vector and the invariant particle energy assume the simple expressions given 
by K — (e _ *,0) and u> = —e~®d v I, while the invariant surface gravity is just given 
by kh — d r Cnj2. As noted above <f> transforms as an ordinary Liouville field, i.e. 
<3> — > <£> + log |<9£>/<9i;|, under v — > v(v), making u) invariant under re-parametrisation of 
the advanced time coordinate. 

Type-II tunnelling paths are easily dispensed for: they cross the horizon along a 
null direction which is outward with zero speed at the horizon, so r = at and in a 
neighbourhood of the horizon, or dr = along the null direction; therefore on a small 
segment crossing the trapping horizon we have 

d r Idr = 0. 

On a real trajectory the energy is real and finite everywhere, which shows that the 
temporal u-contribution does not play any role, too, in what concerns us at the moment 
(the evaluation of Im/). Thus Im J = on type-II paths. 

For a type-I path the ray crosses the horizon along a null ingoing ray but in the 
backward direction; such a ray has v — vq (some constant) and Sv = exactly, therefore 
again the d v Idv term will vanish at the horizon and Sr > 0. From the Hamilton-Jacobi 
equation we see that, for outgoing modes of the kind we consider, 



C{d r I) = 2uj . (4.72) 



Thus one has 



Iml = Im J (d r Idr + d v Idv) — Im J dr — 



J 7 O r C \ H (r - r H - it) k h 

where ujh is u) evaluated on the horizon, the quantity C has been expanded around 
the horizon pole, that is 

C(v, r)=d r C Ar + ... (4.74) 

H 

and Feynman's ie-prescription has been implemented in order to deal with the simple 
pole, kh = d r C\ H /2 is the horizon surface gravity and coincides with our geometrical 
expectations. Unlike the stationary black holes, where uj is a constant of motion, here 
we get the local energy and temperature as measured near the horizon, which are con- 
nected to the quantities at infinity by the corresponding red-shift factors. We see that, 
in a BV coordinate system, the temporal integration does not give any contribution to 
the imaginary part of the action of particles tunnelling through the trapping horizon. 
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It will be very important to discuss the physical meaning of this result but we prefer 
to postpone this discussion to the end of the section, after we will have shown that the 
tunnelling amplitude is not an artefact due to a special coordinate system but holds 
as good as in any other. 



The conformal two-dimensional gauge 

A coordinate system where the temporal contribution to the action plays an essential 
role is the general diagonal form of a spherically symmetric metric, which reads 

ds 2 = e^*' r ) (-dt 2 + dr 2 ) + R 2 (t, r)dn 2 . (4.75) 

In this form, the normal metric is conformally related to the two-dimensional Minkowski 
space-time. The \ function simply reads 

x = e-f[-(d t R) 2 + (d r R) 2 ] , (4.76) 

If R(t, r) is a monotonically increasing function of r (the normal case) then this leads 
to the future trapped horizon condition 

(d t R) H + (d r R) H = . (4.77) 

If instead R{r, t) is decreasing as a function of r, then we have to take the opposite 
relative sign, because we expect that an evaporating horizon should have a negative 
time derivative, (d t R)H < 0. A computation of the expansion scalars will confirm this 
fact. The Kodama vector and the associated invariant energy are 

K = (d r R, -d t R, 0, 0) , (4.78) 

lo = e,-^ {-d r Rd t I + d t Rd r I) . (4.79) 
The dynamical surface gravity reads 

k h = (-d 2 R + d 2 r R) jt . (4.80) 



H 



Due to conformal invariance, the Hamilton-Jacobi equation is the same as in two- 
dimensional Minkowski space-time, namely using double null coordinates d+Id-I = 0, 
and for the outgoing particle we may take 

d+I = d t I + d r I = 0, (4.81) 

since the radial momentum p r — d r I > and dtl < for a real outgoing particle. 
Given that for a type-II path the null expansion condition leads to 6x~ = St — Sr = 0, 
we get dl = d + Idx + + d-Idx~ = on account of the Hamilton-Jacobi equation, and 
there is no imaginary part. On reflection, this result seems to violate our basic tenet 
that there is no amplitude because the photon speed vanishes at the horizon. In fact, 
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the coordinate speed r = ±1 everywhere in conformal gauge and it never vanishes. 
However, the speed of the wave front is 

R = d r Rr + d t R 

and this vanishes on H for an outgoing photon with r = 1, on account of the horizon 
equation ea. (|4.77[) . Thus, what really matters is that the speed of the areal radius 
should vanish at the horizon. 

But there are also the type-I tunnelling paths, for which the null expansion gives 
Sx + = St + Sr = 0. Therefore in this case 

I = (dr d r I + dt dj) = 2 j dr d r I . (4.82) 

J *y J -y 



Furthermore, due to (|4.79p and (|4.8ip . one has 



and we have a pole at the horizon. Making use of near horizon approximation along 
the null direction, from ((¥7771) and (|£5D|) . one has {d r R) H + (d t R)n = 0, St + Sr = 0, 
thus with the understanding that we shall put Sr = r — m , we get 

e-^ [d r R + 8 t R] = e-^ H {d 2 rr R - d 2 rt R - 8 2 tt R + 2 tr R) (r - r H ) + . . . 

= 2k h (r - r H ) + ■ ■ • (4.84) 

and the imaginary part follows. If R(t, r) is a decreasing function of r then we shall 
have to take p r ~ d r I < but again dtl < so the Hamilton- Jacobi equation implies 
now d^I = along the tunnelling path. Moreover since r is decreasing as R increases 
(i.e. toward infinity), on an outgoing (ingoing) null direction we have St + Sr = 
(St — Sr — 0) and the horizon condition becomes 



(d t R - d r R) 



= 

H 



instead of eq. (|4.77l It is now simple to see that a tunnelling type-II path gives no 
amplitude at all, but that a type-I does produce an imaginary part with the right 
magnitude and sign. In fact in this case we obtain first 

I = 1 

e-^'(d r R-d t R)' 

then, recalling <9_I = 0, 

dl = d+Idx+ = 2d r Idr = -2 — — ^ — Trrr: dr = ~ 2 — . 

e-^(d r R-d t R) n H (r-r H ) 

The integral of dl has to be done on the segment of the path crossing the horizon 
along decreasing r, say from ri to r\ with r\ < ra, and this is of course minus the 
usual integral in increasing order; reversing the orientation gives back the usual result 

Ira i,dl=—. 
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The synchronous gauge 

Another coordinate system where the coordinate speed does not vanish at the trapping 
horizon is described by the line element 

ds 2 = -dt 2 + —?—dr 2 + R 2 (r, t)dQ 2 = dy 2 + R 2 (r, t)dil 2 , (4.85) 
B(r,t) 

in which the metric is diagonal, but R is a function of r and t. In this case, one has 

X = -{d t R) 2 +B(d r R) 2 , (4.86) 
If R(r, t) is an increasing function of r, the horizon \h = should be defined by 

(d t R) H = -^B^{d r R) H , (4.87) 

in which we are assuming again a regular coordinate system on the horizon, namely 
that Bh and the partial derivatives are non- vanishing. For R decreasing with r one 
should change the sign to the right of eq. (14.871 The Kodama vector reads 

K = (VBd r R, -y/BdtR, 0, 0) , (4.88) 

and the invariant energy 

lo = y/B(d r Rd t I ~ d t Rd r I) . (4.89) 
The dynamical surface gravity is evaluated to be 

kh = \ (-2d 2 R H + 2B H d 2 ,R H + -^d t R H d t B H + d r R H d r B H ^ . (4.90) 

Making use of the horizon condition, we may rewrite it 

kh = \ (-2d 2 R H + 2B H d 2 R H - -^=d r R H d t B H + d r R H d r B H ^j . (4.91) 

In this case, the Hamilton-Jacobi equation reads simply 

- {dtlf + B{d r I) 2 = (4.92) 

and for an outgoing particle we have to choose the root with d r I > 0. Therefore, in a 
type-II path crossing the trapping horizon we have \/BnSt = 5r, in a neighbourhood of 
TH, because r = \f~B. As a consequence the outgoing temporal contribution will cancel 
the radial one. Type-II paths will then give no tunnelling amplitude. On reflection, 
this happens because although the coordinate speed is non vanishing at the trapping 
horizon, the areal velocity indicating the speed of the wave front 

R = d r Rr + d t R 
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does vanish on TH for an outgoing photon, which has r = vB, on account of ea. (|4.87l 
On a type-I path instead ^/B^At = —A?', and the two contributions add so that 

1 = 2 / drd r I . (4.93) 

The Hamilton-Jacobi equation and the expression for the invariant energy lead to 

d r I = . (4.94) 

Bd r R + VBd t R 

Making the expansion along the null curve, for which \fB~nAt + Ar = 0, in the near- 
horizon approximation, one gets 

Im I = 2 ■ Im / dr = — - , (4.95 

J 1 2n H {r -r H - ie) n H 

leading to the desired amplitude. But notice that in this gauge, the temporal contri- 
bution is essential in order to obtain the correct result. 

The r-gauge 

The fact that type-II paths never contributed to the tunnelling amplitude suggest 
that it is a general fact. Our last example will use a general metric in the r-gauge to 
confirm this. The normal metric here is non-diagonal, but as in BV gauge R — r. We 
have 

ds 2 = dy 2 + r 2 dn 2 , (4.96) 
where the reduced normal metric is now taken in the (redundant) form 

dy 2 = -E{r,t)dt 2 + 2F(r,t)dtdr + G(r,t)dr 2 , F^O. (4.97) 

The horizon is located at the zeroes of 

X (t,r) =y»d i Rd j R = y rr (t,r) = EQ E + - % = (4.98) 

i.e. at Eh = E{t,rjj) — 0, provided Fh ^ 0. The Misner-Sharp mass inside a sphere 
of radius r is 

and its value on the horizon is the black hole mass by definition, M = rjj/2. The 
other ingredient, the Kodama vector, reads 

K= ( . 1 ,0,0,0) , (4.100) 
WF 2 + EG J V ' 

and the invariant particle energy follows 

dj 



(4.101) 

VWTeg 
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The dynamical surface gravity is computed to be 

1 



1 

2*yr 



H 



2F 3 



E'F - -EG 



H 



(4.102) 



where an overdot and a prime denote differentiation with respect to t and r, respec- 
tively. Now the horizon will be future and of outer type provided F > and hh > 0. 
We expect the mass within the sphere to be a decreasing function of r; also, for an 
excreting black hole the time derivative M should be negative as well. From eq. (|4.99[) 
this will be true if E' > and E > 0, although these are only sufficient conditions. It 
follows that for F < 0, kh > always but the horizon is foliated by past marginally 
trapped surfaces (G > of course), while for F > it is future and outer if the 
inequality 

' EG S 
~2F 



E' 



> 



H 



H 



(4.103) 



is true, which we shall assume from now on. The outgoing null rays at the horizon 
satisfy either r — or, if F < 0, r = —2F/G. Only the first type will concern us, since 
photons from the past horizon should always be allowed classically. In fact, from the 
Kodama energy expression (|4.101[) and equation (|4.98p as well, the Hamilton- Jacobi 
equation reads 



X(d r l) 2 - 2- 



ujF 



-.drl-u'G = 0. 



. (4.104) 

VEG + F 2 V ' 

Solving this, one sees that d r I has no pole at H (which is past for F < 0) and therefore 

no imaginary part is found. Returning to our case with F > 0, on a type-II path r = 

at the horizon therefore, as before, the radial momentum vanishes and no imaginary 

term can come from J d r I dr. 

It remains to consider type-I paths. From the metric the null radial expansion across 

the segment crossing the horizon outward gives 8t = — {G/2F)\ij5r, now with F > 0, 

so we can expand the \ function along this ingoing null direction to first order in 

5r = r — r# . We obtain easily 

(T ' r (r-r H ) + --- = 2K H (r-r H ) + 0((r-r H ) 2 ) (4.105) 



X - xSt+x'Sr = ( x' - — 



H 



where = means evaluation on the horizon and eg. (|4. 1021 has been used. Also, dtl — 
—Fhoj from definition eq. (|4.101[) and the horizon condition E = 0. Hence we end up 
with I given by the sum of a real term and a possibly imaginary part coming from the 
horizon: 



/ {drd r I + dtd t I) = I dr 

J "V J *-y 



d r I + -Gh^r 



(4.106) 



I = {drd r I - 

What is remarkable is that in this gauge the temporal part is present but, being regular, 
it does not contribute to the imaginary part of the action. From the Hamilton- Jacobi 
equation eq. (|4.104[) we get, choosing the solution with positive radial momentum, 

ujF 



d r I = 



VEG + F 2 x 



(2 + 0(x)) • 



(4.107) 
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and we see that this has a pole at \ — 0. Making use of this equation, Feynman's 
prescription and ea. (|4.105[) . one has once more 



Iml = / ird r I = Im f ir " f 1 + + ° {x) - - ^4 108) 

We have shown in general that type-II paths will not contribute a tunnelling amplitude 
to future trapping horizons (F > 0), because Sr = to first order in \ an d we are 
assuming the energy to be real and finite. As to type-I paths the situation also is clear, 
they give a amplitude fully controlled by the value of the geometrical surface gravity 
kh- This example also shows the possibility of a quantum tunnelling through past 
trapping horizons of outer type, as for F < there is an amplitude for a photon to 
cross with negative radial momentum, i.e. for a photon to be absorbed (we are using 
the mental picture of a past horizon as one from which it is impossible to enter but 
very easy to escape) . The fact that there is no imaginary part on type-II paths either 
because the coordinate speed of a photon vanishes or because the temporal contribution 
cancels the radial one was the source of much confusion in the past literature. For 
instance, Belinski ventured to say that there is no Hawking radiation from black holes 
just from this fact |172) . since he apparently forgot about the existence of type-I paths. 

We end this section by noticing once more the special role of the trapping horizon. 
The imaginary part being produced on crossing it, it should be natural to think that 
this is the place where most of radiation forms. This has been confirmed in |164) . where 
it is shown that asymptotic observers register a radiation flux that starts increasing at 
a time they see the collapsing shell crosses the surface r = 2M of the model used (a 
special case of Painleve-Gullstrand model with a time dependent mass, see ea. (|4.15[) . 
The stress-tensor for Vaidya space-time has been considered in |1731I174| and should 
be consistent with this result. 

A note — In the previous computations various choices of signs have been applied 
in such a way that it may seem they were chosen somewhat ad hoc in order to get the 
wanted result. This is not so. Once the future sheet of the trapping horizon has been 
chosen, and the sign of the Kodama vector determined so that it is future directed, no 
other sign uncertainties will occur for either outgoing or ingoing particles if one uses 
consistently the Hamilton-Jacobi equation for outgoing particles on every segment of 
the tunnelling path; furthermore, questions such as the right form of the equations 
determining the horizon, the orientation of the space-time coordinates and the sign of 
the metric coefficients are carefully established once for all. 

On the other hand, if there exists a past sheet in the trapping horizon then using the 
tunnelling picture we may as well compute the action along an inward directed curve at 
the horizon (the ambiguity inherent in this and analogous terms is easily resolved if the 
manifold is asymptotically flat). Then there will be again a non- vanishing imaginary 
part, but we can interpret it as a small absorption probability, as we did for the static 
case in relation to white holes. 
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5 Cosmological Horizons, Decays and Naked 
Singularities 

In this section we review the formalism of horizon tunnelling as applied to cosmo- 
logical horizons, decay of unstable particles that in the absence of gravity would be 
otherwise stable and radiation from naked singularities. Especially the subject of 
horizon tunnelling in cosmology attracted much interest recently, leading to a related 
stream of papers. A sample of these articles which seemed relevant to us can be found 
in [175H177] . although we rely mainly on the papers [69, 178 and the references cited 
therein. We shall follow mainly the Hamilton- Jacobi version of the tunnelling picture, 
although the null geodesic method can be still applied. Results in this direction ap- 
peared recently in 179 , where the Parikh-Wilczek and Hamilton- Jacobi methods are 
compared and showed to agree. An early study of the evolution of evaporating black 
holes in inflationary cosmology is in |180| . 



5.1 The FRW space-time 

Consider a generic FRW space-time, namely one with constant curvature spatial sec- 
tions. Its line element can be written as 

ds 2 = -dt 2 + dr 2 + [a(t)r} 2 dft 2 (5.1) 

1 — kr z 

where r is measured in units of the curvature radius and, as usual, k = 0, — 1, +1 labels 
flat, open and closed three-geometries, respectively. In this gauge, the normal reduced 
metric is diagonal and the horizon is implicitly given by the equation x = 0, where 

X (f,r) =l-[o(t)r 

namely 

R H :=a(t)r H = 1 =, with H (t) = ^9 , (5.3) 

provided the space-time energy density g(t) is positive. As always it corresponds 
to the vanishing of the wave front velocity whose radius is R = a(t)r. The surface 
R[i(t) coincides also with the Hubble radius as defined by astronomers for vanishing 
curvature, but we shall call it Hubble radius in any case. The important case k = 
deserves special attention. The horizon is the surface rd(t) — 1; substituting into the 
metric reveals the causal character of the horizon: it will be null if and only if either 
p = g/3 (radiation) or g — —p (vacuum energy), g and p being the energy density 
and pressure, respectively. It will be time-like or space-like according to whether 
l(f? + 3jj)/2g| is less or greater than one. In the former case — g < p < g/3, which 
covers almost the totality of cosmological models, while the space-like character can 
only be achieved if p < —g (dubbed phantom energy) or p > g/3, which includes 



H 2 (t) 



a 2 (t) 



(5.2) 
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stiff matter. These results actually hold for non zero curvature too. In Section 14.11 
we defined future outer trapping horizons. The situation in flat cosmology is easily 
spelled out: the horizon is a surface with 6>_ = 0, 9 + > that for — g < p < g/3 
satisfies also d+0- > 0. In the terminology of Section 14.11 it can be classified as a 
past inner trapping horizon (abbr. PITH), because all metric spheres at larger radii 
are trapped with 0± > 0. It means that given a comoving observer, a spherical 
light beam approaching him from cosmic distance will have increasing area due to 
cosmological expansion. One has to be careful here that the areal radius is R = a(t)r, 
not r, so the area of the wave front is AirR 2 and the horizon can be written as the 
condition HRh = 1. Examples of a FITH are provided by the de Sitter horizon and 
the collapsing Vaidya solution in de Sitter space, ea. (|4.58l 
The dynamical surface gravity reads 

k h = \u y [a{t)r] = - + iff (t) + R H (t) < , (5.4) 

and the minus sign refers to the fact the Hubble horizon is of the inner type. Similarly, 
the Kodama vector is 

K = a/I - kr 2 (d t - rH{t)d r ) (5.5) 
so that the invariant Kodama energy of a particle is equal to 



u = y/l- kr 2 (-d t I + rH(t)d r I) = \J \ - kr 2 Co . (5.6) 

Notice that K is space-like for ra > (H 2 + k/a 2 )^ 1 / 2 , i.e. beyond the horizon, so that 
particles will eventually tunnel from outside to the inner region, r < r#. 

The next ingredient will be the reduced Hamilton-Jacobi equation for a relativistic 
particle with mass parameter m, which reads 

~ (^ J ) 2 + (1 ~ 2 ^ 2) (drl) 2 + m 2 - . (5.7) 
a 2 (t) 

Solving eq. (|5.6l and eq. (|5.7p for d r I and dtl we obtain 

d t I = -Q+ rHd r I (5.8) 



aHCo(ar) ± aJu 2 - m 2 + m 2 (H 2 + A) (ar) 2 

drl — jT\ , (5.9) 

1 - (H 2 + £) (ar) 2 1 ' 

with the sign chosen according to which direction we think the particle is propagat- 
ing. The effective mass here defines two important and complementary energy scales: 
if one is interested in the horizon tunnelling then only the pole matters (since the 
denominator vanishes), and we may neglect to all the extents the mass parameter 
setting m = (since its coefficient vanishes on the horizon). On the opposite, in in- 
vestigating other effects in the bulk away from the horizon, such as the decay rate of 
composite particles, the role of the effective mass becomes relevant as the energy of 
the particle can be smaller than the energy scale settled by m, and the square root 
can possibly acquire a branch cut singularity, eg. (15.91 will be the starting point of our 
considerations since it embraces all semi-classical quantum effects we are interested in. 
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5.2 Horizon tunnelling 

Using ea. (|5.9l we may derive, following [55], the cosmic horizon tunnelling rate. To this 
aim, as we have anticipated, the energy scale is such that near the horizon, we may 
neglect the particle's mass, and note that radially moving massless particles follow a 
null direction. The horizon region being a region with strong gravity we expect indeed 
some "particle creation" . To a pair created near the horizon in the inner region r < m 
corresponds a type-I path, one segment crossing the horizon backward in cosmological 
time. To a pair created outside the horizon corresponds a type-II path, the segment 
crossing the horizon inward being classically allowed. As in Section 14.3.31 only type-I 
paths will contribute to horizon tunnelling. 

Take then a null radial direction crossing the horizon to the inner region: we have 

8t=^B=6r. (5.10) 
Vl — kr z 

The action for particles coming out of the horizon towards the inner region is to be 
integrated on a type-I path: then we must choose the solution with negative radial 
momentum, d r I < and we obtain 

I = J dtd t I + J drd r I = 2 J drd r I , (5.11) 

upon solving the Hamilton-Jacobi equation eq. (|5.7[) with zero mass and using (|5.10p . 
Note that the time derivative of the action contributes to the total imaginary part 
which, as we have seen, is a general feature of the tunnelling method. For d r I we 
use now Equation ea. (|5.9p . which exhibits a pole at the vanishing of the function 
F(r,t) := 1 — (a 2 H 2 + k)r 2 , defining the horizon position. Expanding F(r,t) again 
along the null direction ea. (|5.10p we obtain 

F(r, t) « +An H a{t){r - r H ) + ■ ■ ■ , (5.12) 

where kh given in eq. (|5.4I represents the dynamical surface gravity associated with 
the horizon. In order to deal with the simple pole in the integrand, we implement 
Feynman's it - prescription. In the final result, beside a real (irrelevant) contribution, 
we obtain the following imaginary part 

7TWH 

Iml= -. (5.13) 

This imaginary part is usually interpreted as arising because of a non- vanishing tun- 
nelling probability rate of massless particles across the cosmological horizon, 

r - exp(-2Im7) - e"^'" 11 . (5.14) 

Notice that, since Kh < and luh > for physical particles, (|5.13l) is positive def- 
inite. Due to the invariant character of the quantities involved, we may interpret 
the scalar T = —kh/2tt as the dynamical temperature parameter associated to FRW 
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space-times. In particular, this gives naturally a positive temperature for de Sitter 
space-time, a long debated question years ago, usually resolved by changing the sign 
of the horizon's energy. It should be noted that in literature, the dynamical tempera- 
ture is usually given in the form T = H/2tt (exceptions are the papers |181j b Of course 
this is the expected result for de Sitter space in inflationary coordinates, but it ceases 
to be correct in any other coordinate system. In this regard, the H and k terms are 
crucial in order to get an invariant temperature. The horizon's temperature and the 
ensuing heating of matter was foreseen several years ago in the interesting paper |182) . 



5.3 Decay rate of unstable particles 

We consider the decay rate of a composite particle in a regime where the energy of the 
decaying particle is lower than the proper mass m of the decayed products. A crucial 
point is to identify the energy of the particle before the decay with the Kodama energy. 
We also denote m the effective mass parameter of one of the decay products (recall 
it may contain a curvature terms). The relevant contribution to the action comes 
from the radial momentum given by equation (|5.9[) . If we introduce the instantaneous 
radius ro by 

[a(t)r ] 2 = R 2 := (l - ^] R% , (5.15) 

where Rh is the horizon radius given by Equation ea. (|5.3[) . then the classically forbid- 
den region is < r < ro. From eq. (|5.9l) . we see that for the unstable particle sitting 
at rest at the origin of the comoving coordinates, one has an imaginary part of the 
action as soon as the decay product is tunnelling into this region to escape beyond ro, 

f Ro \l R 2 - R 2 
lmI = mR H / dR V p2 p2 . (5.16) 
Jo Rh ~ 

The integral can be computed exactly at fixed t, and the result is 

Im/ = | J2 ff (m-uj) > 0, (5.17) 

leading to a rate which, assuming a two-particle decay, takes the form 

T = T e- 2nRH( - m - u \ (5.18) 

where Tq is an unknown pre-factor depending on the coupling constant of the inter- 
action responsible of the decay (for instance, for a A</> 3 interaction one should have 
r ~ A 2 ). Of course, each newly produced particle will itself decay, leading possibly to 
the instability mechanism first discussed by Myhrvold |184j in de Sitter space. Since 
the tunnelling process locally conserves energy one should put lj = m/2, so that the 
tunnelled particle will emerge in the classical region at r = ro with vanishing momen- 
tum. Furthermore, the result is again invariant against coordinate changes, since both 
oj and Rr are invariantly defined quantities. 
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A particularly interesting case is represented by de Sitter space. Taking for example 
the line element in the static patch 

ds 2 = -(1 - H 2 r 2 )dt 2 + dr " + rW , (5.19) 

(1 - H^r 2 ) 

for the imaginary part (15.18[) we obtain 

lmI=-^-(m-u,) (5.20) 

a result actually independent of the coordinate system in use. Putting u) = m/2, the 
above result has been obtained by Volovik |183) using the so-called "fluid" static form 
of de Sitter space 

ds 2 = -dt 2 + {dR - H Rdt) 2 + R 2 dfl 2 . (5.21) 
in agreement with the asymptotic approximation of the exact result due to [185 . 



5.4 Particle creation by black holes singularities 

One may also use the tunnelling formalism to investigate whether particle creation in 
the bulk of space-time is possible due to the presence of space-time singularities, for 
example due to static black holes. With regard to this, we consider the exterior region 
of a spherically symmetric static black hole space-time and repeat the same argument. 
Quite generally, we can write the line element as 

ds 2 = ~e 2 ^C(r)dt 2 + C-^rjdr 2 + rW. (5.22) 

From the Hamilton-Jacobi equation, the radial momentum turns out to be 



If we are interested in particle creation we should set u> = 0: in fact, according to the 
interpretation of the Kodama energy we gave before, this approximation simulates the 
vacuum condition. Then 

r 2 1 

drd r I = m I dr (5.24) 



ri 



where the integration is performed in every interval (ri, r2) in which C(r) > 0. eq. (|5.24l 
shows that, under very general conditions, in static black hole space-times there could 
be a production rate whenever a region where C(r) > exists. 

As a first example, let us analyse the Schwarzschild black hole. For the exterior 
(static) solution, one has C(r) = 1 — 2M/r > and <I>(r) = 0, thus the imaginary 
part diverges since the integral has an infinite range. We conclude that the space-like 
singularity does not create particles in the semi-classical regime. In the interior the 
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Kodama vector is space-like, thus no energy can be introduced. A similar conclusion 
has been obtained also for the Big Bang cosmic singularity, the only scale factor leading 
to particle emission being a(t) ~ £ -1 . This is like a big rip in the past. 

The situation is different when a naked singularity is present. Consider a neutral 
particle in the Reissner-Nordstrom solution with mass M and charge Q > (for 
dcfiniteness) given by the spherically symmetric line element 

ds 2 = - {r ~ r - )( ; r ~ r+} dt 2 + - £ dr 2 + rW . (5.25) 

r i (r — r_)(r — r + ) 



Here r± = M ± y M 2 — Q 2 are the horizon radii, connected to the black hole mass 
and charge by the relations 

M= r++r_ Q = ^__ (52g) 

The Kodama energy coincides with the usual Killing energy and 

C{r) = ^-r-){r-r + ) ^ 

The metric function C(r) is negative in between the two horizons, where the Kodama 
vector is space-like, so there the action is real. On the other hand, it is positive within 
the outer communication domain, r > r + , and also within the region contained in the 
inner Cauchy horizon, that is < r < r_. Thus, because of eq. (|5.24[) and assuming 
the particles come created in pairs, we obtain 

T r f~ r , mM , {M + Q\ , 

Imi = —m / — , dr — mq In — . 5.28 

Jo ^/(r--r)(r + -r) 2 \M - Q ) V > 

Modulo the pre-factor over which we have nothing to say, with the usual interpretation 
there is a probability 

/ M ~ n\ mM 

r ~exp(-2ImJ)=^^ T |J e~ 2 ™Q. (5.29) 

Pleasantly, eq. (|5.29p vanishes in the extremal limit M — Q. Being computed for 
particles with zero energy, we can interpret this as an effect of particle creation by the 
strong gravitational field near the singularity. Since the electric field is of order Q/r 2 
near r = 0, there should be also a strong Schwinger's effect. In that case one should 
write the Hamilton-Jacobi equation for charged particles. 

The processes just discussed should bear a bit on the question of the stability of 
the Cauchy horizons. Due to infinite blue-shift of perturbations coming in from the 
asymptotically flat exterior regions both sheets of the Cauchy horizon (Hp and Hp 
in fig'E]) are believed to be classically unstable. Of course, if the naked singularity 
is formed from collapse of charged matter, one asymptotically flat region (say, the 
left one) disappears. Taking into account particle creation, it can be easily seen that 
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escaping particles will reach the future portion Hp of the inner horizon with infinite 
blue-shift, or infinite Kodama energy, as measured by an observer on a Kodama tra- 
jectory. Hence, the future sheet will probably become unstable, this time by quantum 
effects. On the other hand, the particle reaching the singularity will do so with infi- 
nite red-shift, that is with zero energy, giving a negligible back-reaction and substan- 
tially not changing the nature of the singularity (its time-like character, for example). 

Thus, the only process poten- 
tially capable of modifying the 
singularity would be its conjec- 
tured screening by absorption 
of charged particles. Of course 
the mass inflation phenomenon 
is an issue here but, in the con- 
text of perfect spherical sym- 
metry, exact solutions of the 
Einstein-Maxwell system exist 
having both event and Cauchy 
horizons and describing collapse 
of spherical shells. For these 
what we said should apply. The 
Figure 5: A pair of particles is created in the com- present formalism also predicts 
pact region inside the inner horizon. the absence of particle creation 

in the region in between the 
two horizons, despite the metric 
there is dynamical. The possibility that the naked singularity itself radiates away its 
mass will be taken up in the next section. 



5.5 Naked singularities 

A general reference on the physics of naked singularities is [186] , Particle emission 
from naked singularities in higher dimensions has been studied also in 1187). Sticking 
to two-dimensional models for simplicity, consider the following metric f 88 

ds 2 = a~ x dx + dx~ , a = X 2 x + x" - a(x + - x+)8(x + - x+) (5.30) 

where A is related to the cosmological constant by A = — 4A 2 . 

This metric arises as a solution of two-dimensional dilaton gravity coupled to a 
bosonic field with stress tensor T ++ = 2a5(x + — Xq ), describing a shock wave. A 
look at figH] reveals that a — is a naked singularity partly to the future of a flat 
space region, usually named the linear dilaton vacuum. The heavy arrow represents 
the history of the shock wave responsible for the existence of the time-like singularity. 
The Hamilton-Jacobi equation implies either d+I = or d-I = 0, I being the action. 
To find the ingoing flux we integrate along x + till we encounter the naked singularity, 
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using d-I = 0, so that 



'= [dx+d + I= [u>^ = j— +l r (5.31) 

J J 2a J 2{\ 2 X - -a){x+ + ax+/C-ie) 1 ' 

where C = C(x~) := (\ 2 x~ — a) and to = 2ad+I is the familiar Kodama's energy. 
The imaginary part immediately follows on using (x — ie)^ 1 = P- +in5(x), giving the 
absorption probability as a function of retarded time 

r(w) - r e- 2 Im/ = r e~™ /c{x ~ ) , (5.32) 

T being some pre-factor of order one. 



null 
infinity 




The flux is computed by integrating 
the probability over the coordinate fre- 
quency (that is, the variable conjugated 
to the coordinate time) <2> = to/a, with 
the density of states measure du>/2n, 
which gives 



'-++ 



2 



af 



27T 3 CT 2 



(5.33) 



shock wave 



Similarly, in order to find the outgoing 
flux we integrate along x~ starting from 
the naked singularity, this time using 
d + I = 0. A similar calculation gives 



ImJ 



Figure 6: The naked singularity formed by 
the shock wave. 



TTUJ 

2X 2 x+ 



(5.34) 



then, integrating the probability over the 
coordinate frequency, we obtain 



\\x+) 2 

The outgoing flux is 2(T ++ — T ). The conservation equations 

ct<9 + T__ + <9_((7T + _) = 0, <jd-T ++ + a+(crT + _) = 



(5.35) 



(5.36) 



determine the components only up to arbitrary functions B(x~) and A(x + ), something 
which corresponds to the freedom in the choice of a vacuum. For instance, requiring the 
fluxes to vanish in the linear dilaton vacuum fixes them uniquely. As it is well known, 

T-\ is given by the conformal anomaly: T = 4aT + - = R/24n (for one bosonic degree 

of freedom). Matching to the anomaly gives the pre-factor r = 7r 2 /24, of order one 
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indeed. These results agree with the one-loop calculation to be found in |188) . Note 
that the stress tensor diverges while approaching the singularity, indicating that its 
resolution will not be possible within classical gravity but requires quantum gravity 
instead [15511190] . 

We return now to the Reissner-Nordstrom solution. Could it be that the naked 
singularity emits particles? In the four-dimensional case one easily sees that the ac- 
tion has no imaginary part along null trajectories either ending or beginning at the 
singularity. Formally this is because the Kodama energy coincides with the Killing 
energy in such a static manifold and there is no infinite red-shift from the singularity 
to infinity. Even considering the metric as a genuinely two-dimensional solution, this 
would lead to an integral for / 



I = 



(r-r+)(r-r-) ^ + 



where x ± = t ± r*, with 

,,2 



In 



r + — r_ \ r + J r + — r_ 
But close to the singularity 



r 2 = 



3r + r 



+'- / + --\ 
2~ -s ) 



2/3 



(5.37) 



(5.38) 



(5.39) 



not a simple pole, rather, an integrable singularity. It is fair to say that the Reissner- 
Nordstrom naked singularity will not emit particles in this approximation. This seems 
to be coherent with quantum field theoretical results. With the customary u — x~ and 
v = x + , there is a map u — > v = G(u) which gives the family of ingoing null geodesies, 
characterised by constant values of v, which after reflection in the origin emerge as the 
the family of outgoing null geodesies with constant u. According to [191j . the radiated 
s-wave power of a minimally coupled scalar field is given in terms of the map G(u) by 
the Schwarzian derivative 



w = 



1 

247 



3 fGT_ 
2\G> 



G'" 
~G T 



(5.40) 



The (u, v) section of the Reissner-Nordstrom metric is conformally flat, hence the 
above map is trivial (i.e. linear) and W = 0. 



6 Conclusions 

We have completed our journey into the realm of quantum gravitational effects that can 
be understood as a result of tunnelling processes. In this journey, we visited black hole 
radiation, de Sitter and cosmological radiation, decays and naked singularities, but we 
have not included tunnelling between different gravitational field configurations, nor 
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topology changing processes, nor sonic black holes or other analogue objects. Thus 
we may generically qualify the subject as covering semi-classical effects in a given 
background field. From this point of view it seemed remarkable that it could be 
extended to cover back-reaction effects to first order in the ratio E/M, where E is the 
energy scale of the tunnelling process and M the black hole mass or some generalisation 
thereof. We would like to conclude by summarising the achievements, the drawbacks 
and the prospects. 

Achievements. We list the results which seem to us more robust and have a chance 
to last longer. And also what is to be expected, reasonably, by pursuing this approach 
still further. Results that were obtained in Section [5] and Section [3] comprise: 

1. The ability of the method to include first order back-reaction effects in order to 
comply with energy conservation and the predicted correlations among successive 
emission events is one of the key results. 

2. In the stationary case, it is possible to provide a foundation of tunnelling methods 
using analytic continuation of the classical action alone, avoiding in this way the 
seemingly impossible task of crossing the horizon. The analytic continuation 
could be extended to cover white holes and to associate with them a well defined 
geometric temperature, a point not always appreciated in literature. 

3. In the same vein, one can prove the complete equivalence of the null geodesic 
method with the Hamilton- Jacobi method for stationary fields. We also gave 
a "covariantised" version of the null geodesic method, which allows to avoid 
questions such as covariance or canonical invariance. 

4. Many special cases have been considered, including rotating black holes, Godcl 
space-times, higher dimensional solutions, de Sitter and anti-de Sitter spaces, 
Taub and Taub-NUT solutions; also, we considered some supergravity solutions 
and fermion tunnelling, including gravitinos. The temperature of these solutions 
were computed and showed to agree with those obtained by geometric meth- 
ods, confirming the universal character and kinematics of the Hawking effect. It 
should also be noticed that the Hawking temperature computed from the tun- 
nelling method does not receive higher order corrections in h, contrary to past 
suggestions. 

In Section [4J using a local notion of horizon and dynamical surface gravity for time 
varying spherically symmetric black holes, the main results comprise: 

1. The application of the tunnelling method to non stationary, spherically symmet- 
ric black holes and the fact that it naturally selects a trapping horizon of the 
future non degenerate type, either inner or outer. It can also be applied to past 
horizons under the same condition and by interchanging emission with absorp- 
tion. In the limiting situation where the surface gravity vanishes it is consistent 
with the idea that extremal black holes have zero temperature. 
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2. One of the principal consequences of the tunnelling method is the observation 
that the radiation seems to originate near the local trapping horizon, not the 
global event horizon. 

3. Using the Kodama-Hayward theory of spherically symmetric space-times, and 
noting that one can possibly associate an observable temperature to a black 
hole only if it is a coordinate scalar, we showed that the tunnelling probability 
depends on the ratio, 2ttE/k, of the Kodama invariant charge taken as energy 
and the geometrical surface gravity as defined by Hayward. 

4. One can include non interacting fermions in the formalism, as demonstrated in 
several papers. Also noticeable is the fact that any mass term gets strongly 
suppressed by the horizon pole, which dominates the rate well over mass contri- 
butions. 

5. The time-like nature of the trapping horizon of an excreting black hole is con- 
sistent with the physical interpretation of the formalism. It is true that it can 
be crossed either ways, but it remains the fact that there is a class of paths for 
which an imaginary part exists, and these paths precisely correspond in Feynman 
diagram language to particle creation. 

In Section [S]we addressed other themes - cosmology, decays and naked singularities. 
Since FRW spaces are spherically symmetric, the general theory outlined in Section 
14.21 applies. We showed that, with obvious modifications, the tunnelling picture works 
equally well. In particular the type-II paths never contribute, while type-I paths do. 
We may consider these findings for cosmological horizons as a generalisation of the 
well known facts discovered for de Sitter space-time, in complete analogy with the 
parallel extension of the theory from stationary to dynamical black holes. We also 
considered the decay of particles which in the absence of gravity would not occur, and 
find that the Hamilton-Jacobi equation can manage them through the presence of a 
branch point singularity in the radial momentum, that would not be present in absence 
of an external field. From this point of view the tunnelling method has universally 
valid features. Again, its value is to be found where it provides approximate results 
when exact calculations are in general impossible, at the same time being consistent 
with the exact calculations when available. Such is the case of de Sitter space. We 
showed also that the presence of strong fields in a certain space-time region, even in 
the absence of horizons, can result in particle production by means of a tunnelling 
process "from nothing" . It can be considered as the gravitational analogue of the 
Schwinger effect in a strong external electric field. Finally, it is of interest that the 
tunnelling method can also be applied to naked singularities, although in this case 
there is nothing to tunnel through, the singularity being the boundary of space-time. 
We showed nevertheless that the emerging picture is quite in agreement with quantum 
field theory in two-dimensional models, making one confident of its general validity. 

As regard the drawbacks, in the stationary case one can perhaps complain that 
the method gives only approximate results. But for black holes with really intricate 
metrics exact calculations are illusory even in stationarity, so we feel that this is not 
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a serious drawback. Even if there is little doubt that the method is correct, still it 
remains of a hypothetical nature. However, we think that also in case of a failure due 
to some internal inconsistency, that should be nonetheless an important message. In 
fact, tunnelling processes are predictions of quantum theory, so a failure of them in 
gravity theory would certainly signal something interesting. Despite all the desirable 
properties discussed above, the tunnelling method is and remains an essentially semi- 
classical procedure carrying along with it all the limitations inherent to its nature. 
In particular, it applies only to free particles while in principle quantum field theory 
methods certainly have a wider scope. 

Prospects. The tunnelling method is not a closed subject. For instance, in striking 
contrast with the stationary case, an important missing point of the dynamical case is 
the absence of an extension of Kodama-Hayward's theory to dynamical axis-symmetric 
black holes, to be used as a tool for studying tunnelling. We mention here the papers 
[19211193] where a tentative theory is developed. We feel that this is perhaps the most 
important missing point. The first back-reaction corrections are also of great interest, 
as they are not so well developed in dynamical situations. Processes where small black 
holes are emitted by large ones should also be within the range of problems where the 
tunnelling picture could work, and that also is an important "to do" . 

In cosmology, the physical interpretation of Hayward's dynamical surface gravity 
as a temperature is not as clear as in de Sitter or black hole case, but it is strongly 
favoured by the tunnelling picture. In particular, an Unruh-DeWitt detector moving 
in de Sitter or black hole space-time always thermalises at the corresponding Gibbons- 
Hawking temperature even if it is not on a geodesic path, while for cosmology there is 
not an analogous result. One may conjecture that the cosmological temperature given 
by eq. (|5.4I is the one that a point-like detector registers when moving on a Kodama 
trajectory, whose four-velocity is the Kodama vector field. In general, such a trajectory 
is accelerated relative to geodesic observers so one expects a mixed contribution from 
the horizon and the acceleration together. These expectations were actually confirmed 
years ago in de Sitter and anti-dc Sitter space and resulted in beautiful formulas 
[19411195] . It would be very interesting to apply tunnelling methods also to Unruh's 
like detectors (but see |51) ) . Finally, two more intriguing prospects are the application 
to field theories where fields, rather than particles, tunnel across forbidden regions; 
and to sonic black holes [201] . or dielectric black hole analogues |202[ . since both can 
be studied in laboratory experiments. But these arguments would easily cover another 
review. 
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